THE NONLINEAR STABILITY OF THE TRIVIAL SOLUTION TO THE 
MAXWELL-BORN-INFELD SYSTEM 



JARED SPECK 

Abstract. In this article, we use an electromagnetic gauge-free framework to establish the existence of small-data 
global solutions to the Maxwell-Born-Infeld (MBI) system on the Minkowski space background in 1 + 3 dimensions. 
Because the nonlinearities in the system satisfy a version of the null condition, we are also able to show that these 
solutions decay at exactly the same rates as solutions to the linear Maxwell-Maxwell system. In addition, we show 
that on any Lorentzian manifold, the MBI system is hyperbolic in the interior of the field-strength regime in which 
its Lagrangian is real- valued. 
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1. Introduction 

The Maxwell-Born-Infeld (MBI) system is a nonlinear model of classical electromagnetism that was introduced 
in the 1930's by Born and Infeld jBI34j , with [Bor33j a precursor by Born. In this article, we study the source- 



free (i.e. the right-hand sides of (1.0.1a) - (1.0.1b) are 0) MBI system in the fixed spacetim^j] (M, g). We will 
assume throughout the article that (M,g) is equal to M 1+3 equipped with the usual Minkowski metric, which has 

components g^ v = diag(-l, 1, 1, 1) in an inertial coordinate system (x°, (x 1 ,^ 2 ,^ 3 )) d = f (t,x). Nonetheless, much of 
our discussion regarding the structure of the MBI system remains valid in an arbitrary spacetime. As is explained 
in detail in Section [4j the MBI equations can be expressed as 

(1.0.1a) dT = 0, 

(1.0.1b) dM = 0, 

where d denotes the exterior derivative operator, the Faraday tensor J 7 , which is a two-form, is the fundamental 
unknown, the Maxwell tensor Ai, which is also a two-form, is defined by 

(1.0.2) M = r ( 1 MBI) (^ + ^), 

* denotes the Hodge dual, d = f \{g l )^ K {g l ) nX J r c,r)^ r K\, ^(2) ils 1 )^^ ^^F^T K \ are the electromagnetic 
invariants, (g^ 1 )^ are the components of the inverse of the spacetime metric g, and £(mbi) d = (l + ^(1) ~ tfi)) 1 ^ 2 ■ 



^By spacetime, we mean a 4-dimensional time-oriented manifold M together with a Lorentzian metric g of signature (-,+,+,+). 
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Born and Infeld's contribution to the above system was their provision of the constitutive relation (1.0.2), while 
equations ( |1 O.lah - ( |1.0.1b[ ) were postulated^ in the 1860's by Maxwell |Max98aj . |Max98bj . We recall that in 
contrast to (1.0.2), Maxwell adopted the linear constitutive law Ai = *T '. Hence, we refer to the nonlinear system 



(1.0.1a), (1.0.1b), (1.0.2) as the "Maxwell-Born-lnfeld" equations, and the linear system (1.0.1a), (1.0.1b), M = *J r 



as the "Maxwell-Maxwell" equations. We summarize our main results here; they are thoroughly stated and proved 
in Sections [13] and [HI 

Main Results. The trivial solution to the MBI system on the 1 + 3 dimensional Minkowski space back- 
ground is globally stable. More specifically, if the initial data for the MBI system are sufficiently small 



as measured by the weighted Sobolev norm defined in (9.1.1), then these data launch a unique classical 



solution to the MBI system existing in all of Minkowski space. Moreover, these small-data solutions decay 
at the same rate as solutions to the linear Maxwell-Maxwell equations. In addition, the MBI system is 
hyperbolic^ in the interior of the field-strength regime in which its Lagrangian is real- valued. In particular, 
the system is locally well-posed in the aforementioned weighted Sobolev space. 

Remark 1.0.1. Certain large fields can cause I(mbi) to become complex- valued. For such fields, MBI theory is 
not even well-defined. However, as is discussed in Remark |6.8.2 the MBI equations are well-defined and hyperbolic 
for all finite values of the state-space variables (B,D), which are introduced in Section 6.8 In particular, the 



MBI system is well-posed for sufficiently regular initial data belonging to the interior of the region of state-space 
in which the equations are well-defined. 

Remark 1.0.2. Although we only discuss global existence to the future, our results apply to just as well to the 
past. Our notion of "future" is determined by assumption that the vectorfield c\ is future-directed. 

Recently, several scientific communities have expressed renewed interest in the MBI system for a variety of 
reasons. As an interesting example, we cite the works [Kie04aJ, [Kic04bJ, in which Kiessling has proposed a model of 
classical electrodynamics with point charges that has the promise of self-consistencjQ it is expected that the theory 
is well-defined without truncation, regularization, or renormalization. This theory couples a first-order guiding 
law for the point charges, whose corresponding relativistic guiding field satisfies a Hamilton-Jacobi type PDE, to 
the MBI field equations. In contrast to the case of the linear Maxwell-Maxwell equations, the electromagnetic 
potential^] of the solutions to the MBI system with non- accelerating point charge sources in Minkowski space can 
be chosen to be Lipschitz continuous; it is expected that this continuity property should remain true even for 
accelerating point charges, which would then allow for a well-defined coupling to the Hamilton-Jacobi theory. As 
a second example, we note that the MBI system has mathematical connections to string theory, for its Lagrangian 



(see (4.2.1)) appears in connection with the motion gauge fields (arising in the study of attached, open strings) on 



a D-brane (see e.g. |Gib03j ). 

As is true for Kiessling, our interest in the MBI system is further motivated by results contained in [Boi69] and 
|Ple70j , which show that it is the unique^] theory of classical electromagnetism that is derivable from an action 
principle and that satisfies the following 5 postulates (see also the discussion in [BB83J, |Kic04a] ) : 

(1) The field equations transform covariantly under the Poincare group. 

(2) The field equations are covariant under a Weyl (gauge) group. 

(3) The electromagnetic energy associated to a stationary point charge is finite. 

(4) The field equations reduce to the linear Maxwell-Maxwell equations in the weak field limit. 

(5) The solutions to the field equations are not birefringent (we will soon elaborate upon this notion). 

We remark that the Maxwell-Maxwell system satisfies all of the above postulates except for (3), and that the MBI 
system was shown to satisfy (3) by Born in [Bor33j. 



2 MaxwelFs formulation of electromagnetism was not presented using the framework of the Faraday tensor, nor that of the familiar 
electric field E and magnetic induction B\ rather, he used the structure of quaternions to write down a system of 20 equations in 20 
unknowns. The familiar "vector" formulation in terms of E and B was developed by Heaviside [Nah02 . 

^By "hyperbolic," we mean that there is a local energy estimate available that can be used to prove that initial data have a non-trivial 
development, and that furthermore, the system has finite speed of propagation. 

4 The motion of point charges in linear Maxwell-Maxwell theory is mathematically ill-defined under the usual Lorentz force law. 
^Recall that an electromagnetic potential is a one form A such that T = dA. 

6 More precisely, there is a one-parameter family of such theories indexed by |3 > 0, where (3 is Born's "aether" constant. 
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We would now like to further discuss postulate (5). Physically, it is equivalent to the statement that the "speed 
of light propagation" is independent of the polarization of the wave fields. Mathematically, it can be recast as a 
statement about the characteristic subset of the field equations. To flesh out this notion, we need to discuss some 



technical details. First, we remark that equation (|4.3.2b|), which reads H^^V u,Fk\ = 0, is equivalent to (jl.O.lb ) 



- (1.0.2) modulo equation (1.0.1a), where the tensorfield H^ vkX is defined in ( 4.3. la[ ). Now for each covector 



£ e T*M, the cotangent space of M at p, we consider the quadratic form x^iO d = H^^^^x- Because of the 
properties (4.1.10a) - (4.1.10c), which are also possessed by H^^, it follows that £ is an element of -?V(x(£))> 



the null space of x(0- The characteristic subset of T*M, which we denote by C*, is defined to be the set of all 
non-zero £ such that N(x(0) is strictly larger than span(£); i.e., 

(1.0.3) c; = f o 6 t;m \ n( x (0)/spM0 * 0}- 

As discussed in detail in |Chr00l Chapter 6], the set C* governs the local speeds of propagation of solutions to 
the MBI systenfl It is easy to see that C* is a conical set in the sense that if £ e C*, then any non-zero multiple 
of £ is also an element of C*. In general, this conical subset may have several different "sheets." However, in the 
case of the MBI system, there is a degeneracy resulting in the presence of only a single sheet (i.e., there is only 
one "null cone" associated to the MBI system); this is the mathematical characterization of "no birefringence." 
As we alluded to above, the Maxwell-Maxwell system also possesses this property. Moreover, in the case of the 
Maxwell-Maxwell system, C* exactly coincides with the gravitational null cone {£ e T*M | (5 _1 ) kA £ k £a - 0}. 
However, in a general nonlineaij^] theory, and specifically in the case of the MBI system, C* does not coincide 
with the gravitational null cone. Although we do not directly prove this fact in this article, we plan to discuss 
this issue in detail in a future publication, in which we will give a complete discussion of the geometry of the 



MBI system; see also the discussion in the proof of Proposition 8.4.2 and in [BB83j. We do however, as an aside, 
investigate a related issue that would be relevant if one wanted to couple the MBI system to the equations of 
general relativity. Namely, we prove that MBI system's energy-momentum tensor satisfies the dominant energy 



condition; see Lemma 8.1.1 Physically, this means that the speeds of propagation associated to the MBI system 
are no larger than the speeds associated to the gravitational null cone; i.e., the "speed of MBI light is less than or 
equal to the speed of gravity." Mathematically, this means that C* lies outside^] of the gravitational null cone. 

1.1. Comparison with related work. The core of our proof is based on a blend of ideas presented in (CK90 
and [ChrOOj; in [CK90], Christodoulou and Klainerman use methods similar to the ones used in this paper to 
analyze the decay properties of solutions to the linear Maxwell-Maxwell equations in Minkowski space, while in 
|Chr00j . Christodoulou provides a framework for deriving positive "almost conserved" quantities for nonlinear, 
hyperbolic PDEs that are derivable from a Lagrangian, of which the MBI system is an example. In short, using 
the methods of [ChrOOj . we are able to construct certain almost- conserved (in the small-solution regime) energies 
that have coercive properties nearly identical to those of the conserved quantities constructed in [CK90J. 

The aforementioned works and the present one are applications of a collection of geometric-analytic techniques 
that are applicable to a large class of hyperbolic PDEs derivable from a Lagrangian. These techniques are 
often collectively referred to as the vectorfield method. The term "vectorfield" refers to the fact that in typical 
applications, coercive quantities are constructed with the help of special vectorfields connected to the symmetries 
(or approximate symmetries) of the system. Originally introduced by Klainerman |Kla85] , [Kla86j in his analysis 
of small-data global solutions to nonlinear wave equations, the vectorfield method has grown into its own industry. 
As examples, we provide a non-exhaustive list of topics for which the vectorfield method has proven fruitful: 

• Global nonlinear stability results for the Einstein equations [BZ09j . |CK93| . |DH06j . |KN03j . [LR05] . 
pnoj . [RS09] . |SpelO| . 

• Small-data global existence for nonlinear elastic waves [Sid96| . 

• The formation of shocks in solutions to the relativistic Euler equations |Chr07j . 



More precisely, it is C p , the characteristic subset of T P M, the tangent space of M at p, that corresponds to the local speeds of 
propagation. C p is dual to C p in a sense defined in [ChrOOj . 
^More specifically, we mean quasilinear. 

^The dual picture perhaps more intuitively corresponds to the notion of MBI light traveling "more slowly than gravity:" C p lies inside 
of the gravitational null cone in T P M, which is defied to be {X 6 T p M \ g KX X K X x = 0}. 
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• Decay results for linear equations on curved backgrounds [AB09], [BIuOSJ . |DR.05| . |DR.08|, |DR.09| . |HollOj . 

• The formation of trapped surfaces in vacuum solutions to the Einstein equations |Chr09] , |KR09| . 

• Local existence and non-relativistic limits for the relativistic Euler equations without the use of symmetriz- 
ing variables | Spe 09a| , |Spe09b| , [SS10| . 

1.2. Difficulties in working with a four-potential. In various contexts during the study of linear Maxwell- 
Maxwell theory, authors commonly analyze the components of a four-potential A and its derivatives, rather than 
the Faraday tensor itself ( |Jac99l chapter 6] is a classic reference, and [Loi08j . jLoi09] are examples in the context 
of the Einstein-Maxwell system). Recall that a four-potential is a one- form A such that T = dA; the existence of 



such a one-form is guaranteed by (1.0.1a) and Poincare's Lemma. A is not unique, for any "gauge" transformation 
of the form A -> A + V7, where 7 is a scalar-valued function, preserves the relation T = dA. A well-known method 
of capitalizing on this gauge freedom is to work in the Lorenz gauge, which is the added condition 



(1.2.1) X7 K A K = 0. 

The advantages of the Lorenz gauge are discussed below. Of course, the viability of the gauge condition (1.2.1), 
which can be arranged to hold initially, depends on the fact that it is preserved by the flow of an appropriate 
version of the Maxwell-Maxwell equations (e.g. the system ( 1.2.5[ ) below) 



We would now like to discuss some subtle issues concerning the difficulties that may arise in an attempt to 
work with the Lorenz gauge when the electromagnetic equations are quasilinear. As in the remainder of the 
article, we assume in this section that (M,g) is Minkowski spacetime, and furthermore, that we are working in 
an inertial coordinate system. However, these assumptions have no substantial bearing on the issues at hand, 
for the same issues arise in any other spacetime (M,g) equipped with any coordinate system. We begin with a 
brief summary of the framework used for discussing an arbitrary nonlinear covariant theory of electromagnetism 
that is derivable from a Lagrangiarf^} If we choose to describe such a theory through the use of four-potentials 
A, then the Lagrangian *Jzf = *Jz?[vA] can be written as a function of VA. A very detailed elaboration of this 
discussion can be found in [ChrOOj; here, we only introduce the facts that are relevant to the issues at hand. The 
Euler-Lagrange equations (equations of motion) for such a theory can be written as 

(1.2.2) ^V C V^ = 0, {v = 0,1, 2, 3), 

where 

pp.* Cfi 

1 ' 1 0(Vc^O0(V„A") ' 

Note that h has a symmetry property that will be important for the construction of energies; it is invariant under 
the following simultaneous exchange of indices: 

(1.2.4) /X <->• z/, C^V- 

Before discussing the difficulties that arise in the quasilinear case, let us recall the advantages of using the Lorenz 
gauge in the linear Maxwell-Maxwell theory. In this case, there is a well-known, remarkable simplification that 



occurs in Lorenz gauge: the equations (1.2.2) can be written as a system of completely decoupled wave equations 



for the components of A. That is, in Lorenz gauge, the components of A are solutions to the following system: 

(1.2.5) S^GT^VcV^ = 0, (u = 0,1,2,3). 

Consequently, we have that 

(1-2.6) h^ = giil/ {g- l f\ 



We are slightly departing from the usual convention by referring to the Hodge dual of the Lagrangian, which we denote by *Jzf , as 
the Lagrangian; Jzf is a four-form, while *J£ is scalar-valued. 
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Because of the full decoupling at the quasilinear level, we can multiply both sides of (1.2.5) by the "seemingly 
non-geometric" quantity VoA/ (with the index v downstairs!), integrate over IR 3 , and integrate by part^jto show 
that the following energy £ is conserved for solutions to (1.2.5): 



(1.2.7) 



C,v=o 



In the language of [ChrOO] , the special structure of hfy in ( 1.2.6 ) is called separability] the existence of the conserved 
coercive quantity (1.2.7) is because of this additional structure, which is not typically present in the equations of 



a quasilinear theory. 

Let us contrast this to the case of the MBI system (or any other quasilinear perturbation of linear Maxwell- 
Maxwell theory derivable from a covariant Lagrangian). In the case of the MBI system in Lorenz gauge, it can be 
shown using (1.2.1) that the MBI equations can be written in such a way that 



(1.2. 



+ h (v 



where hjH, which has the symmetry property ( 1.2.4[ ), is a term that is of quadratic order in \/A in the small-solution 
regime. The corresponding system of PDEs is therefore 



(1.2.9) 



<7„*(<r 1 ) c ' 7 vcv^ + /$V ( V„^ = 



(i/ = 0, 1,2,3). 



Unfortunately, in general, it is not possible to simply multiply both sides of (1.2.9) by VqA Ui integrate over 



and integrate by parts; the tensorfield in (1.2.8) is not separable in general, nor in the particular case of the MBI 



system Note that this difficulty does not arise in the study of a single quasilinear wave equation; 

e.g., small quasilinear perturbations of the linear wave equation in Minkowski space preserve hyperbolicity and 
the availability of a basic L 2 energy estimate. 

One may attempt to resolve this difficulty by using the framework of energy currents developed in (Chr00| . A 
natural quantity that arises from an application of this framework is I(t), which is defined by 

(1-2-10) I(t) d =° f J° (MBI+Lorenz) (t,x) d 3 x, 

where the energy current J^Bi+Larenz) is defined by 



(1.2.11) 



u (MBI+Lorenz) 



L cf -^( V ^)(V A A ) + ^ ^S(V C ^)(V^ A ), 



(/x = 0, 1,2,3). 



As explained in detail in [ChrOOj and in Section [8j the current ( |1.2.11[ ) can be constructed by contracting a certain 
tensor, namely the canonical stress, against the vectorfield dt- The details of this construction do not concern us 
here. We remark only that the quantity I(t) is what one first tries to construct in an effort control solutions to 
the MBI system, i.e., during a proof of local well-posedness. On the one hand, it can be shown that ^jl(t) has 
one of the properties that is essential in order for it to be of use in analyzing the solution VA, namely that its 
time derivative can be bounded in terms of the L 2 norm of V A itself. More specifically, it can be shown that 



(1.2.12) 



d_ 

dt 



J(i) <C(||V/ l |[^)|| V A||| 2 . 



We remark that a quick way to see ( |1.2.12 ) is to use the equations (1.2.2) and the symmetry property (1.2.4) to 
show that \V nJ( MBI+Lorenz -) \ < C(|| V/i||l~)|V^4| 2 ; (1.2.12) then follows from the divergence theorem. Alternatively, 
one may multiply both sides of ( 1.2.9| ) by the "geometric" quantity VqA u (with the v index upstairs!) and integrate 
by parts with the help of the symmetry property (1.2.4), arriving at (1.2.12). 



These steps can alternatively be carried out using an energy current framework, similar to the energy current estimate (1.2.121 
described below. 
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However, we quickly run into a difficulty: ( 1.2.8 ) and (1.2.11) imply that in the small-solution regime, J (mbi+ Lorenz) 
is indefinite in V A : 



(1-2-13) J° { MBI + Loren Z ) = « £ 9*x( Vf^) ( V^ A ) + 0(|VA| 4 ). 

Therefore, is not a coercive quantity, and in particular, it is of no use in controlling the L 2 norms of solutions 



to (1.2.9). 



These difficulties are not fatal in the sense that the fundamental unknown is the Faraday tensor T = dA, and as 
explained in Section [9j we can construct suitable positive energies by working directly with T. More specifically, 
our energies control the combinations V ^A u - V U A^ of any four-potential satisfying T = dA, but they do not 
control the individual components V p,A v . Furthermore, there is an important advantage to working directly with 
T ■ our conditions for global existence depend only on physical quantities, and not on auxiliary mathematical 
quantities such as the values of a four-potential A. We remark that we do not know whether or not an alternate 
argumenip 2 ] could produce a positive quantity that is suitable for controlling the I? norm of the components V^A^ 
of solutions A to the MBI system (or, for that matter, any other covariant, quasilinear, non-separable system of 
electromagnetic equations derivable from a Lagrangian) in Lorenz gauge. That is to say, it is not clear whether 
or not the hyperbolicity of the equations is visible at the level of the components of VA in Lorenz gauge. If the 
answer turns out to be "no," then this would mean that the Lorenz gauge is not viable. We believe that this is an 
interesting question worthy of further investigation. 



1.3. Comments on the analysis. In this section, we summarize the main ideas of our proof. We first remark 
that all of the discussion in this section assumes that we have fixed an inertial coordinate system (t,x) on M, which 
is a global coordinate system in which the spacetime metric has the components = diag(— 1, 1, 1, 1). Throughout 
this article, we work directly with the Faraday tensor T , and thus avoid the aforementioned difficulties associated 
with choosing a gauge for the four-potential. To analyze J-, we use the framework of |CK90j and decompose it into 
its Minkowski null components. Before discussing the notion of null components, we first introduce the following 
foliations of Minkowski space: the family of ingoing Minkowski null cones C~ d = f {(r, y) \ \y\ + r = s}; the family 
of outgoing Minkowski null cones d = f {(r, y) \ \y\ - r = q}; and the spacelike hypersurfaces St d = f {(r, y) | r = t}, 

which intersect the null cones in spheres S r j d = f {(t, y) | r = t, \y\ = r}. All of these families of surfaces, will play an 
important role in this article. 

Now in a neighborhood of each spacetime point, there exists a null frame {L, L,e\,e2}, where L d = f dt - d r is an 

ingoing geodesic vectorfield tangent to the corresponding cone C~ , L d = f dt + d r is an outgoing geodesic vectorfield 
tangent to the corresponding cone normalized by the condition g(L,L) = -2, and the orthonormal vectorfields 
ei, e2 are tangent to the corresponding sphere S r t and normal to both L, L. At each point p where it is defined, the 
null frame forms a basis for T p M. The null components of the two-form T are then defined to be the following pair 
of one-forms a = afj 7 ], a = a[T] tangent to the S r ,t, and the following two scalar quantities p = p[F], a = c[JF] : 

a a = Fal, 
a = ^12, 

where we have abbreviated Tal - ^Ll^kXi F12 = ^\&\^Fk\-, etc; see Section M for more details. 



12 The authors in |CH03| claim to have overcome these difficulties, but their chain of reasoning in going from equation [CH03I Eqn. 
(3.8)] to equation CH03, Eqn. (3.11)] is difficult to follow; in particular, in equation CH03, Eqn. (3.10)], it is not clear whether 
the v index for the four-potential A is supposed to be "upstairs" or "downstairs," a distinction which is essential for establishing the 
positivity of their energies. 
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1.3.1. Linear analysis. The following decay properties, which can be expressed with the help of the null coordinates 
q d = r-t, s d = f r + t, where r d = f \x\, were shown in |CK90| for solutions to the linear Maxwell-Maxwell system arising 
from data with suitable decay properties at infinit}p*| 

• The worst decaying component is a, which decays like (1 + s) _1 (l + \q\) ~ 3 / 2 . 

• The fastest decaying component is a, which decays like (1 + s)~ 5 / 2 . 

• p and a each decay at the intermediate rate (1 + s)~ 2 (l + \q\) ~^ 2 . 

• Any derivative tangential to the outgoing cones C* (i.e. Vl> V eA ) creates additional decay of order (l+s) -1 , 
while the transversal derivative Vl creates additional decay of order (l+|g|) _1 , which is weaker than (l+s) -1 . 



In Section 14 (see also Proposition 11.0.9), we show that small-data solutions to the MBI system have exactly 
the same decay properties. Since the analysis of the linear theory also plays a key role in our analysis of the 
MBI system, we first discuss the basic strategy for establishing the aforementioned decay of solutions to the linear 
Maxwell-Maxwell system 



(1.3.1) AF = 0, d*F=0. 

We recall that for any two-form T, the corresponding Maxwell-Maxwell energy momentum tensor is 
/i q 9^ /^(Maxwell) 4s' i: kt _„ tt ttkA 



and that if T is a solution of ( [L3~I| ), then V ^Q^ MaxweU ) = Q, (v = 0,1,2,3). Furthermore, using the timelike 



conformal Killing fielq^ K, which has components K° = 1 + t 2 + (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , = 2tx 3 , (j = 1,2,3), we can 

, , i , jii ^iilMaxwelh^rpK- t-, ^.(Maxwell) ■ , . / _i \ K \ ^.(Maxwell) 

construct the energy current J? Maxwel t) = ~^ * ' -° ecause Q]iv J s symmetric, (g ) Q kX ~ 

0, and K is a conformal Killing field, it thus follows that 

(1.3.3) VuJft, , n = 0. 

v ' A 1 {Maxwell) 

Additionally, Qjj^ axwel1 ' has the following positivity property: for every pair of future-directed causal vectors X, Y, 
we have that Q^ axwel1 ^ x K Y x > 0. In particular, choosing d = f 5q, d = f ', it can be shown that (see Lemma 



8.4.1) 



(1-3.4) J° ( Max W eii) = ^{(1 + Q 2 )h\ 2 + (1 + « 2 )l«| 2 + (2 + Q 2 + s 2 )(p 2 + a 2 )}, 

where a, a, p, and a are the null components of T . If we define the energy £ > by 

(1-3-5) £ 2 (t) = J J^ Maxwell ^(t,x)d 3 x, 



then it follows from (1.3.3) and the divergence theorem that £{t) is constant in time if it is initially finite: 



(1.3.6) |(£ 2 (t)) = 0. 



The various weights in (1.3.4) are the first hint that different null components of T have different L°° decay 
properties. However, the full proof of decay requires that we commute the Maxwell equations with various 
conformal Killing fields and apply the global Sobolev inequality. Let us explain what we mean by this. Given any 



solution T of (1.3.1) and any conformal Killing field Z, it can be shown that LzT is also a solution to the linear 
Maxwell-Maxwell equations. Here, LzJ~ is the Lie derivative of T with respect to the vectorfield Z. Iterating this 
process, we conclude that C^J- is a solution, where / is a multi-index, and CJ Z i s shorthand notation for iterated 
Lie derivatives with respect to vectorfields Z e Z. In this article, the relevant set of conformal Killing fields Z 
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The finiteness of ||(B,Z))|| H 3 is sufficient, where (B,D) is the electromagnetic decomposition of the data for T described in Section 



6.8 



and the weighed Sobolev norm H\ is defined in Definition 



9.1.1 



"^Recall that a conformal Killing field is a vectorfield X that satisfies V M X„ + Vi/-X/x = 4>x9fiv for some scalar- valued function <j> x 
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consists of: the 4 translations = d^, (/i = 0, 1,2,3); the 3 rotations &(jk) = %jdk ~ x kdj, (1 < j < k < 3); the 3 

def , ^ r, / . , „ „\ , -,■ r. , , ry dei 



Lorentz boosts &roj) - ~^dj -Xjdt, (j = 1,2,3); and the scaling vectorfield S = x K d K . Furthermore, as in (1.3.6), 
the weighted L 2 norms of the various null components of the C^J- are constant in time. Now in order to derive 
L°° decay, we need to connect these weighted L 2 norms of C^J- to weighted L°° norms of T . This is exactly what 
the global Sobolev inequality provides; see Proposition 11.0.9 for the details. 

Let us also discuss the heuristic mechanism for the fact that derivatives tangential to the C* (i.e. Vl, V eA ) 
have better decay properties than the transversal derivative Vz,. As examples, we consider the outgoing vectorfield 

L d = dt + d r and the ingoing vectorfield L d = f dt - d r , where d r denotes the radial derivative. Simple algebraic 
computations lead to the identities 



(1.3.7) 



L = 



S-u a Q 



(Ott) 



L = — 



S + u a Q 



(Ott) 



where w 1 d = f x i /r, and q, s are the null coordinates mentioned above. Therefore, if we have achieved good control of 



the quantities Vs-?"" and Va n .,^, then the formulas (1.3.7) suggest that we can achieve even better control of the 



outgoing derivative V because of the favorable denominator s' 1 . On the other hand, the transversal derivative 
Vl-F features a less favorable denominator q~ l . More specifically, the q^ 1 term fails to provide decay in the "wave 
zone" r « t, while in the entire region {t > 0} we have that s _1 < min{i _1 , r -1 }; i.e., decay in s implies decay in r 
and t. 

1.3.2. Nonlinear analysis. We now outline the key differences between the proof of decay of solutions to the linear 
Maxwell-Maxwell equations, and the proof of the global existence of and decay of solutions to the MBI system in the 



small-data regime. To analyze solutions to the MBI system, the "working form" of which is given below in (4.3.2a) 



(4.3.2b), we will use the same Minkowski null decomposition of the Faraday tensor described above. In particular, 



in order to derive our desired estimates, we do not need to use the characteristic geometry of the MBI 
system; in using the "wrong" Minkoswkian geometry (which has the advantage of relative simplicity), we are 
deviating from the correct MBI geometry (which is governed by the reciprocal Maxwell-Born-lnfeld metric {lr 1 )^ u 
defined in ( 8.4. 14[ )) by small error terms that are controllable. Now like the Maxwell-Maxwell system, the MBI 
system has a corresponding divergence-free energy- momentum tensor , which is given below in (8.1.4); the 



availability of this tensor is a well-known consequence of the fact that the MBI Lagrangian *J£(mbi) ( see (4.2.1)) 
depends covariantly on only the metric g and the field variables T . This tensor can be used in conjunction with the 
vectorfield K to estimate the weighted L 2 norm of the solutions J- to the MBI system. However, to estimate the 
weighted L 2 norm of C^T, we need a different tensor, which is described by Christodoulou in detail in [Chr00| : the 



canonical stress Q u , which is defined below in (8.3.4). Now in the case of the linear Maxwell-Maxwell system, the 



fi(Maxwell) 



constructed 



canonical stress corresponding to the C^J- coincides with the energy momentum tensor Q„ 
out of the jC^F, but in a general nonlinear theory, the two tensors differ. The important point is that the CJ Z J~ are 



solutions to the linearized equations (8.2.1a) - (8.2.1b), which are derivable from a linearized Lagrangian ££ (see 



(8.3.1)) depending on the metric g, the linearized variables T d = C^T, and also the background T. Although the 
dependence of Jz? on the background precludes the availability of a divergence-free tensor for the linearized system, 
we may nevertheless use Christodoulou's framework to construct the tensor Q u . Although Q^ u is in general not 
even symmetric, nor is Qv divergence-free, the role that Q% plays in the analysis of the linearized equations is 
roughly analogous to the role played by the energy momentum tensor in the original equations: V^Qt, though 
non-zero, is of lower order (in terms of the number of derivatives), and furthermore, Qv possesses some positivity 
properties under contractions against certain pairs (£,,X) of timelike (covectors, vectors). 

Once we have Q% y we can again use the vectorfield K construct energies £jsr[J-(t)], which are a sum over |7| < N 
of the energy of C^T, that are analogous to the energies (1.3.5) defined in the Maxwell-Maxwell case; the precise 
definition is given in (9.0.7) below. However, in the MBI system, Spfi^i^)] is n °t constant. Additionally, in the 



nonlinear problem, the q, s weighted factors appearing in the expression £^[3^(1;)] are not manifestly uniform, but 
instead depend on the solution T itself. For these reasons, it is convenient to introduce a norm \ F(t) \c z ;N 



whose q, s weights don't depend on T\ see (9.0.6). In order to compare the two quantities, we establish inequality 
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(9.0.11b), which shows that in the small-solution regime, £N[!F(t)] w^f- T{t) ^c z ;N ■ The crux of the global 



existence proof is the following: even though E^Tit)^ is not constant, we are nevertheless able to derive an a- 
priori bound for \ T{t) \c z \N which shows that it remains uniformly small on any time interval of existence for the 
solution. According to the continuation principle of Proposition 13.0.4, such an a-priori bound for -jf J-{t) ^c z ;N 
implies global existence. 

Now in order to estimate \ J~(t) ^c z ;N, we need to handle the numerous "error" terms arising in the ex- 
pression for f t (£ 2 N [F(t)]). The first source of error terms was alluded to above, namely that the divergence of 
Q^v is non-zero. The second source comes from the fact that the C^J- are solutions to the linearized equations 
with inhomogeneous terms: many inhomogeneous "error" terms arise from commuting the operator CJ Z through 
the equation (4.3.2b); see Proposition 8.2.1 This commuting is accomplished through the use of modified Lie 



derivatives Cz, which are equal to ordinary Lie derivatives plus a scalar multiple of the identity; see Definition 



5.0.2 and Lemma 7.0.4. A careful analysis of the special structure of the error terms (which are discussed in the 



next section), in conjunction with the global Sobolev inequality, leads to the a-priori estimate (12.0.4), which is 
valid in the small-solution regime, and which we restate here for convenience: 



(1.3. 



t Ht) \l z ;N< C{ t ^(0) e Cz ;N + Y~^2 + ^T) %l 2 . N dr). 



Applying Gronwall's inequality to (1.3.8), we thus conclude the desired result: -jf T{t) \c z \N is globally bounded 
in time, if \ J~(0) \c z \N is sufficiently small. In addition, we remark that the aforementioned decay properties 
of the solution J- are a by-product of the previous analysis. More specifically, the decay properties of J- follow 
directly from the global bound on jf- F(t) \c z \N and the global Sobolev inequality (Proposition 11.0.9). 

We have a few final comments to make concerning the smallness of the data. The initial data consist of a 
pair of one-forms (D,B) that are tangent to the Cauchy-hypersurface So, and that satisfy the constraints (which 
are familiar from linear Maxwell-Maxwell theory) V a D a = V a -B a = 0. Here, V denotes the Levi-Civita connection 
corresponding to the first fundamental form g of So (see Section [31). As is described in Section 6.8 .F(O) can be 

constructed out of (D,B). However, the quantity -jf !F(Q) %c z ;N, the smallness of which is required to close the 
global existence argument, involves derivatives of T that are normal to So (i.e., time derivatives). However, by 
repeatedly using an appropriate version of the MBI system, the normal derivatives of T along So can be written 
in terms o f th e tangential derivatives (i.e. spatial derivatives) of (D,B). Consequently, as is explained in detail 
in Section 9.1, it is possible to devise a smallness condition involving only the data (D,B) and their tangential 
derivatives, from which the smallness of -jf J-(0) \z z \N necessarily follows. This allows for a "proper" formulation 
of the small-data global existence condition of Theorem [T] in terms of quantities inherent to the data. 

1.3.3. The error terms. Let us now make a few remarks concerning the many error terms that arise in our study 
of ^(<?^[.F(t)]), since the study of these error terms is at the heart of our analysis. In the small-solution regime, 
the MBI system is a cubic quasilinear perturbation of the linear Maxwell-Maxwell system. It is well-known that 
for linear hyperbolic PDEs whose solutions possess the decay properties of solutions to the Maxwell-Maxwell 
system, cubic perturbations^ do not destroy the existence of small-data global solutions. In fact, a much shorter 
proof of small-data global existence could be provided by using the vectorfield dt in place of the vectorfield K in 
our construction of the energies. However, in order to show that small-data MBI solutions have the same decay 
properties as solutions to the linear Maxwell-Maxwell system, we make full use of the vectorfield K, together with 
an algebraic property of the MBI system: its nonlinearities satisfy the null condition. The null condition, which 
was first identified by Klainerman |Kla86] in the context of nonlinear wave equations, is a collection of algebraic 
properties that are satisfied by special nonlinearities. Roughly speaking, when a nonlinearity satisfies the null 
condition, the worst kind of terms (from the point of view of decay) are not present. More specifically, in the case 
of the MBI system, the expression for ^(£^[.F(i)]) involves quartic terms in T and its iterated Lie derivatives 
jC^-^") multiplied by weights in q and s arising from the vectorfield K and its covariant derivative S/K . Because 
these terms are fourth order, we do not need to perform a fully detailed null decomposition in order to prove our 
desired estimates. That is, there is room for imprecision; we only prove estimates that are sufficient recover the 
full decay properties possessed by solutions to the linear Maxwell-Maxwell system. Let us summarize the version 



'We are assuming that the perturbations involve only 1 or fewer derivatives, and that the perturbed system is also hyperbolic. 
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of the null condition that we show is satisfied by the MBI system (see Section 10 for the details): for those terms 
involving a weight of s or 1 + s 2 , at most two of the four factors correspond to the worst decaying components 
alCgJ-]. It is also true that for those terms involving a weight of q or 1 + q 2 , at most three of the four factors 
correspond to the worst decaying components alC^J-]. However, we do not make use of the availability of the one 
"good" factor since our estimates close without it. 

1.3.4. The large-data well-posedness of the MBI system. Finally, we would like to make a few remarks about the 



local existence proof that is briefly sketched in Section 13 This result is interesting in itself because it shows 



the following fact, which is arguably not manifest: the MBI system's initial value problem is well-posed in every 
field-strength regime in which its Lagrangian is real- valued, i.e., in every regime in which the theory is well-defined. 



The crucial estimate in this regard is contained in Proposition 8.4.2, which shows that it is always possible to 
construct an energy current for the linearized equations with positivity properties that are sufficient to deduce 
local existence in the weighted Sobolev space of relevance for our global existence result. This fact is strongly 
related to the internal geometry of and the hyperbolicity of the MBI system (i.e., the characteristic subsets), which 
will be explored in detail in an upcoming article by the author and his collaborators. We remark that the energy 
current we use for deducing the local existence result is constructed by contracting the canonical stress against a 
suitable "multiplier" vectorfield Xi oca i, so that it differs from the current used in our small-data global existence 
proof; the vectorfield K may not be a suitable multiplier for deducing large-data local existence. 

For an alternative proof of the large-data well-posedness of the MBI system's initial value problem, one may 
consult [BreOlj (see also jSer04j ). In this work, Brenier "augments" the MBI system by taking as his 10 unknowns 



the non-trivial components of the electromagnetic quantities (see Section 6.8) B, D, P, and h. Along the "MBI 
submanifold," P coincides with the Poynting vector (P = B*D), and h coincides with the 00 component of the MBI 
energy-momentum tensor Qn^BI) ( see (8-1.4)), but for general augmented MBI solutions, P and h are independent 



unknowns. To compensate, the additional evolution equations |Bre04l Eqn. (1.9)] and [Bre04, Eqn. (1.10)], which 
are redundant for solutions belonging to the MBI submanifold, were added to the MBI system (i.e., so that there 
are 10 equations for the 10 unknowns). From the point of view of hyperbolicity, the most important feature of this 

augmented system is that the function S(B,D,P,h) d = f 1+ ^ + ^ + ^ , which coincides with a constant multiple 
of the quantity h = Q^mbi) ^ or solutions constrained to the MBI submanifold, satisfies the properties of a smooth, 
strictly convex entropy function of the augmented variables (B,D,P,h). Thus, using the general framework of 
hyperbolic conservation laws (see e.g. [DaflO]), it follows that there exists a change of state-space variables in 
which the augmented MBI system becomes symmetric hyperbolic. For symmetric hyperbolic systems, there exists 
a well-established theory of well-posedness based on energy estimates (see e.g. [CH89J, [DaflOj . |Fri54j . |Lax06j , 
|Maj84| , [Ser99] ). 

1.4. Outline of the article. The remainder of the article is organized as follows: 

• In Section 2 we collect together much of the notation that is introduced throughout the article. 

• In Section 3, we recall some basic facts from differential geometry. 

• In Section 4, we provide a detailed introduction to the MBI system. 

• In Section[5j we discuss the collections of Minkowski conformal Killing fields that play a role in our analysis. 
We also introduce modified Lie derivatives, which have favorable commutation properties with the MBI 
equations. 

• In Section [6j we introduce the null frame and the null decomposition of a tensor. We then decompose 
the MBI system relative to a null frame. We also introduce several electromagnetic decompositions of the 
Faraday and Maxwell tensors. 

In Section [7| we provide some commutation lemmas that will be used throughout the remainder of the 



article, especially in Section 11 



In Section [8j we discuss the energy-momentum tensor associated to the MBI system and the canonical 
stress tensor associated to the equations of variation. 

In Section [9j we introduce the norms, seminorms, and energies that will be used in the proof of our main 
theorem. 



In Section 10, we perform a partial null decomposition of the nonlinear error terms that appear in the 



expression for the time derivative of the energy. It is here that the null condition is revealed. 
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In Section 
L°° bounds. 
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we discuss the global Sobolev inequality, which connects weighted L 2 bounds to weighted 



In Section 12, we prove the a-priori bound (1.3.8), which is the most important inequality in the article. 



In Section 13, we briefly discuss local existence for the MBI system. We also recall the availability of a 
continuation principle, which provides criteria for the existence of a global classical solution. 



In Section 14, we combine the results of Sections 12 and 13 in order to establish our main theorem. 



2. Notation 

In this section, we collect together for convenience much of the notation that is introduced throughout the article. 



2.1. Constants. We use the symbol C to denote a generic positive constant that is free to vary from line to line. 
In general, C can depend on many quantities, but in the small-solution regime that we consider in this article, C 
can be chosen uniformly. Sometimes it is illuminating to explicitly indicate one of the quantities O that C depends 
on; we do by writing Cq or C(0). If A and B are two quantities, then we often write 

A<B 

to mean that "there exists a C > such that A < CB." Furthermore, if A < B and B < A, then we often write 

AmB. 

2.2. Indices. 

• Lowercase Latin indices a,b,j,k, etc. take on the values 1,2,3. 

• Greek indices K,X,fJL,f, etc. take on the values 0,1,2,3. 

• Uppercase Latin indices A, B etc. take on the values 1,2 and are used to enumerate the two orthogonal 
null frame vectors tangent to the spheres S r t- 

• Indices are lowered and raised with the Minkowski metric g^ u and its inverse (g 1 ) fiu ■ 

• Repeated indices are summed over. 

• We sometimes use parentheses to distinguish indices that are labels from coordinate indices; e.g., the "0" 
in T( ) is a labeling index. 

2.3. Coordinates. 

• {x^ 1 } u=o,i,2,3 are the spacetime coordinates; in our fixed inertial coordinate system only, we use the 

. . . dof n / 1 9 '3 \ 

notation t = x , x = \x ,x , x ). 

• Relative to our inertial coordinate system, r d = f |x| d = \J(x 1 ) 2 + (x 2 ) 2 + (x s ) 2 denotes the radial coordinate. 

• q d = f r - t, s d = f r + t are the null coordinates. 

2.4. Surfaces. Relative to the inertial coordinate system: 

• Cg d = f {(r, y) | \y\ + t = s} are the ingoing null cones. 

• Cq d = f {(t, y) | \y\ - 1 = q} are the outgoing null cones. 

• T>t d = f {(t, y) \ t - t} are the constant time slices. 

• SV,i d = f {(t, y) | t = t, \y\ = r} are the Euclidean spheres. 

2.5. Metrics and volume forms. 

• g denotes the standard Minkowski metric on M 1+3 ; in our fixed inertial coordinate system, g^ u = diag(-l, 1, 1, 1). 

• e denotes the standard Euclidean metric on M 1+3 ; in our fixed coordinate system, e^ = diag(l, 1, 1, 1). 

• e _1 denotes the inverse of the standard Euclidean metric on M 1+3 ; in our fixed inertial coordinate system, 
(e- 1 )^ = diag(l, 1,1,1). 

• g denotes the first fundamental form of S^; in our fixed inertial coordinate system, g = diag(0, 1, 1, 1). 
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• g denotes the first fundamental form of S r ^\ relative to an arbitrary coordinate system, 
&nu~ 9fj.f + \{j- t \iL. v + LluLv), where L, L are defined in Section 2.9 

• C/j,vk\ - Idetg^^l/ivnX] denotes the volume form of g; [0123] = 1 = -[1023], etc. 



£ VK \ = (-llvk\{T(q))^ denotes the volume form of g, where T( ) is defined in Section 
ii-iu - ^£fiun\L. h 'L x denotes the volume form of g . 



21 



2.6. Hodge duality. For an arbitrary two-form T^ v ■ 

• *^nv = lg^jiguu^ UKX T K x = -\\det g\~ 112 g^g vV [fwKX\T K \ denotes the Hodge dual of with respect to 
the spacetime metric g. 

2.7. Derivatives. 

In an arbitrary coordinate system {^ M } M =o,i,2,3, = -J^, V M = V a ■ 

dxf 1 

V denotes the Levi-Civita connection corresponding to g. 

V denotes the Levi-Civita connection corresponding to g. 
f denotes the Levi-Civita connection corresponding to g . 

In our fixed inertial coordinate system, d r = u; a d a denotes the radial derivative, where uj j = x J 'jr. 
V x denotes the differential operator X K V K . 

If X is tangent to £t, then Vj^ denotes the differential operator X K V_ K . 
If X is tangent to S r ,t, then 'fx denotes the differential operator X K f K 



Vfn^U denotes the n th covariant derivative tensorfield of the tensorfield U. 



— (n)^ denotes the n th Sj-intrinsic covariant derivative tensorfield of a tensorfield U tangent to the hyper- 
surfaces S t . 

f( n ) U denotes the n th <SV,t -intrinsic covariant derivative tensorfield of a tensorfield U tangent to the 
spheres S r> f. 

div U = g X V\U K denotes the intrinsic divergence of a vectorfield U tangent to the hypersurfaces E| . 

( curl U) u = e UK X \7 K U X are the components of the intrinsic divergence of a vectorfield U tangent to the 
hypersurfaces T, t . 

dfv U =g^ \7\U K denotes the intrinsic divergence of a vectorfield U tangent to the spheres S r j- 
cv/rl U = j K x ^K.U x denotes the intrinsic curl of a vectorfield U tangent to the spheres SV,*- 
Cx denotes the Lie derivative with respect to the vectorfield X. 

[X, YY = (C X YY = X K d K Y^ - Y K d K X^ denotes the Lie bracket of the vectorfields X and Y. 

For Z e Z, Cz - £>z + 2c^ denotes the modified Lie derivative, where the constant cz is defined in Section 



• C^U, and £jU, V^U respectively denote an \I\ th order iterated Lie, iterated modified Lie, and iterated 
covariant derivative of the tensorfield U with respect to vectorfields belonging to the set A; f® U is an 
iterated intrinsic (to the spheres S r> t) covariant derivative of U with respect to rotation vectorfields. 

2.8. Minkowski conformal Killing fields. 

Relative to the inertial coordinate system {x M } A1= o,i,2,3 = (t,x) ■ 

T(n) = c^, (/x = 0, 1, 2, 3), denotes a translation vectorfield. 
Q(jk) = Xjdk - Xkdj, (l<j<k< 3), denotes a rotation vectorfield. 
$7( j) = -tdj - Xjdt, (j = 1,2,3), denotes a Lorentz boost vectorfield. 

5 = x K d K denotes the scaling vectorfield. 

K(n) = -2x M 5 + g K \x K x x d fM , (fi = 0, 1, 2, 3) denotes an acceleration vectorfield. 
K - i'T(o) + T( ) denotes the Morawetz vectorfield. 

6 = { r (/i)}o<At<3- 
C = {^(12)^(13)^(23)}- 

T, O, and Z are the Lie algebras generated by ^,0, and Z respectively. 
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• For Z e Z, ( 'tt^u = V ' pZ v + V V Z^ = czg^v, denotes the deformation tensor of Z, where cz is a constant. 

• Commutation properties with the Maxwell-Maxwell term [(g' 1 )^ '(g~ 1 ) vX - (g' 1 ) f " X (g' 1 ) L ' K ]v ^kX ■ 

&z{[(9- 1 )' ut (g- 1 r x - (g-'rHg-'r^^x} = [(g-'r^g- 1 )^ - Gr^GrYlv ^ j z t kX . 

2.9. Null frames. 

• L d = f dt - d r denotes the null vectorfield generating the C~ and transversal to the . 

• L d = f dt + d r denotes the null vectorfield generating the . 

• e A , A =1,2 denotes orthonormal vectorfields spanning the tangent space of the spheres SV,t- 

• The set C d = f {L} contains only L. 

• The set 7 = {L, ei,e2} denotes the frame vector fields tangent to the C+. 

• The set U d = {L,L,e\,e2} denotes the entire null frame. 

2.10. Null frame decomposition. 

• For an arbitrary vectorfield X and frame vector U tU, we define Xjj d = f X K U K , where X^ d = f g^ K X K . 

• For an arbitrary vectorfield X = X K d K = X L L + X-L + X A e A , where 

vL _ 1 v vL _ 1 v ~vA _ \ r \ r v 
A - ~2 A L> A— - -^A^, A - -A-Ai -A- A ~ -*-e A - 

• For an arbitrary pair of vectorfields X, Y : 
g(X,Y) = X K Y K = -\X l Yl- \X l Y l + 5 AB X A Y B . 

If T is any two-form, its null components are 

. p = \T Xr L k L x . 



1 IkX rr~ 



2* 

b 

2.11. Null Forms. For arbitrary two-forms T, Q ■ 

• Q W (F, G) = T kX Q kX = -5 AB a A [T]a B [G] - 5 AB a A [g]a B [T] - 2p[T]p[Q] + 2a[T]a[Q]. 

• Q(2){F,Q) = *? kX QkX = t AB a A {T\a B \G\ + t AB a A [g}a B [F] - 2a[T]p[Q] - 2p[f]a[g]. 

2.12. Electromagnetic decompositions. Given a two-form T and its associated MBI Maxwell tensor 
■M./XU = ^jf5/)(*^/n/ + ^(2)^>^)) its electromagnetic components relative to an arbitrary coordinate system are 

• = F^kT^ . 

• = - ^kT^j . 
. D^ = -*M^T« Q) . 

. £ M = ^ K S K . 

2.13. Norms and energies. For an arbitrary tensor U of type ( n ), and A e 0, Z} : 

• I^ILn = £i^I^I 2 - 

• l^ ;JV -E|/|^|V^| 2 . 

• l^ o;JV = f E|7|,iv|^ U\\ 

For an arbitrary type (°) tensor F, and V, W e {/2,T,ZY} : 

• |F|vw = Zvzv,wzw \V K W X F KX \. 

\V K W x W\7 y F KX \. 

For an arbitrary two-form T with null components a, a, p, &; and A £ : 

• \jr | 2 = (1 + g 2 )|a| 2 + (1 + s 2 )|d| 2 + (2 + g 2 + s 2 )(p 2 + d 2 ). 
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• + T \\ A . N = E|/|<tv + \ 2 . 

• f Ht) !L;iV =7 R 3 + Ht,X_) \l z , N d 3 X. 

For an arbitrary tensorfield U defined on the Euclidean space S with Euclidean coordinate system x ■ 

• W\\ 2 HN = Zn^oM 1 + N 2 ) (5+n) |V (n) ?7(x)| 2 d 3 x is a weighted Sobolev norm of [/. 



• ll^llcrjv d = f E^=o su PxeE ( 1 + l^| 2 )^ +n ' ) lY( n ^(x)| 2 is a weighted pointwise norm of U. 

For arbitrary two-forms T and J 7 : 

• H^^Vu^/tA) where H^ ukX depends on J 7 , is the principal term in the equations of variation (8.2.1b). 

• Ctv = H^ kX T kX T uC - ^H^^T^T^ is the canonical stress tensor. 



• J£[F] = -Q^yK is the energy current (used during the proof of global existence) constructed from the 
variation P, the background JF, and the Morawetz-type vectorfield K - |{(l + s 2 )L+(l + q 2 )L\- 

• £^[.F(i)] = E|/|<jv fu 3 jpi^zF] d 3 x is the square of the order N energy of T . 

2.14. Function spaces and the regularity of maps. 

• Hf is the set of all distributions / such that ||/||#jv < oo. 

• is the set of all functions / such that ||/|| c jv < oo. 

• C fc ([0,T) x M 3 ) denotes the set of /c-times continuously differentiable functions on [0,T) xR 3 . 

• If X is a function space, then C fc ([0,T),X) denotes the set of A;-times continuously differentiable maps 
from [0, T) to X. 



3. Geometry 

In this section, we recall some basic facts from differential geometry that will be used throughout the article. 



3.1. Inertial coordinate systems, the spacetime metric, and the Riemannian metric. In Minkowski 
space, there exists a family of global coordinate systems, which we refer to as inertial coordinate systems, in which 
the metric g^ v and its inverse (g -1 )^ have the following components: 



(3.1.1) ^ = (p- 1 )^ = di«g(-l,l ) l,l)- 

It will be convenient to carry out calculations and to define various tensors relative to an inertial coordinate 
system. Therefore, we fix a single inertial coordinate system {2^)^=0,1,2,3 on Minkowski space. For 
the remainder of the article, when we decompose tensors with respect to an inertial frame, it will 
always be relative to the frame corresponding to this fixed inertial coordinate system. When working 
in this coordinate system, we often use the abbreviations 

(3.1.2a) x od = cf i, x^(x\x 2 ,x 3 ), 

(3.1.2b) d p*JL, d t ^d -T {0) . 



We recall the following partitions of T p M and T*M induced by g. 
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Definition 3.1.1. Vectors X e T p M are classified as timelike, null, causal, spacelike as follows, where g{X, X) = 
9,xX K X x : 

(3.1.3a) g(X,X)<0 (timelike), 

(3.1.3b) g(X,X) = (null), 

(3.1.3c) g(X,X)<0 (causal), 

(3.1.3d) g(X,X)>0 (spacelike). 

Furthermore, relative to our fixed inertial coordinate system, X is classified as future-directed or past-directed 
as follows: 

(3.1.4a) X°>0 (future-directed), 

(3.1.4b) X° < (past-directed). 

Covectors £^ are defined to have the same classification as their metric dual X^ d = f (g 1 )^^- We sometimes 
refer to £ as the g-dual of X in order to emphasize that this notion of duality depends on g. 

In order to measure the size of various tensor, it is convenient to introduce a Riemannian metric on M 4 . A 
natural choice is the Euclidean metric e, which has the following components relative to an arbitrary coordinate 
system: 

(3-1-5) e M „ = f g^u + 2(T( )) At (T( 0) ) I/ . 

In the above formula, T( ) is the "time translation" vectorfield, which is defined to coincide with dt in our inertial 
coordinate system. Therefore, relative to this coordinate system, the metric e and its inverse e _1 have the following 
components: 

(3.1.6a) e^ v = diag(l, 1, 1, 1), 

(3.1.6b) (e-y^diagtl, 1,1,1). 

We now define the aforementioned tensorial norm. 
Definition 3.1.2. If U is a tensor of type f n ), then we define the norm | • | > of U by 

(3.1.7) M^Ke- 1 )^ 1 "^- 1 )^^ 

3.2. Lie derivatives and covariant derivatives. Given any pair of vectorfields X, Y, we recall that relative to 
an arbitrary coordinate system, their Lie bracket [X, Y] can be expressed as 

(3.2.1) [X,YY = X K d K Y^-Y R d K X li . 
Furthermore, we have that 

(3.2.2) C X Y=[X,Y], 

where C denotes the Lie derivative operator. Given a type tensorfield U, and vectorfields Y^\, ■■■Yr m \, the 



Leibniz rule for C implies that (3.2.2) generalizes as follows: 



(3.2.3) {CxU){Y {1)r -,Y {m) ) = X{U(Y w ,---,Y^ 

i=l 
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Remark 3.2.1. The Lie derivative operator does not commute with the raising and lowering of indices via the 
metric g. Thus, in order to avoid confusion, we use the convention that Lie derivatives are applied to two-forms 
J-^v with both indices down. In particular, the quantity CzF is understood to be a two-form with the indices 

down, and C z ^ v = f {g- l Y K {g- l ) vX CzF K \. 

There is a unique affine connection V, which is known as the Levi-Civita connection, that is torsion- free and 
compatible with the metric g. These properties are equivalent to the requirement that the following identities hold 
for all vectorfields X, Y, Z : 

(3.2.4) V X Y-V Y X = [X,Yl 

(3.2.5) Vx{g(Y, Z)} = g(V x Y, Z) + g(X, V Y Z). 

Furthermore, given a type (°) tensorfield U, and vectorfields 5?i),"">^?m)> the Leibniz rule implies that 

m 

(3.2.6) (vx[/)(y (1 ),-^y (m) ) = x{[/(y (1) ,--.,y^ 

i=i 



We remark that relative to an arbitrary coordinate system, (3.2.5) is equivalent to 



(3.2.7) VaS^ = 0, (A, /i,i/ = 0,1,2,3). 

Furthermore, in our inertial coordinate system on Minkowski space, if U is any type ( n ) tensorfield, then 
V aU IM ., Ul '" Vn = d a U,., ,. uv " Un . In the above formulas and throughout the article, we use the notation 

(3.2.8) Vx d =° f A^V K . 

The Riemann curvature tensor R(-,-)- is defined by the requirement that the following identities hold for all 
vectorfields X, Y, Z : 

(3.2.9) V X Vy2 - VyVjZ = R(X, Y)Z + V [x ,y]Z. 

In Minkowski space, R(X,Y)Z = 0. 

The following standard lemma gives a convenient formula relating Lie derivatives and covariant derivatives. 

Lemma 3.2.1. |Wal84| 

Let X be a vectorfield, and let U be a tensorfield of type ( Then CxU can be expressed in terms of covariant 
derivatives of U and X as follows: 

(3.2.10) {c x u)^ m y&xii)^ + v^^" + - + ^.., ra ,r-"v Mra r 



□ 



It follows that 

(3.2.11) Cxg^u = (X) vr^, 
where 

(3.2.12) Wvr^^V^ + V^ 
is the deformation tensor of X. 
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3.3. Volume forms and Hodge dual. There is a canonical volume form e^ UK \ associated to the metric g. 
Relative to any local coordinate system, we have that 

(3.3.1) e KXlM „ = \det(g)\ l ' 2 [K\fiu], 

(3.3.2) e KA ^ = -\det(g)\- 1/2 [K\^iy], 

where [kA/xz^] is totally antisymmetric with normalization [0123] = 1. It can be checked that the covariant derivative 
of the volume form vanishes: 

(3.3.3) V^a^O, (ft k, A, /i,i/ = 0,1,2,3). 

The Hodge dual operator, which we denote by *, plays a fundamental role throughout our discussion. 
Definition 3.3.1. If T is any two-form, then its Hodge dual *T is defined as follows: 

(3.3.4) *T^ U = f ^e KX ^F KX . 

3.4. HtiS r t, and the first and second fundamental forms. 

Definition 3.4.1. The following two classes of spacelike submanifolds of Minkowski space, which we define relative 
to the inertial coordinate system {x At } At= o,i,2,3 will play a role throughout the remainder of the article: 

(3.4.1) £i d = cf {(r,y)|r = t}, 

(3.4.2) S rit d ^{(r,y)\T = t,\y\ = r}, 

where \y\ d = f yj (y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 - We refer to the T, t as "time slices," and the S r j as "spheres." 

The future-directed normal to the St is the time translation vectorfield T(o), while the S r ,t have two linearly 
independent null normals. We denote the one pointing in the "outward" direction by L, and the one pointing in 
the "inwards" direction by L. The vectorfields L and L, which are defined on M/0, are unique up to multiplication 
by a scalar function. We choose the normalization so that they have the following components relative to our 
inertial coordinate system: 

(3.4.3a) U 1 = (1, -w\ -oj 2 , -w 3 ), 

(3.4.3b) W = (l,u 1 ,w 2 ,u 3 ), 

where w J = x 3 jr. With d r d = \x a d a denoting the radial vectorfield, L, L can be expressed as 



(3.4.4a) L = dt- d, 

(3.4.4b) L = d t + d,, 



We remark that beginning in Section 6.1 L and L will play a key role in the null decomposition of the MBI system. 
We now recall the definitions of the first fundamental forms of E$ and of S r t- 

Definition 3.4.2. The first fundamental forms of Hf,S r ,t are the Riemannian metrics on Y>t,S r j respectively 
induced by the spacetime metric g. In an arbitrary local coordinate system, g,$ can be expressed as follows: 

(3-4.5) g^ = f g^ + (T {0) ) ^(T {0) ) u , 

(3-4.6) ^ = f g^ + \{L^L U + L^Lj). 
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We remark that the tensors j" = 5^ + (T^) fl (T^) 1 ' and <j)^ = 5^ + ^(L^L 11 ■\-L fl L v j orthogonally project onto 

Et and SV,t respectively. Furthermore, the volume forms of g and fi, which we respectively denote by t VK \ and ^, 
can be expressed as follows relative to an arbitrary coordinate system: 

(3-4.7) € ukX = e^xT^, 

(3.4.8) ^ = l -e^ UKX L K L x . 

Definition 3.4.3. Let U be a type ( n ) spacetime tensor. We say that U is tangent to the time slices Ej if 

(3.4.9) ET.. ll Ur " Uri = g pl -g Tim gJf 1 -gJf 1 U~ „ v ^ v ™. 

Equivalently, U is tangent to the Ej if and only if any contraction of U with T( ) results in 0. 
Similarly, we say that U is tangent to the spheres S r> t if 

(3.4.10) U uv " Vn =h pl ■■■ <h pm i>^ ••• 6~ U1 U~ ~ Vl '" Vn . 

Equivalently, U is tangent to the spheres S r t if and only if any contraction of U with either L or L results in 0. 

We also recall the following relationships between the Levi-Civita connections V,^ corresponding to g,</i and 
the Levi-Civita connection V corresponding to g, which are valid for any tensor U of type (^) tangent to the 
T, t ,Sr t t respectively: 

(3-4.12) fx =tftf - ifctf - VjU^J 1 ^- 



As in ( |3.2.8 ), throughout the article, we use the notation 



(3.4.13) V x = f X K V K , (if X is tangent to E t ), 

(3.4.14) fx d =X K f K , (if X is tangent to S r ,t). 

We recall the definitions of the second fundamental form of the Et, and null second fundamental forms of the 
S r ,t- 

Definition 3.4.4. The second fundamental form of the hypersurface E^ is defined to be the tensorfield 

(3.4.15) v M (r (0) )„. 

The null second fundamental forms of the S r t are defined to be the following pair of tensorfields: 

(3.4.16) V^L U , V M L„. 

In the next lemma, we illustrate one of key properties of the second fundamental forms. 

Lemma 3.4.1. The second fundamental form V^(Tf \) v is a symmetric type (°) tensorfield that is tangent to the 
time slices Similarly, the null second fundamental forms V^L U , V ^L u are symmetric type (°) tensorfields that 
are tangent to the spheres S Tt t- 

Proof. The fact that V^L U , V u,L v are tangent to the S r t follows from contracting them with the vectors L, L, which 
form a basis for the orthogonal complement (in M) of the tangent space of S r ^ : and using ( |6.1.2a ) - (6.1.2b). For 



the symmetry property, let X, Y be vectorfields tangent that are tangent to S r j- Then [X, Y] is also tangent to 



S r> t- Therefore, using the fact that V u,L v is tangent to the S r t, (3.2.4), and (3.2.5), we deduce that 
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(3.4.17) 



XfyVv^v = g( VxL, Y) = Vx g(L, Y) -g(L, V X Y) 



= -g(L,V Y X)-g(L,[X,Y]) 



= -S/Yg(L,X)+g(S/ Y L,X) 
= g(V Y L,X) = Y^X u V l ,L u . 

The proofs for V^L U and V /i (T( ))v are similar. 

Remark 3.4.1. Lemma 7.0. 11| provides very simple expressions for the null second fundamental forms. 
Remark 3.4.2. By Lemma 



□ 



3.4.1 



use the abbreviations ^ L v ^g^ VjjL 1 ' and ^ U V-p,L u , which should cause no confusion. 

To conclude this section, we recall the following basic facts concerning the metrics g and g . 



we have that g^ \Jf L L v = V ^L v ', and similarly for L. Therefore, we sometimes 



dcf , 



Lemma 3.4.2. Let g and $ be the first fundamental forms of g defined in Definition 3.4-2. Let V, be their 



corresponding Levi-Civita connections, as defined in (3.4.11), (3.4.12) respectively. Then 



(3.4.18) 
(3.4.19) 



Y A 5 



0. 
0. 



(A,^ = 0,l,2,3), 
(A,m,^ = 0,1,2,3). 



Proof. Lemma 3.4.2 follows from the expressions (3.4.11) - (3.4.12) and Lemma 3.4.1 



□ 



4. The Maxwell-Born-Infeld System 



In this section, we first discuss the equations of motion for a generic covariant theory of classical electromagnetism 
that is derivable from a Lagrangian. We then introduce the Maxwell-Born-Infeld Lagrangian and derive several 



versions of its equations of motion. The final version, namely equations (4.3.2a) 
throughout the remainder of the article. 



(4.3.2b), will be the one we use 



4.1. The Lagrangian formulation of nonlinear electromagnetism. In this section, we recall some facts from 
classical nonlinear electromagnetic field theory in a Lorentzian manifold (M, g) of signature (-,+,+,+). We restrict 
our attention to theories of nonlinear electromagnetism derivable from a Lagrangian ££ . The fundamental quantity 
in such a theory is the Faraday tensor J 7 ^, a two- form (i.e., an anti-symmetric tensorfield) that is postulated to 
be closed: 



(4.1.1) dT = 0, 

where d denotes the exterior derivative operator. This equation, which is the first of two equations that will define 
a particular nonlinear theory, is known as the Faraday-Maxwell law. In local coordinates, it can be expressed in 
the following two ways 

(4.1.2a) <9 [A J>] = 0, (X,fi,v = 0, 1,2,3) 

(4.1.2b) V [A ^]=0, (A, fi t v = Q,l, 2, 3), 

where [■■■] denotes antisymmetrization. 
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In any covariant theory of classical electromagnetism, *«Sf is a scalar- valued function of the two invariants of T, 
which we denote by and ^( 2 ); i.e., *JSf = *«Sf (^m [J 7 ], ^ 2 ) t^ 7 ])- They can be expressed in the following ways: 



(4.1.3a) 
(4.1.3b) 



*( 2 )=* (a) F] "^(p- l )""(/r l ) A "^,x*^ 



r l r M (5 -l)A- 
-1 -1 \ Af - 



-*(T**T) = 



l-El' 



1 

8~ 



1 



XT aT) = E K B K , 



where a denotes the wedge product, and E, B are the electromagnetic one-forms defined in Section 6.8 As we 
will discuss in Section 6.4 the invariants and ^ 2 )> viewed as quadratic forms in T, have a special algebraic 
structure. More specifically, we will see that from the point of view of the decay estimates of Proposition 11.0.9 



the worst possible quadratic terms are absent from and *7(2)- This is the one of the fundamental reasons that 
small-data solutions to the MBI system have the same decay properties as solutions to the linear Maxwell-Maxwell 
equations. 

We now introduce the Maxwell tensor A4, a two tensor whose Hodge dual *A4 is defined by 



(4.1.4) 

Furthermore, we postulate that A4 is closed: 
(4.1.5) 



*M 



de£ 



dT L 



dM = 0. 



Taken together, (4.1.1) and (4.1.5) are the equations of motion for the theory arising from the Lagrangian Jzf '. We 



remark that (4.1.1) and (4.1.5) are respectively equivalent to 



(4.1.6a) 
(4.1.6b) 



= 0, 



= 0, 



(f = 0, 1,2,3), 
(u = 0, 1,2,3). 



Equations (4.1.6a) are sometimes referred to as the Bianchi identities. We furthermore remark that the solutions 



to (4.1.1), (4.1.5) are exactly the stationary points (under closed variations dT - with support contained in 



compact subsets <£) of the action functional 

(4.1.7) AdT] d = cf f £sM ^ {1) [T],^ 2) [T])d^ 



9i 



and that (4.1.6b) are the Euler-Lagrange equations of *Jzf. In the above formula, d[i g d = f Idetgl 1 / 2 ^ 4 ^ is the measure 



associated to the spacetime volume form (3.3.1) 



The Euler-Lagrange equations (4.1.6b) can be written in the following form: 



(4.1.8) 
where 



h^ KX V„T KX = 



(i/ = 0, 1,2,3), 



(4.1.9) W" vkX = 

The tensorfield h^ UKX , which has the properties 



1 d 2 *££ 

2 dT au dT K 



(4.1.10a) 
(4.1.10b) 
(4.1.10c) 

is of fundamental importance throughout this article. As is explained in Section [8j its algebraic and geometric 
properties are intimately related to the hyperbolic nature of the MBI system. In particular, the symmetry prop- 



erties (4.1.10a) - (4.1.10c) are needed to ensure that the canonical stress tensor, which is defined in Section 8.3 



has lower-order divergence. 
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We state as a lemma the following identities, which will be used for various computations. We leave the proof 
as an exercise for the reader. 

Lemma 4.1.1. The following identities hold: 

(4.1.Ha) f^l^MGT 1 )-, 

(4.1.11b) d (g ~ 1)WA =-(g- 1 T»(g- 1 ) Xl/ , 

(4.1.11c) = \detT\\det g]' 1 , 

(4.1. lid) (g^r^^ux - (g^r^^Fvx = 

(4.1. He) (g^V^^uX = \ 2 )9^ 

(4.1.11f) ^ = -g K x^ K F u \ 

og^u 

9^(2) 1 

(4.1. llh) ^ = 2^-^+^T^, 

%(2) 

(4.1.111) = 27*", 

(4-l.Hj) = *P iV , 

(4.i.nk) = (g-'r^g- 1 )^ - (g-y^r, 

(4.1.111) V— = e^ KA . 

4.2. Derivation of the MBI equations. The Lagrangian for the MBI model is 

(4.2.1) *J? (MBI) d = ef ^ - ^(1 + (3%) LF] - P 8 ^ 2) m) 1/2 = ^ - ^(det s ( 5 + F)f'\ 

where |3 > denotes Bom's "aether" constant. For the remainder of the article, we set (3 = 1 for simplicity] 
however, the analysis in the case |3 + 1 easily reduces to the case (3 = 1 by making change of variable T = $ 2 J- in 
the equations. For future use, we introduce the abbreviation 

(4-2.2) W)= f ( 1 + V)-^(2)) 1/2 > 

which implies that 

(4.2.3) *££ (MBI) = 1 ~ i(MBI) ■ 



Using definition (4.1.4) and Lemma 4.1.1, we compute that in the MBI model, !M M can be expressed as follows: 



(4.2.4) = -r {MB i } {^ v - l^x^rn = -q MBI) {f - HtfJ*"). 



Taking the Hodge dual of (4.2.4), we have that 



(4.2.5) M^ = e7 Mm) (*^ + ^ i2) ^). 



From (4.2.4), it follows that the Euler-Lagrange equations (4.1.6b) for the MBI model are 
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(4.2.6) - \*^ v V^ K x*F KX ) - \q 2 MBI) {T^ - \t k xT kX ^ u )v^\t kX ^ x - 1(^ kA ^ a ) 2 ) = 0. 

Furthermore, it follows from (4.2.6) that the tensorfield h^ VKlX from ( 4.1.9[ ) can be expressed as 



(4.2.7) hr* = \li$ i Bi)[(9- 1 )' m (9- 1 r x - (g-'YHg- 1 )™] - t ( 3 MBI) ^F KX + ^qUi)^^ + *j* v f kX ) 

~~ ^(MBJ) + *(2) t (MBI) ) * J~ ~ c (MBiy(2) e p 

4.3. H^ vkX and the working version of the MBI equations. To ease the calculations, it is convenient to 
perform two simple modifications of the tensorfield from (4.2.7) obtaining a new tensorfield; the modi- 

fications will not alter the set of solutions to the MBI system. First, we drop the ^lj B j^(2) eflUKX term from 



(4.2.7). This is permissible because its contribution to the Euler-Lagrange equations is ^mbi^C 2 ) 6 ^^^ v-F k \ = 
^(MBI)^(^)^ )j*^ U ~ 0> on account of equation (4.1.6a). Second, we multiply the remaining terms in (4.2.7) by 
i(MBi)- We denote the resulting tensorfield by H^ vkX . Furthermore, it is convenient to split H^ VK ^ into a main 
term, which coincides with the tensorfield in the case of the Maxwell-Maxwell equations, and a quadratic error 



term, which we denote by H A 



(4.3.1a) 
(4.3.1b) 



jjIivkX def 



1 



livnX 



)-h(g- 1 r K (g- 1 y x -( 9 - 1 r x (g- 1 r] + H, 



[ivkX 



H 



1±vkX def 



■t. 



-2 

(MBI) 



rnupnx + ^ 2) r^ BI) [^ v *T KX + *t^t kX ) - (l + ^2) r (MBi)Y^ U ^ K> 



It follows that the system (4.1.6a), (4.1.6b), (4.2.4) is equivalent to the following version of the MBI system: 



(4.3.2a) 
(4.3.2b) 



0, 



H^ kX X7^T k x = 



(A, fx,v = 0,1,2,3), 
= 0,1,2,3). 



5. CONFORMAL KILLING FIELDS AND MODIFIED LlE DERIVATIVES 



In this section, we recall the definition of conformal Killing fields. This collection of vectorfields, which has 
the structure of a Lie algebra under the Lie bracket operator X, Y [X, Y] , comprises the generators of the 
conformal symmetries of the spacetime (M,g). We focus on the case of Minkowski space, which has the maximum 
possible number of generators (15). In particular, we introduce several subsets of the Minkowski conformal Killing 
fields, each of which will play a role throughout the remainder of the article. More specifically, they appear in the 
definitions of the norms and energies (see Section [9]) that are used during our global existence argument. Finally, 
for a special collection of Minkowski conformal Killing fields Z we define modified Lie derivatives Cz-, which are 
equal to ordinary Lie derivatives plus a scalar multiple of the identity. This definition is justified by Lemma 



7.0.4 



which shows that the operator Cz has favorable commutation properties with the linear Maxwell-Maxwell 
equation [(<f ^ '(if 1 )"* - (<f ^Or 1 )"*] V^ kA = 0. 



Definition 5.0.1. A Killing field of the metric g^ v is a vectorfield X such that 
(5.0.3) (x) vr^ = 0, 

while a conformal Killing field X satisfies 

(5.0.4) Wvr^ = <f> xgiu , 
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for some scalar- valued function 4>x(t,x)- In the above formulas, the deformation tensor ^ '7T^ U is defined in (3.2.12 ) 



The conformal Killing fields of the Minkowski metric are generated by the following 15 vectorfields (see e.g. 
IChrOSj ): 

(1) the four translations Tr^, (fi = 0, 1,2,3), 

(2) the three rotations Qfjty, (1 < j < k < 3), 

(3) the three Lorentz boosts ^(oj)> (j = 1)2,3), 

(4) the scaling vectorfield S, 

(5) the four acceleration vectorfields K^\, (// = 0, 1,2,3). 

Relative to the inertial coordinate system {x At } At= o,i,2,3, the above vectorfields can be expressed as 

(5.0.5a) T (fi) = a M) 

(5.0.5b) O(^) = x^dy - x v d^ 

(5.0.5c) S = x K d K , 

(5.0.5d) = -2 Xfl S + 5K AX K x A a M . 

When working in our fixed inertial coordinate system, we use the notation T( ) = dt interchangeably. In this article, 

we will primarily make use of the vectorfields in (1) - (4), together with K d = f +T^, which has the following 
components relative to the inertial coordinate system: 

(5.0.6a) K° = 1 + t 2 + (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , 

(5.0.6b) K 3 =2tx j , = 1,2,3). 

We remark that the translations, rotations, and boosts are Killing fields, while relative to the inertial coordinate 
system, we have that 

(5.0.7a) (5) 7r^ = 2o 



(5.0.7b) ^WK^ = -4x x g, 

(5.0.7c) = 4tg 



(IV 



The subset ', which consists of all the translation^} the subset O, which consists of the rotations, and the 
subset Z, which consists of all generators except for the accelerations, and which have cardinalities 4, 3, and 11 
respectively, will play a distinguished role throughout this article: 

(5.0.8a) ^ d = f {T (M) }o<^<3, 

(5.0.8b) o d ={n (12) ,n i23) ,n {13) }, 

(5-0. 8c) Z d = f {T( M ) , , S} <^<v<3. 

We denote the Lie algebras generated by ^,0, and Z by T, O, and Z respectively. Note that for each vectorfield 
Z e Z, there is a constant cz such that 

(5.0.9a) CzQiiv = czgfiu, 

(5.0.9b) C z {g- l Y v = -c z {g- x Y v . 



Also note that on the left-hand side of (5.0.9b), the indices of g are raised before differentiation occurs. 



It will be convenient for us to work with modified Lie derivatives^' 1 ] Cz- 



This is not to be confused with the subset T of frame field vectors, which is defined in (6.f .4a l 



Note that these are not the same modified Lie derivatives that appear in |BZ09] . [CK93 



. and ZipOO . 



The Nonlinear Stability of the Trivial Solution to the Maxwell-Born-lnfeld System 
Definition 5.0.2. For each vectorfield Z e Z, we define the modified Lie derivative Cz by 



25 



(5.0.10) 



del' 



where cz denotes the constant from (5.0.9a). 



The crucial feature of the above definition is captured by Lemma 7.0.4 below, which shows that for each Z e Z, 
the operator Cz can be commuted through the linear Maxwell-Maxwell equation 

[Gr 1 )^- 1 )^ - (g'^ V)" ^^ = 0, resulting in the ide ntity [{g' 1 )^ {g' 1 )^ - (<T 1 )' iA Gr 1 ) 1 " , ]v ^ z T kX = 0. 
As is shown in Lemmas 10.0.5 and the proof of Lemma 12.0.2, a similar result (involving nonlinear error terms) 
also holds for the MBI equation ( |4.3.2b ). 

We now introduce some notation that will allow us to compactly express iterated derivatives. If A is one of the 



sets from (5.0.8a) - (5.0.8c), then we label the vectorfields in A as Z Ll , Z im , where m is the cardinality of A. 



Then for any multi-index / = (ti,---, ife) of length k, where each l% e {1,2, ■■■,m}, we define 



Definition 5.0.3. 






(5.0.11a) 


pi def n 

L A = Lzn ° ■ 


■■ o C Z <- k 


(5.0.11b) 


pi def p 

L A = Lzn o • 


■■ o C Z l k 


(5.0.11c) 


^1 def „ 


■■ o Vz'fe 


(5.0.11d) 


to = fz>i °- 


• Vz'fc , 



etc. 

Under this convention, the Leibniz rule can be written as 

(5.0.12) C I Z (UV)= £ (^zU)(C z V), 

h+h=l 

etc., where by a sum over I\ + I2 = I, we mean a sum over all order preserving partitions of the index / into two 
multi-indices; i.e., if I = (ti, ik), then I\ = , ij a ), I2 = (^ a+1 , ii k ), where ii, i& is any re-ordering of the 
integers 1, k such that i\ < ■■■ < i a , and i a+ \ < ■■■ < ik- 

We end this section with the following lemmas, which provide expressions for the commutators and 
[Vx,*] acting on two-forms. 

Lemma 5.0.1. If X is a vectorfield, and T is a two-form, then 



(5.0.13) 



Proof. Lemma 5.0.1 follows easily from definition 3.3.4 and the fact V x^uukX = 0, which is a simple consequence 



of (3.3.3). 



□ 



Lemma 5.0.2. [CK90J Eqn. 3.25] If X is a vectorfield, and T is a two-form, then 



(5.0.14) 



Proof. The relation (5.0.14) follows from the using the expressions (3.3.4) and (3.2.10) for *J- and Cx, together 



with the well-known identities 



(5.0.15) 

(5.0.16) 
(5.0.17) 



£x€k\[iu 



>^ ^ f3 e K.Xp,u> 



XjTU 



cA 



to express Cx*Fuv in terms of the components of *T . We leave the details to the reader. 



□ 
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We are particularly interested in the case that X is a conformal Killing field. Under this assumption, we have 
the following simple corollary of the lemma. 

Corollary 5.0.3. [CK90, Corollary of Proposition 3.3] If X is a conformal Killing field and T is a two-form, then 



(5.0.18) 



Proof. Corollary 5.0.3 follows from (5.0.4) and Lemma 5.0.2 



□ 



6. Tensorial Decompositions 

In this section, we have three main goals. First, we will introduce the null frame decomposition of tensorfields. 
Related to this decomposition is the notion of a null form Q(-, •), which is a quadratic form that acts on a pair of 
type ( 2 ) tensors, and that has a special algebraic property: the "worst" possible quadratic combinations, from the 
point of view of the null decomposition, are absent. Since the two invariants of the Faraday tensor are each 
multiples of a corresponding null form Qu\(J-, F), the net effect is that every nonlinear term in the Euler-Lagrange 
equation (4.2.6) of the MBI system can be expressed as functions of null forms in T,X/T. Later in the article, 



with the help of the additional null form structure present in the the expression (9.2.2) for the divergence of our 
energy currents, together with Proposition 11.0.9, we will be able to prove sharp decay estimates for the null 



components of solutions J- to the MBI system. Next in this section, we decompose the MBI system relative to the 
null frame. From the point of view of proving the sharpest possible decay estimates, the most important equation 
is (6.7.7b), which shows that the "worst" derivative of the a null component of T can be expressed in terms of 
"good" derivatives of other null components of plus cubic error terms. Finally, we introduce electromagnetic 
decompositions of T and Ai, where A4 is the Maxwell tensor from (4.2.5 ). These decompositions will be useful for 



proving various identities and inequalities concerning and for expressing the smallness condition in our global 
existence theorem directly in terms of the data (B, D), which are one-forms inherent to the Cauchy hypersurface Sq. 



6.1. The Null frame. Before proceeding, we introduce the subsets Cg,C s of Minkowski space, which will play 
a role in the sequel. 

Definition 6.1.1. In our fixed inertial coordinate system (t,y), we define the outgoing Minkowski null cones C*, 
and ingoing Minkowski null cones C~, as follows: 

(6.1.1a) c;^{(T,y)\\y\-T = q}, 

(6.1.1b) C;^{(T,y)\\y\ + T = s}. 

In the above formulas, y d = (y 1 , y 2 , y 3 ) , and \y\ d = f ((y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 ) 1 ^ 2 - 

A Minkowski null frame is a locally defined collection of four vectorfields on M/0 that we will denote by 
L, L,e%,e2, and that have the following properties: 

• At each point p, the set {L, L, e\,e-i\ spans T p M 2 M. 

• For each seK, and for all nonzero p e C~ , L\ p is a future-directed, ingoing null geodesic vectorfield tangent 
to C~. 

• For each gel, and for all nonzero p e Cg, L\ p is a future-directed, outgoing null geodesic vectorfield tangent 
to C+. 

• For all t 6 M, for all r > 0, r = \x\, and all points p e S r< t, L\ p and L p are normal to S r> t- 

• g(L,L) = -2. 

• At each point p e S r< t, e\\ p and e2\ p are tangent to SV,f. 

• 9{eA,e B ) = Sab- 



The Nonlinear Stability of the Trivial Solution to the Maxwell-Born-lnfeld System 

27 

In formulas (3.4.3b) - (3.4.3a), we provided the components of L,L relative to our inertial coordinate system. 
Relative to an arbitrary coordinate system, the aforementioned properties of L, L, e\, and e 2 can be expressed as 
follows: 



(6.1.2a) 


\7 L L = \JlL 


(6.1.2b) 


VlL = VlL 


(6.1.2c) 


L K L K = -2, 


(6.1. 2d) 


e A^n = e A L K 


(6.1. 2e) 


9K\e-A e> B = &AB- 



Additionally, it follows from (3.2.4) and (6.1.2b) that L and L commute as vectorfields: 
(6.1.3) [L,L] = 0. 



In the analysis that will follow, we will see that the decay rates of the null components (see Section 6.3) of 
J- will be distinguished according to the kinds of contractions of J- taken against L, L, e%, and e 2 . With these 
considerations in mind, we introduce the following sets of vectorfields: 



£ d = cf {L}, 



T = {L,ei,e 2 }, 



U A = {L,L,ei,e 2 }. 



(6.1.4a) 

We will often need to measure the size of the contractions of various tensors and their covariant derivatives 
against vectors belonging to the sets C,T,U. This motivates the next definition. 

Definition 6.1.2. If V, W denote any two of the above sets, and F is an arbitrary type (2) tensor, then we define 
the following pointwise seminorms: 



(6.1.5a) 
(6.1.5b) 



\F\ 



|VF| 



vw 



del' 



E 

VeV,WeW 



del' 



vw 



\v k w x f kX \, 



Y, \v k w x u^v^f kX \ 

Uell,VeV,WeW 



Observe that if F is a two tensor, then \F\ w \F\uu, where \F\ is defined in (3.1.7). 



6.2. Null frame decomposition of a tensorfield. For an arbitrary vectorfield X and frame vector U eW, we 
define 



(6.2.1) 



X V d ^X K U K , where X^ d =g^X K . 



The components Xjj are known as the null components of X. In the sequel, we will abbreviate 



(6.2.2) 



V d i: f V 



It follows from (6.2.1) that 



„ def „ 

Va = V eA , etc. 



X = X K d K = X L L + X L L + X A e A , 



X = — Xt, , 

2 - 



x- = — X L , 
2 



x = X A . 



(6.2.3a) 
(6.2.3b) 

Furthermore, it is easy to check that 

(6.2.4) g(X,Y) = X K Y K = ~X L Y L - \xlY L + S AB X A Y B . 

The above null decomposition of a vectorfield generalizes in the obvious way to higher order tensorfields. In the 
next section, we provide a detailed version of the null decomposition of two-forms since they are the fundamental 
unknowns in any classical theory of electromagnetism. 
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6.3. The detailed null decomposition of a two-form. 

Definition 6.3.1. Given any two-form J-, we define its null components a,a,p,o~ as follows: 

(6.3.1a) = ^ T V \F^ , 

(6.3.1b) a M T uX L X , 

(6.3.1c) p^ tl -F KX L^L\ 

(6.3.1d) cr d = ef ifV KA . 

It is a simple exercise to check that a, a are tangent to the spheres S r j • 

(6.3.2a) a K L K = 0, q k L k = 0, 

(6.3.2b) ct K L K = 0, a K L K = 0. 

Furthermore, relative to the null frame {L, L,ei,e2}, we have that 

(6.3.3a) a A = J^al, 

(6.3.3b) «A = -^al, 

(6.3.3c) P=^ll, 



(6.3.3d) a = Fi 



12- 



In terms of the norms introduced in Definition 6.1.2 it follows that 



(6.3.4a) \T\ m |.F| WW « |a| + |a| + |p| + 

(6.3.4b) \T\cu - |a| + \p\, 

(6.3.4c) l^lrr * M + H- 

The null components of *T can be expressed in terms of the above null components of T . Denoting the null 
components 1 ^ of *T by ®a, ®a, @ p, ®cr, we leave it as a simple exercise to the reader to check that 



(6.3.5a) @ a A = -a £ba, 

(6.3.5b) ®oa = a B ^BA, 
(6.3.5c) @ P = cr, 

(6.3.5d) ®<r = -p. 

6.4. Null forms. 

Definition 6.4.1. Let F, G be any pair of type (g) tensors. We define the null forms Q(i)(-, •), Q(2)('i ■) as follows: 

(6.4.1a) Q {1) (F,G)^F kX G kX , 
(6.4.1b) Q (2) (F,G)^*F kX G kX . 
It is easy to check that Q (i) (F,G) = Q (i) (G,F) for i = 1,2. 



18 



We use the symbol © in order to avoid confusion with the Hodge dual; i.e., it is not true that *(^[-^]) = M*-? 7 ]- 
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Remark 6.4.1. Observe that the two invariants (|4.1.3a|) - (|4.1.3b|) of the Faraday tensor J- are 
(6.4.2) 



Based on the above remark, it is clear that we will be primarily interested in the case in which F, G are two- 
forms. The next lemma describes the fundamental algebraic properties of null forms; i.e., the absence of the "worst 
possible" quadratic terms. 



Lemma 6.4.1. Let F,Q be a pair of two-forms, and let Q(j)(-,-) be one the null forms defined in (6.4.1a 



( |6.4.1b) , and let \ ■ |vw be the norms defined in ( 6.1.5a[ ) Then for i = 1,2, the following pointwise estimate holds: 

(6.4.3) Q(i)GF,S) S \H0\cu + \F\cu\G\ + \?\tt\Q\tt- 

Proof. Let afj-"], crLF], p[T] and a[G], a[Q], a[Q], p[Q] denote the null components of T and Q respectively. 

Then using the relations (6.3.5a) - (6.3.5d) (for the case i = 2), we compute that 

(6.4.4) F kX Q rX = -5 AB a A [T]a B [G] - 5 AB a A [G]a B [F] - 2p[T]p[G] + 2a[T]a[G], 

(6.4.5) *T kX G k x = i AL *a A [F]a B [g] + { AB a A [Q]a B [T] - 2a[T]p[Q] - 2p[T]a[G]. 



from which (6.4.3) follows. 



□ 



6.5. Intrinsic divergence and curl, and the cross product. In this section, we recall the definitions of the 
intrinsic divergence and curl of vectorfields U that are tangent to the submanifolds St and S r t- 

Definition 6.5.1. If U is a vectorfield tangent to St, then its intrinsic divergence and curls are defined relative 
to an arbitrary spacetime coordinate system as follows: 



(6.5.1) 
(6.5.2) 



<tivU d =V K U K = g A V A C/ K , 

— K 

( curl U) v d - ~ VK " ttX 



(!/ = (), 1,2,3), 



where the volume form e VK \ of Ef is defined in (3.4.7). Relative to the Euclidean coordinate system x = {x 1 } x 2 ,x^) 
on Ht, we have that 



(6.5.3) 
(6.5.4) 



div U = V a U a , 
( curl Uy d =^ a h V U\ 



0' = 1,2,3). 



If U is a vectorfield tangent to the spheres S r ,t, then its intrinsic divergence and curls are defined to be the 
following scalar quantities: 



(6.5.5) 
(6.5.6) 



dfvU^f K U K ^ X V X U K , 



curll^Av^ 



where ^n U , the volume form S r j, is defined in (3.4.8). Relative to a null frame, we have that 



(6.5.7) dfrU d =6 AB f A U B , 

(6.5.8) cir\U d =tABfAU B . 

Note that in the above definitions, contractions against the frame vectors e\,e2 is taken after covariant differ- 
entiation; e.g., fA U B d = e A e B f K U\. 

Definition 6.5.2. If U and V are vectors tangent to Sj, then we define their cross product, which is also tangent 
to St, as follows, relative to the Euclidean coordinate system x = (x 1 ,x 2 ,x 3 ) on St : 



(6.5.9) 



(U x vy d =e> ab u a v b . 
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6.6. The null components of K and VK. In this section, we provide the null components of the vectorfield 



K, which is defined in (5.0.6a), and its covariant derivative VK. K is central to our global existence argument, 



because it is a fundamental ingredient in the energies we construct; see (9.0.5) and (9.0.7). We don't provide any 



proofs in this section, but instead leave the simple computations as an exercise for the reader. 



We first note that using definitions ( |3.4.4b[ ) and ( |3.4.4a[ ), K can be expressed as follows: 
(6.6.1) 



1 



K = -{(l + s 2 )L + (l + q 2 )L}. 



From (6.1.2c), (6.1. 2d), and (6.6.1), it easily follows that 

(6.6.2a) 
(6.6.2b) 
(6.6.2c) 



Finally, we compute that 



(1 V), 

2^ 



K L = - 

K k =-{l + s z ), 
K A = 0. 



(6.6.3a) 




(6.6.3b) 


VLKL d =L K L X S7 K K x = 4q, 


(6.6.3c) 


\7 A K B d =e K A e x B \7 K Kx = 2t5 AB , 


(6.6.3d) 


V L K L d = L K L X X7 K K X = 0, 


(6.6.3e) 


S7 L K A d = L K e x A S7 K K x = 0, 


(6.6.3f) 


VlKl A = L k L x X7 k K x = 0, 


(6.6.3g) 


X7 L KA d = L K e x A X7j( x = 0, 


(6.6.3h) 


V^L= f elL A V K K A = 0, 


(6.6.3i) 


v^L d = of 4L A v K K A = o. 



6.7. The null decomposition of the MBI system. In this section, we decompose the MBI system into equa- 



tions for the null components of T . We begin by noting that after simple computations, the MBI system (1.0.1a) 



(1.0.2) can be expressed in the following equivalent form: 



(6.7.1a) Vx^nv + V^a + V„J> = 0, 
(6.7.1b) VaT^ + V/^a + V/^ 

" 2^A/B/){( V Al) " 2 ' '(2) V A2))*-7>i, + (V^(i) - 2^( 2 )V^(2))*^A + (Vi/fyl) - 2^( 2 )V ! y^(2))*^ r A M } 

~~ 2^(MB/)' ? ( 2 ){( Va ^( 1 ) ~~ 2 ^(2)Va^(2))-7>" + ( V^(l) - 2^(2)V /1 ^(2))^A + (Vi^(l) - 2^(2)V ! y^( 2 ))^ r A M } 
+ (Va^(2))^/» + (V^( 2 ))^A + (Vr/^2))^ = 0. 

In our calculations below, we will make use of the identities 



(6.7.2) 



VaLl = -r 1 e A , 



VaL = r l eA, 



which can be directly calculated in our inertial coordinate system using (3.4.3a) - (3.4.3b). We will also make use 
of the identity 
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(6.7.3) 



fA e B = Va&b + ^g(VAeB,L)L + ^g(\7 A eB,L)L 
= Vacs - -g(e B ,V aL)L - -g(e B , V aL)L 
= Va&b + ]^~ 1 8ab(.L - L), 



which follows from (3.4.12) and d6.7.2) 



Furthermore, if U is a type (° J tensorfield, and Xr^, 1 < i < m and Y are vectorfields, then by (3.2.6), we have 
that 



(6.7.4) Vy(f/(X( 1 ),X(2), •••,X( m ))) = (VyJ7)(X( 1 ),X( 2 ),---,X( m )) + C/(VyX( 1 ),X( 2 ),-">^"(m)) 

+ ••• + {7(X( 1 ),X(2),---, VyX( m )). 

Similarly, if [/ is tangent to the spheres Sy^, and e B(1) , e B(m) e {ei,e 2 }, then 
(6.7.5) 

fA {U B(iy -B (m) ) d =ie A {U{e B(1) ,e B(2y -,e B(m) ) 

= (^A ^)(es (1) ,es C2) ,-,e B(m) ) + [7(^4 e B(1) , e B(2) e B(m) ) + ■■■ + C/(e B(1) , e B(2) ,-, e B(ra) ). 



Applying (6.7.4) and (6.7.5) to J 7 , and using (6.7.2) plus (6.7.3), we compute that the identities contained in the 



following lemma hold. 

Lemma 6.7.1. [CK90, pg. 161] Let T be a two-form, and let a, a, p, a be its null components as defined in 



Definition 6.3.1. Then the following identities hold: 



(6.7.6a) 




(6.7.6b) 


^aFbl 


(6.7.6c) 


^a*Fbl 


(6.7.6d) 


Va*Fbl 


(6.7.6e) 


\vaFll 


(6.7.6f) 




(6.7.6g) 


^JaFbc 



.-1/ 



„-!/ 



1 _i 
-r L \ 
2 



\va*Fll =^a o-+ -r 1 (-£ B Aa B + £ba«b), 



VaFbc = £bc{ fAcr+]-r 1 (~iDAa D + ^daod)}- 



Note that in all of the above expressions, contractions are taken after differentiating; e.g., Va^bl = e^sL^^ ^J~k,\- 



□ 
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Lemma 6.7.2. The MBI system can be decomposed into principal terms and "cubic error terms" as follows, 
relative to a Minkowski null frame: 



(6.7.7a) 

- ^(MM){(V^(l) - 2^ (2) VL^(2))aA - HabUb *(1) - 2^(2) f B *( 2 ))<7 + (V^(i) - 2* (2) Vl^(2) )«a} 



^(2){(Vz>(i) - 2^( 2 )Vi^(2))^ J 4BQ!B - Hab($B ^(l) " 2^( 2 ) ^(2))p~ (Vi^(i) " 2^( 2 ) Vl^(2) )/ylBa B J 



(6.7.7b) 



_ 7^(MB/){(Vi^(l) - 2^( 2 )Vl>(2)W - Hab{17b ^(1) - 2^( 2 ) ^(2))^ + (Vl^(i) - 2 < 7( 2 )Vz^(2))au} 



~ 7^(MB/)^(2){(Vl^(1) - 2^( 2 )Vx^(2))^ABQ!B - Hab($ B ^(1) - 2^( 2 ) fB ^(2))p~ (Vi^(l) " V/^(2) 
+ ^{( V ^(2))^B«B - HabUb ^(2))P- (^L^(2))^ABO_ B j = 0, 

(6.7.7c) 

-e#w a - Vlp + 2r _1 p - ^ fB/) { " (Vz>(i) - % (2) Vl>( 3 ) )p - ( - 2* (2) *( 2 ) )0s} 

- ^(mb/)^(2){(VlV) - 2^ (2) Vl^( 2 ))^ + &ib(?a V) - 2^ (2) ^ ^(2))«b} + {(Vl>( 2 )V + Unit A ^(2))«b} = 0, 
(6.7.7d) 

cifrl a + Vl<j - 2r~ l o = 0, 
(6.7.7e) 

d{va-V L p-2r- 1 p--r ( l [BI) [ - (Vx^(i) - 2* (2) Vl^( 2 ))p + <5 AB (Va - 2^ (2) f A ^(2))a B } 

~ ^(MB/)^(2){(ViV) - 2^ (2) Vl^(2))^ + UbHa - 2-7(2) f A fyj))***} + {(Vz^(2))^ + IabW 'a ^(2))«b) = 0, 

(6.7.7f) 

c?/H a + Vlc + 2r _1 (j = 0, 



where the differential operators dfv andcdrl are defined in (6.5.7) - (6.5.8). 

Remark 6.7.1. If we discard the nonlinear terms, then the resulting system is the null decomposition of the 
Maxwell-Maxwell system. 

Proof. Contract the vectors L A L M e^ against equation (6.7.1a), and then against equation (6.7.1b). Adding the 
resulting equations and using (6.7.6a) - (6.7.6g) gives (6.7.7a), while subtracting the resulting equations and using 
- (6.7.6g) gives (6.7.7b). Equations (6.7.7c) - (6.7.7d) follow from contracting L X e^e B against (6.7.1b) and 

- (6.7.7f) follow from contracting L x e^e B 

□ 



6.7.6a 



6.7.1a 



respectively, and using (6.7.6a) - (6.7.6g). Similarly, (|6.7.7e 



against (6.7.1b) and (6.7.1a) respectively, and using (6.7.6a) - (6.7.6g 



6.8. The electromagnetic decompositions of T^ u and A4^ u . In this section, we provide the familiar decom- 
position of the Faraday tensor T into the electric field E^, and the magnetic induction field B^. We also introduce 
a related decomposition of the Maxwell tensor A4 into the electric displacement D^, and the magnetic field H^. 
For computational purposes, it will also be convenient to introduce two additional one- forms, £^ and 23^, which 
are related to and B^. 
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The decompositions require the introduction of i, the interior product operator, the action of which on T is 
defined by the requirement that the following relation should hold for all vectors X, Y ■ 

(6.8.1) ixHY) = F(Y, X) = T kX Y k X x . 
In an arbitrary coordinate system, ixJ~ can be expressed as 

(6.8.2) (ix?)^V. 

Definition 6.8.1. In terms of the Faraday tensor T and the Maxwell tensor Ai (an expression for Ai in the MBI 



system is given in (4.2.51)), the one-forms E, B, D, H,£, and 53 are 



(6.8.3a) E d i { i Tw F, B d i ! -t T{0) *F, D d = f -i T(Q) *M , H d = f -l T(0) M 

(6.8.3b) £ d =i s T, 33= f i 5 ^, 



where T( ) and S are the time translation and scaling vectorfields defined in (5.0.5a) and (5.0.5c). 



In components relative to the inertial coordinate system {x /i } /J= o,i,2,3 5 we have that 
(6.8.4a) E li = F va , -*T„ , D„= -*M„ Q , H^-M^o, 

(6.8.4b) ^ = x K T^ K , <B M = x K *^ K . 

The anti-symmetry of Ai implies that E,B,D, and H are tangent to the hyperplane Yif In the inertial coordinate 
system, this is equivalent to Eq-Bq- Dq = Hq = 0. We may therefore view these four quantities as one-forms that 
are intrinsic to Ef. 

Remark 6.8.1. Our definition of B coincides with the one commonly found in the physics literature, but it has 
the opposite sign convention of the definition given in [CK90j . 

The identity 

(6.8.5) g KX X K X x ^ u = {%xF)»X v - (ixF)^ + X K (i x *T) x e KXflu 

shows that T is completely determined by ixF and ix*J~ at any spacetime point where g K \X K X x t 0. In particular, 
E and B completely determine and it can be checked that in the inertial coordinate system, 

(6.8.6a) Fjo = Ej, Tj\. = e ijfe S l , Bj = -ef^Fab- 

(6.8.6b) M j0 = -H j , M jk = e ijk D\ D <j = ^e/ 'M ab . 

Now in linear Maxwell-Maxwell theory, the relations E - D, B = H hold. In contrast, in the case of the MBI 
system, tedious computations lead to the following relations (see e.g. [BB83J, |Kie04a]): 

E+(E a B a )B 

(6.8.7a) D = 



(6.8.7b) H = 



(l + \B\2-\E\i-(E a B-y) 1/2 ' 
B - (E a B a )E 



(l + \B\ 2 - \E\ 2 - (E a B a ) 2 ) 



1/2' 



34 

Furthermore, we have that 
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(6.8.8a) 
(6.8.8b) 



E = 
H = 



D + Bx(DxB) 



(1 + \B\ 2 + \D\ 2 + \Dx £12)1/2' 

B-Dx(DxB) 
(l + \B\ 2 + \D\ 2 + \DxB\ 2 y/ 2 ' 



In the above formulas, x denotes the usual intrinsic cross product on Sj; see (6.5.9). 



Working in our inertial coordinate system, we set v = in equations (4.1.6a) - (4.1.6b) and use (6.8.4a) to 
deduce the constraint equations for B and D : 



(6.8.9a) 
(6.8.9b) 



div D = 0, 
div B = 0. 



Setting A = 0, fi = a, v = b in (6.7.1a), then contracting against €j ab and using (6.8.6a), we deduce that 



(6.8.10a) 



d t B = - curl E. 



Similarly, we use an equivalent (modulo equation (6.7.1a)) version of (6.7.1b), namely V \Ai^ u + V^.Mj,a + Vv-M-Xfj, 



0, set A = 0, fj, = a, v = b, contract against ej ab , and use (|6.8.6b|) to deduce that 



(6.8.10b) 



dfD = curl H. 



Now from (6.8.8a), (6.8.8b), (6.8.10a), and (6.8.10b), it follows that B and D satisfy the following evolution 
equations: 



(6.8.11a) 
(6.8.11b) 



dtB = - curl 
dtD = curl 



D + Bx{DxB) 



{l + \B\ 2 + \D\ 2 + \DxB\ 2 yi 2 

B-Dx(DxB) 
(l + | J B| 2 + | J D| 2 + | J Dx J B| 2 ) 1 /2 



(6.8.9a 


)-( 


6.8.9b 


) plus ( 


6.8.11a 


) - (6.8.11b 



In regards to the terminology "constraint equations" used above, we make the following remark: it follows 



that if (6.8.9a) - (6.8.9b) are satisfied along So, and if D,B are classical solutions to the evolution equations 
(6.8.11a) - (6.8.11b) existing on the slab [T_,T+] x M 3 , then (6.8.9a) - (6.8.9b) are satisfied in the same slab; i.e., 



the well-known identity div o curl = implies that Vtdiy B = Vjdiy D - 0. 



Remark 6.8.2. Using definition (4.2.2) and the relation (6.8.8a), and performing some tedious calculations, it 
follows that 



(6.8.12) 



f 2 

l (MBI) 



(1+1 


B 


2 ) 


2 


1 + 


B\ 2 + 


D 


12 


+ 


B x D\ 2 



Consequently, as shown by Proposition 8.4.2 and Proposition 13.0.4, the regime of hyperbolicity for the MBI 
system has a particularly nice interpretation in terms of the state-space variables (B, D) ■ the MBI equations are 
well-defined and hyperbolic for all finite values of (B,D). 



7. Commutation Lemmas 
In this section, we prove some basic commutation lemmas that will be used in the following sections. 
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Lemma 7.0.1. Let X and Y be vectorfields, let U be any tensorfield, and let V denote the Levi-Civita connection 
corresponding to the Minkowski spacetime metric. Then 



(7.0.1) 



Proof. Simply use (3.2.9) and the fact that Minkowski space is fiat. 



□ 



Lemma 7.0.2. Let O and Z be the subsets of Minkowski conformal Killing fields defined in (5.0.8b) - (5.0.8c 
Then for any Z e Z, we have that 



\VZ\<1, 



V(2)^l = 0. 



(7.0.2) \Z\ < s, 

Furthermore, for any O e O, we have that 
(7.0.3) \0\<r, |vO| = const, |V (2) O| = 0. 

Proof. The simple computations are easily performed in the inertial coordinate system. 



□ 



Lemma 7.0.3. Let V denote the Levi-Civita connection corresponding to the Minkowski metric, and let I denote 
a multi-index for the set Z of Minkowski conformal Killing fields. Let C z be the iterated modified Lie derivative 



from Definitions 5. 0.2 and 5. 0.3, Then 



(7.0.4) 



[v,/£] = o, 



[v,c z ] = o. 



Ln an arbitrary coordinate system, equations (7.0.4) are equivalent to the following relations, which hold for all 
type (^) tensorfields U ■ 



(7.0.5) 



Proof. The relation (7.0.4) can be shown via induction in |/| using ( 3.2.10[ ) and the fact that V (2)Z = (i.e. 

p3|). □ 



Lemma 7.0.4. Let L denote a multi-index for the set Z of Minkowski confor mal Killin g fields, and let T be a 



5.0.2 



and 



5.0.3, Then 



two-form. Let C z be the iterated modified Lie derivative from Definitions 
(7.0.6) ^{[(g-^ig- 1 )^ - (<f ^"V 1 )""] V^a} = [(<f ^""(if ^ - (<f ^"V 1 )""] V „C z T kX . 



Proof. Let Z e Z. By the Leibniz rule, (5.0.9b), and Lemma |7.0.3[ we have that 



(7.0.7) ^{[(g-^ig- 1 )^ - Gr^V 1 )™] V.a} = -2c z [(g- 1 )» K (g- 1 y x - Gr^V")™] V M ^A 

+ [(g-y^g- 1 )^ - (g-YHg-'r^aCz^x. 

It thus follows from Definition 15.0.21 that 

(7.0.8) ^{[(g-'r^g- 1 )^ - (g- l Y x (g- l r K ]v^K>] = [(g- 1 )^' 1 )^ - Of ^"V 1 )*"]? ^c z t kX . 



This implies (7.0.6) in the case \I\ = 1. The general case now follows inductively. 



Lemma 7.0.5. If O e 0, then the vectorfields, L,L,0 mutually commute: 
(7.0.9) [L,L] = 0, [L,O] = 0, 



[L,O] = 0. 



□ 
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Furthermore, with q = r - t, s = r + i denoting the null coordinates, and S = x K <9 K denoting the scaling 
vectorfield, the following differential operator commutation relations hold when applied to arbitrary tensorfields: 



(7.0.10a) 
(7.0.10b) 
(7.0.10c) 
(7.0.10d) 



[Vl,Vl] = 0, 
[Vl,Vo] = 0, 
[Vi,£ O ] = 0, 
[?Vl,£s] = 0, 



[?Vl, sVl] 


= 0, 


[Vl,Vo] 




[Vl,£o] 




[sVl,£ 5 ] 


- 0. 



Proof. (7.0.9) follows from simple computations. ( 7.0.10a[ ) then follows from (7.0.1), (7.0.9), and the fact that 
Vls = V Lg = 0. fl7.0.10b[) also fo llows from ( friol) ) and ( |7.0.9( . ( |7.0.10cp follows from fl3.2.10[ ), ( |7.oT3| >, and ( |7.0.10b[ ). 
( 7.0.10d[ ) follows from flV.O.lOaj ), identity 2S = sL - qL, and the fact that by fl3.2.10[ ), £ s t7 = V S U + (m-n)U for 
a type (,™ ) tensorfield £/. □ 



Lemma 7.0.6. Let O e O, and let e K x^u, 
Then 



and £uu be as defined in (3.3.1), (3.4.6), and (3.4.8) respectively. 



(7.0.11a) 
(7.0.11b) 
(7.0.11c) 
(7.0.11d) 
(7.0.11e) 
(7.0.11f) 



Vl^ = Vl ^ u = 0, 
Vl/p = VUtw = o, 
fx = 0, 

£o ^ = o, 

C>oiiiv - 0. 



Proo/. T he relation Q7.0.11a[) follow s from definiti on fl3.2.7[ ), ( |3.4.6D , a nd (|6.1.2aD - fl6.1.2b| ). Using also ( |3.3.3[ ), 



(7.0.11b) follows similarly. (7.0.11c) follows from (3.3.3), the formula (3.4.12), and the fact that the null second 



fundamental forms of the spheres S r j, are tangent to S r j (he., Lemma 3.4.1). (7.0. lid) follows from (5.0.15) and 
the fac t that (0) 7r^ = (i.e., O is a Killing field). ( |7.0.11eD and ( |7.0.11f| ) follow from definitions pXIj ) and ([3X6]), 
(177091), and d7.0.11d|). □ 



Corollary 7.0.7. Let T be a two-form and let a, a, p, a be its null components. Let O e O. Then CoaiJ 7 ] = 
a[£o^ r ], £o«[-^"] = a [^-o^], ^-opi^] = p[CoF], and Coo~[^F] = cr[Co^]- An analogous result holds the operators 
Vl and Vl; i.e., Cq,Vl, and Vl commute with the null decomposition of T . 



(7.0.11f) 



Proof. Corollary 7.0.7 follows from Definition 6.3.1, (6.1.2a) - (6.1.2b), (7.0.9), (7.0.11a) - (7.0.11b), and (7.0. lie) 



□ 



Lemma 7.0.8. Let X and Y be vectorfields, and let J- be a two-form. Let ix3~ be the interior product defined by 
the requirement that ix3~(Y) = T(Y,X) holds for all pairs of vectors X,Y. Then 



(7.0.12a) 
(7.0.12b) 



ixZyF - ZyixF = i[x,Y]F, 



Proof. The relation (7.0.12a) follows from the Leibniz rule and the fact that -CyX = [X, Y~\. The relation (7.0.12b) 
follows from the Leibniz rule. □ 



We leave the computations required to prove the next lemma up to the reader. 
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Lemma 7.0.9. [CK90, pg. 139] The 15 generators {T^), ^(^), S, -f^) }o</i<;/<3 of the Minkowski conformal Killing 
fields, which are defined in Section^ satisfy the following commutation relations: 



(7.0.13a) 


T( ,,\ i T/jA 


= o, 


(m,z; = 0,1,2,3), 

\ / ' 3 ' 7/3 


(7.0.13b) 


\ 1 ( ^\ , 1 It ...A 




(k, u, ^ = 0, 1, 2, 3), 


(7.0.13c) 


L (W' J 




(ft, //, f = 0, 1, 2, 3), 


(7.0.13d) 


[^( M )' T (^)] 


= 2^5 + 2£2( AtI/ ), 


( M) i/ = 0,1,2,3), 


(7.0.13e) 


[^(«A))^(^)] 




(k, \,fj,,v = 0, 1,2,3) 


(7.0.13f) 






( M ,i/ = 0,1,2,3), 


(7.0.13g) 




= 9khK(v) - 9tiuK^, 


(M,v = 0,1,2,3), 


(7.0.13h) 




= K {^ 


(^ = 0,1,2,3), 


(7.0.131) 


[if (At) ,if (!/) ] 


= 0, 


(0,^ = 0, 1,2,3). 



□ 

The following simple corollary follows directly from Lemma |7.0.9 

Corollary 7.0.10. Let T and Z denote the Lie algebras of vectorfields generated by the sets 3~ and Z respectively. 
Then for fj, = 0, 1, 2, 3, we /iaue 



(7.0.14a) 
(7.0.14b) 



[T (M) ,Z]cT, 
[5,Z]cT. 



□ 



Lemma 7.0.11. The null second fundamental forms \l ^Lf and V ^L v (see Definition 3.4-4) can be expressed as 



V U L? = -r 



(7.0.15a) 
(7.0.15b) 

Furthermore, the intrinsic covariant derivatives of VL_ and VL vanish: 



■i j. v 



(7.0.16a) 
(7.0.16b) 



t{M) L =V(M-i) VL = 0, 



^(M) ^ =V (M-l) VL = 0, 



(M>2), 
(M > 2). 



Proof. (7.0.15a) - (7.0.15b) follow from simple computations. (7.0.16a) follows then follows from the Leibniz rule, 
( |3.4.19[ ), and the fact that f r = 0. □ 

Lemma 7.0.12. Let U be a type (^) tensorfield tangent to the spheres S r j- Then for all integers k,l>0, V^V^C/ 
is a tensorfield tangent to the spheres S r t, and 



(7.0.17) 



fx (V' L V|C/ W . 



-jfer-VrV*" 1 



L v L 



'Mm 



Proof. The fact that V^V^l/ is tangent to the spheres follows from contracting any of its indices with either L or 
L, and using (6.1.2a) - (6.1.2b) to commute the contractions through the derivatives; the result is 0. The relation 
(7.0.17) follows from the Leibniz rule, (3.4.12), and Lemma 7.0.11 □ 
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Corollary 7.0.13. Let k,l,m>0 be integers, and let U be a tensorfield tangent to the spheres SV,t- Then 



(7.0.18) 



k'+l'<k+l 

l'<Lk'<k,m'<m 



_ { (fe-fc') + (i-i')}| fc' i' 



Proof. Inequality ( 7.0.18 ) follows from repeated use of the Leibniz rule, Lemma 7.0.12 , and the fact that ~f r = 0. □ 



The following simple lemma gives pointwise bounds for the tensorfields L,L, <//, and their full spacetime covariant 
derivatives. 



Lemma 7.0.14. Let M > be any integer. Let L,L be the null vectorfields defined in (3.4.3a) - (3.4.3b), and let 



$ be the first fundamental form of the S r j defined in (3.4.6). Then 



(7.0.19a) 
(7.0.19b) 

and 



V 



(Af) 
7 (M) 



L\<r 
L\<r 



-M 



-M 



(7.0.20) 



V 



(M) 



< f 



-M 



Proof. Since in the inertial coordinate system, Lf 1 = (-1, w 1 , uj 2 , w 3 ), L M = (l,a/,u; 2 ,u; a ), and u J d = f x J /r, it is 



easy to check directly that (7.0.19a) - (7.0.19b) hold. Inequality (7.0.20) then follows from (3.2.7), definition 
(3.4. 6P, (17.0. 19ap, fl7.0.19b|), and the Leibniz rule. □ 



Lemma 7.0.15. Let OzO. Then 



(7.0.21) \i{M)0\<r 1 - M . 
Proof, follows from repeated use of ( |3.4.12[ ), ( |7.0.3 ), and Lemma 7.0.14 



□ 



an enumer- 



Lemma 7.0.16. Let U^...^ be a type (°) tensorfield tangent to the spheres S r) t- Let O l ,l - 1,2,3 be 
ation of the three rotational vectorfields belonging to the set O. For any rotational multi-index L = (ii,--,t m ) of 

length m, where Li e {1,2,3}, we define |0 7 C/| = \0^--- 0^™U Kl ... Km \. Then the following pointwise estimate holds: 



(7.0.22) 



Y |0 7 ^|^r m |C/|. 



Proof. It is easy to check that at any point in p e SV,t; a pair of the rotations, say Ll and O t2 , has the following 

two properties: i) \0 Ll \ 2 ,\0 L2 \ 2 > r 2 /3; ii) | j, (O tl ,O t2 )| 2 < \\0 Ll \ 2 \0, 2 \ 2 (where $ (X,Y) d % A X K Y X ). This latter 
property implies that the (smallest) angle between Ll and O i2 , viewed as vectors in the two-dimensional plane 
T p S r> t, is at least 60 degrees. Therefore, any covector £ e T*S r> t, has a corresponding ^ -dual vector X e T p S r j 

that makes an angle < 60 degrees with one of tO L1 , tO L2 , and it follows that either | (0 Ll ,X)\ 2 > ^\X\ 2 = y£| 2 
or I 4 (0 L2 ,X)\ 2 > r -^\X\ 2 = ^|£| 2 . Thus, E|/|=i |C 7 £| 2 £ r 2 \£\ 2 . 

On the other hand, the reverse inequality E|/|=i |C 7 £| 2 ~ r2 |£| 2 trivially follows from (7.0.3). We have thus 



shown (7.0.22) in the case m= 1. The general case ( 7.0.22| ) follows from applying similar reasoning to the elements 



{e* <8> ■■■ ® e„* } of a <h -orthonormal basis for T*S r t 



1 Tp S r j 



m copies 

Lemma 7.0.17. Let N > be an integer. Then for any tensorfield U, we have that 
(7-0.23) \U\ Co -N - |^| Vo ;JV- 



□ 
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(7.0.24a) 
(7.0.24b) 

(7.0.24c) 



N 

\U\ fo ;N«Z rn \f(n) U\, 
n=0 

P\fo;N ~ \U\co;N, 
N 

\U\ Co ;N*Z rn \f(n) U\. 
n=0 



The norms in the above inequalities are defined in Definition \9. 0.1 
Proof. We use the notation denned in Lemma 



relation (3.2.10), together with the inequalities 



7.0.3). 



7.0.16 Inequality (7.0.23) can be proved inductively using the 



To prove (7.0.24a), we first note that Lemma 7.0.16 implies that for each integer N > 0, the inequalities in 



(7.0.24a) are equivalent the following inequalities: 



(7.0.25) 



P\ fo ;N = f E I foU\» E P 1 t Oil) U\. 
\I\<N \I\<N 



To prove (7.0.25), we argue by induction, the base case N - being trivial. For the induction, we assume that the 



inequality is true in the case N. Let / = (ii, fcfc) be a rotational multi-index with \I\ = k < N + 1. Then by the 
Leibniz rule and (7.0.21), we have that 



(7.0.26) 



K rZ 

\to n - Uo Lk u-o%-o% f ai - fa k u\<Y, E (III o H \)\ f (k - p) u\ 



foV 



o'f w u 



p=la!+---+a fc =p i=l 

0<ai<i-l 

fc-1 



p'=0 



where in the last step, we have used (7.0.24a) under the induction hypothesis. Summing over all \I\ = k < N + 1, 
we conclude that 



l^o;iv + i- E P 1 tm u \\ * l^l^oiiv ^ E P 1 u l 

\I\<N+1 



'\<N 



where in the last step, we have used (7.0.25) under the induction hypothesis. From (7.0.27), the induction step 



for (7.0.25) easily follows. 



To prove (7.0.24b), we also argue by induction using the identity (3.2.10) (which is valid if V is replaced with 
1 f), the base case being trivial. For the induction, we assume that inequalities hold in the case N. Let I be any 
rotational multi-index with 1 1\ = k < N + 1. Repeatedly applying the Leibniz rule and identity (3.2.10), and using 
( |7.0.21 ) plus (7.0.24a), we deduce the following inequalities: 



(7.0.27) 



C'oU- f u\<Y, E ( EI I tM Cy)l f(k- P ) u\ 



p=lai+---+a k =p i=l 
0<Oi 

fc-1 



p'=0 



Summing over all |J| = k < N + 1, we conclude that 



(7.0.28) 



\U\c a ;N+l ~ \U\f a -N+l % \U\f a -N £ \U\c a ;N 
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where in the last step, we have used the induction hypothesis. The induction step easily follows from (7.0.28), 
which completes the proof of (7.0.24b). 



The estimate (7.0.24c) then follows trivially from (7.0.24a) and (7.0.24b). 



□ 



Corollary 7.0.18. Ifk,l,M > are integers, and U is any tensorfield tangent to the spheres S r j, then we have 
the following estimate: 



(7.0.29) 



M 



M 
m=0 



vM f (m ) u\. 



The norm on the left-hand side of (7.0.29) is defined in Definition \9.0.1\ 



Proof. Corollary 7.0.18 follows from Corollary 7.0.13 and Lemma 7.0.17 



□ 



8. The Energy-Momentum Tensor and the Canonical Stress 

In this section, we discuss the building blocks of our energies. We will begin by defining the MBI system's 
energy-momentum tensor, which we denote by Q^ M bi) > and recalling its key properties. We remark that Q^mbi) 
is the usual tensor associated with energy estimates. However, in order to derive energy estimates for the deriva- 
tives of a solution, we will need a different tensor, namely the canonical stress Q v . The bulk of this section is 
therefore devoted to an analysis of the properties of Q%. 



8.1. The energy-momentum tensor Q^ u . The energy- momentum tensor corresponding to an electromagnetic 
Lagrangian *Jz? is defined as follows: 



3.1.1) 



dg, 



(ji,v = 0,1, 2, 3). 



fa' 



It follows trivially from the definition that Q^ u is symmetric: 



(8.1.2) Q"" = Q"i*. 

We recall the fundamental divergence-free property: for energy-momentum tensors Q^ u constructed out of a 



solution of the equations of motion (4.1.6a) - (4.1.6b) corresponding to *Jzf, we have that 



3.1.3) 



VuQ^ = o, 



= 0,1,2,3). 



In the particular case of the MBI model, (4.2.3) and Lemma 4.1.1 imply that 



3.1.4) 



^ (MBI) 



The next lemma is not needed for any of the main results presented in this article, but it is of interest in itself. It 
shows that Q^f BI ^ satisfies the dominant energy condition. 



Lemma 8.1.1. Let J- be any two-form for which the quantity £(mbi)j which is defined in (4.2.2) is real; i.e., any 



T for which 1 + ^[F] - ^ 2) [-? r ] > 0. Then the MBI energy-momentum tensor of T satisfies the dominant energy 
condition; that is, for any pair of future-directed causal vectors X,Y we have that Q(mbi)(X,Y) > 0. 

Proof. Let p e M, and let X, Y be any future-directed causal vectors belonging to T p M. Then in the plane spanned 
by X, Y, there exist a pair of null vectors L' , Jj such that g(L_' , L') = -2, and such that X = aV + hfJ , Y = cL' + dLf, 
with a, b,c,d> 0. Let us complement L' , V with a pair of orthonormal vectors e[, e' 2 belonging to the g-orthogonal 
complement of span{L',L'}. The set {L' , L' , e' 1; e' 2 } is therefore a null frame at p. Let a',a',p',cr' denote the 
null decomposition of T with respect to this frame. By bilinearity, it suffices to check that Q(mbi){Ll \ Ll) ^ 
0, Q(mbi){L' i L') > 0, and Q(mbi)(L.' , L') > 0. We leave the following two simple calculations to the reader: 
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(8-1-5) Q(MBI)(liM) = ^(MB/)I^'| 2 ) Q(MBI)(L\L') = £(MBI)\ a '\ 2 - 

For the remaining term, we first calculate that 
(8-1-6) ^mbi)Q{mbi)(L',L') = g^ A a ,A + 2p' 2 + 2^ 2) + 2l {MBI) (l {MBI) - 1). 

We now express 2l {MBI) (£ {MBI) - 1) = - tf 2) ), where f(v) = 2(1 + v) - 2(1 + t;) 1 / 2 , and /(0) = 0, /'(0) = 1, 

f"(0) = 1/2. Setting v = - ^ 2 2 ^ , using the convexity of /, and using the identity = -a^a ,j4 - p' 2 + er' 2 , it thus 
follows that 

(8.1.7) 1(MBI)Q(MBI)(I/,L') > g! A a ,A + 2p' 2 + 2^ 2 2) + - ^ {2) ) > p' 2 + a' 2 + ^ 2 2) . 

This concludes our proof of Lemma |8.1.1| 

□ 



8.2. Equations of variation. The definition (9.0.7) of our energy involves integrals of weighted squares of the 
components of C^J- over the spacelike hypersurfaces £j. In order to prove our global existence theorem, we need to 
understand the evolution the energy, which in turn requires that we investigate the evolution of the C^J-. These 



quantities are solutions to the equations of variation, which are the linearization of the MBI system (4.3.2a) - 
(4.3.2b) around the background T. More specifically, the equations of variation in the unknowns T^ u are defined 
to be 



3.2.1a) 
3.2.1b) 



Va>7>z/ + V m ^va + V^A/i - 3 a 



(A, //,*/ = 0,1,2,3), 
(^ = 0,1,2,3), 



where the tensorfield H^ unX is defined in (4.3.1a), and Zx^v and J 1 ' are inhomogeneous terms that need to be 
specified. In our applications below, the C^T will play the role of T . In order to understand the evolution of the 
C^J-, we of course need to understand the nature of the inhomogeneous terms appearing in equations of variation 
that they satisfy. The next proposition provides detailed expressions for these inhomogeneities. 



to the equations of variation (8.2.1a) - (8.2.1b) with inhomogeneous terms given by 



Proposition 8.2.1. If J- an is a solution to the MBI system (4.3.2a) - (4.3.2b), then Tn U = C^Fuv is a solution 



3.2.2a) 



3.2.2b) r {1) = Hr x v»(£z^x) - ^(tir^^x), 



(A, ^ = 0,1, 2, 3), 
(u = Q, 1,2,3). 



In the above formula, the tensorfield H 



[ivkX 



which depends quadratically on T, is defined in (4.3.1b), while the 



iterated modified Lie derivatives £J Z are defined in Section^^ 



Proof. Recall equation (4.3.2a), which states that T is a solution to V\ K J~ iiV ~\ = 0, where [•••] denotes antisym- 
metrization. Using property (7.0.4), it therefore follows that 



3.2.3) 



= CzV[ K T^ = V^CzT^y 



This proves (8.2.2a) 



To prove (8.2.2b), we first recall equation (4.3.2b), which states that J- is a solution to 



(8.2.4) h^ kX v^ kX d = e£ |i[Gf ^Or 1 )"* - GrYV 1 )™] + Ff^lv^A = o. 

Applying C J Z to ( |8.2.4[ ) and using (17.0.61), we conclude that 
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(8.2.5) 
It therefore follows that 



3.2.6) 



which proves ( 8.2.2b[ ) 



□ 



Although the energy momentum tensor (8.1.1) is useful for estimating Tav, 

the 



1.3.2 



8.3. The canonical stress 

it is not quite the right object for estimating its derivatives C> Z T ■ As explained in detail in Section 
reason is that the C Z J~ are solutions not to the MBI system itself, but rather to the to the equations of variation 
(8.2.1a) - (8.2.1b), whose linearized Lagrangian, which is defined in (8.3.1) below, depends on the background 



T . Nonetheless, we will be able to construct the canonical stress, which encodes the approximate conservation 
laws satisfied by solutions to the equations of variation, and which we denote by Q%. We remark that a general 
framework concerning the properties of the canonical stress was developed by Christodoulou in |ChrOO| : here, we 
only recall its basic features. As we will see, Q„ has the following two key properties: 

• Its divergence is lower-order (in terms of the number of derivatives) . 

• It has positivity properties related to, but in general distinct from those of the energy-momentum tensor 

The first property is explained in detail at the end of this section, while the second is discussed in limited fashion 
in Section [8T4l 

In order to understand that origin of the canonical stress, let us first define the linearized Lagrangian 3f, which 
is, despite its name, a quadratic form in the variations T with coefficients depending on the background T. 

Definition 8.3.1. The linearized Lagrangian ^ffj 7 ;^"] corresponding to the Lagrangian =£?[•] and the background 
T is defined as follows: 



3.3.1) 



where the ^"-dependent tensorfield hS 71 is defined in (4.1.9) 



The significance of ££ \f\ J-] is that its corresponding equations of motion are the equations of variation. More 
specifically, if we consider T to be the unknowns, then the "principal part of the Euler-Lagrange equations (assuming 
the stationarity of the linearized action ^.^[J 7 ] d = f f«r sM S£ \T\ T\ dfi g , under closed variations of F) corresponding 



to J£\T;T\ are the linearization (around J-) of the Euler-Lagrange equations (4.1.8) corresponding to ££ 



Definition 8.3.2. Given a linearized Lagrangian Jz? = S£\T\T\ we define the corresponding canonical stress 
as follows: 



3.3.2) 



where the tensorfield h is defined in (4.1.9) 



In Section 4.3, we modified the tensorfield h corresponding to the MBI system, obtaining a new tensorfield H. 



Therefore, for the purposes in this article, it is convenient to construct the MBI canonical stress using H : 
(8.3.3) 



01 = f h^ kX t kX t k - -^h^ x t^ kX , 
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where H is defined in (4.3.1a). For the remainder of the article, the term "canonical stress" is used to 
refer to the tensor defined in (8.3.3). Note that Q^ v depends quadratically on J 7 , and it also depends on the 
background T. 

In general, Q^ u is not symmetric, nor is Q„ traceless. However, in the case of the MBI system, it can be checked 
that Q^y is in fact traceless. More specifically, using (4.3.1a), we compute that 



5.3.4) 



linear Maxwell-Maxwell terms 



'(2) t (MB/), 



1 



from which it also follows that 



5.3.5) 



Q\ = o. 



We also note that the first two terms on the right-hand side of (8.3.4) are the components of the energy-momentum 



tensor for the linear Maxwell-Maxwell equations (in the unknown J-). 
The failure of symmetry of is captured by the following expression: 



5.3.6) Q M „ - Q 



1 



r 2 

1 r 2 

2"(2) l (MBJ) 



■> kA| •> •> vC, ~ •> v •> Kj- 



(2) C (MBI) 



We conclude this section by proving a lemma that illustrates the first key property of Q^ u , namely that its 
divergence is lower order. 



Lemma 8.3.1. Let T be a solution to the equations of variation (8.2.1a) - (8.2.1b) corresponding to the background 



J 7 , and let Q? v be the tensor defined in (8.3.4). Let Zxaui 3 U be the inhomogeneous terms on the right-hand sides 
of (18.2. la!) - (18.2. lbh. Then 



8.3.1 



follows from using the properties (4.1.10a) - ( |4.1.10c ), which are satisfied by , and 

□ 



(8.3.7) 

Proof. Lemma 
simple computations 

8.4. Positivity properties of Q v . The canonical stress has positivity properties that are analogous to, but 
distinct from those of the energy-momentum tensor. For a complete discussion of these properties, which are 
related to the geometry of the equation^} see |Chr00| . Here, we only discuss the positivity properties that 
are relevant to our small-data global existence proof and our sketch of a large-data local existence proof. For 
purposes of constructing an energy suitable for the global existence proof (see Section [9]), the main quantity of 

interest will be Q(^°\ K) d = f Q K X ^^ K , where is the g-dual of the timelike translation Killing vectorfield 



T( ) defined in ( |5.0.8a ), and K = |{(1 + s 2 )L + (1 + q 2 )L} is the vectorfield defined in (6.6.1). As we will see 



in the next lemma, in the small-field regime, the resulting expression is positive definite in the null components 
of the variation T . However, different components carry different weights. Ultimately, the different weights will 
translate into the fact that the various null components of a solution J- and its derivatives have distinct rates 



In general, the geometry of the electromagnetic equations (i.e., the characteristic cones) is distinct from the geometry of spacetime; 
however, in the case of Maxwell-Maxwell theory, the two geometries coincide. 
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of dec ay; see the global Sobolev inequality (Proposition |1 1.0.9 ) . For large-data local existence (see Proposition 
13.0.4), we will use a vectorfield Xi oca i (see Proposition 8.4.2) in place of K. The reason for this modification is 



that the positivity of Q(^°\K) = f R X may break down in the large-data regime. On the other hand, 

the vectorfield Xi oca i, which is constructed with the aid of the reciprocal Born-Infeld metric (b' 1 )^", maintains its 
positivity in all regimes in which the MBI equations are well-defined. Although Xi oca i does not provide good q, s 
weights for the null components of !F, a bound of the form Q(^°\ X[ oca [) > C|^"| 2 for some positive constant C is 
sufficient to prove local existence. 

We begin with a lemma that addresses the positivity properties of Q(^°\K). 

Lemma 8.4.1. Let T \T be arbitrary two-forms, and let a,a,p,&, be the null components of T. Let Q be the 
canonical stress (8.3.4) associated to J-,J-, and let T^o), K = ^{(1 + «s 2 )L + (1 + q 2 )L} be the conformal Killing 
fields defined in (5.0.5a) and (6.6.1) respectively. Let £^ ' be the g-dual ofTr^y There exists an e > such that if 
1.7-1 < e, then 



3.4.1) 



Q(e (0) , K) = Q, U T^K V = + q 2 )\a\ 2 + (1 + s 2 )|d| 2 + (2 + q 2 + s 2 )(p 2 + a 2 )} 



+ l"| 2 o((l + s) 2 \F?cu) + (1 + s) 2 (|«| 2 + p 2 + ^ 2 )0(|^| 2 ) 
+ (1 + \q\) 2 (\a\ 2 + \a\ 2 + p 2 + a 2 )0(|^| 2 ). 
In the above expression, 0{U) denotes a quantity which is < CU in magnitude for some positive constant C. 



Remark 8.4.1. As we shall see, in the small-solution regime, the dominant term on the right-hand side of (8.4.1 ) 



is l{(l+q 2 )\a\ 2 + (l + s 2 )\a\ 2 + (2 + q 2 + s 2 )(p 2 + & 2 )Y This motivates definition ( |9.0.3| ) below. Furthermore, at first 

sight, one might worry that the (1 + s 2 ) factor could cause the "error terms" O(-) to become large. This worry 
is alleviated by Corollary 1 1 1 . . 1 



Proof. Using (8.3.4), we decompose Q into its linear "Maxwell" part Q( MaxweU ) ; anc j the remaining "error terms:" 

Qtiu = Q ( ^ Iaxwell '> + -£( 2 fB/) | - T^T^T^Tkx + -g^ui^^Kx) 2 ) 

>> 2 } 



3.4.2) 



1 



4^ 
1 



2 



3.4.3) 



^(Maxwell) def j- £ j- _ 1 f ^ 



As in the proof of Lemma 8.1.1 leave it as a simple exercise for the reader to check that 



3.4.4) 



& Maxwell \L,L) = \a 



• |2 



Q 



(Maxwell) 



(L,L) = |df 



Q 



(M axwell) 



(L,L) = (p 2 + a'). 



Since Q^ Maxwell \T (0 y K) = Q( Maxwell \±(L + L), (1 + s 2 )L + (1 + q 2 )L_), it easily follows from gZg that 



3.4.5) 



Q( Maxwell \T (0) ,K) = !{(! + , 2 )|d| 2 + (1 + S 2 )|d| 2 + (2 + g 2 + s 2 )(p 2 + a 2 )}, 



which gives the principal term (i.e., the term in braces) on the right-hand side of (8.4.1). Note that in the above 



formulas, we have abused notation by identifying covectors with their g-dual vectors. 

For the error terms, we first note that ^mbi) * s an or der 1 factor when T is sufficiently small, so that we 



may ignore it. With this fact in mind, we then evenly divide the 8 error terms on the right-hand side of (8.4.2) 
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6.4.1 



into two classes: those that contain the g^ factor, and those that do not. For the first class, we note that 
\g^T^ 0) K u \ = ||2 + q 2 + s 2 \ < (1 + s 2 ). Therefore, using Lemma 

(8.4.6) 



it follows that 

2^1- 1 2 , -2 , -2> 



\g, v T^K\T KX T RX f\ < (1 + S 2 ){|^|LI«! 2 + l^| 2 (l«! 2 + P 2 + v 2 )}, 



and similarly for the other 3 terms of this type. 

For the first term belonging to the second class, we use the null decomposition (6.6.2a) - (6.6.2c) of K to 
conclude that 

(8.4.7) IT^ZV/^I < (1 + q 2 )\T\\T\ + (1 + s 2 )\T\\F\cu- 

Therefore, by a second application of Lemma 6.4.1 and the inequality \ab\ < a 2 + b 2 , it follows that 
(8.4.8) 

\T^K v T£f u ^ x T K> \ < (1 + g 2 )\T\ 2 \T\ 2 + (1 + s 2 )\T\ 2 cu \T\ 2 + (1 + s 2 )\T\ 2 \T\ 2 cu + (1 + s 2 )\T\ 2 \T\ 2 TT 



< (1 + g^l^ffldf 4 + |d| 2 + p 2 + o*) + (1 + s 2 )|^|a| 2 + (1 + s")|.7f (|df + p* + ^) 



• 2\ 



The remaining 3 terms of this type can be estimated similarly. Inequality (8.4.1) thus follows. 



□ 

The next proposition provides the fundamental estimate that is needed to deduce the "large-data" local existence 



result of Proposition 13.0.4 



Proposition 8.4.2. Let Jif d = {T \ £(mbi)[F] d - (l + ^(i) [•?""] - L f%\ [-T""]) 1 ^ 2 > 0} denote the interior of the subset 
of state-space for which the Maxwell-Born-lnfeld Lagrangian is well-defined, and let & be a compact subset of J^. 



Let T € A, let T be any two-form, and let Q„ be the canonical stress (8.3.4) for the MBI equations of variation 



corresponding to the background J- and the variation T. Let Xi oca i be the vectorfield defined below in (8.4.15) 
Then there exists a constant C > 0, depending only on R, such that the following inequality holds: 



m {Q \x local ) 



def 



(8.4.9) 

In the above expression, the covector £W is the g-dual of the the time translation Killing field T^, which is defined 
in (5.0.5a). 



Remark 8.4.2. Let us view Q^u^X^^ as a function of the spacetime point p and of the covector £ e T*M, 



i.e., as a function on T*M. Then since inequality (8.4.9) is strict, there is an open neighborhood j\f c T*M of 



(p,£(°)) and a uniform constant Cj\f > such that Q^u£,^Xi ocal - Cjv|^| 2 holds for all (p, £) e j\f. Using the methods 
described in |Chr00| Chapter 5], such an inequality implies the hyperbolicity of the MBI equations of variation 
and finite speed of propagation for their solutions. 



Remark 8.4.3. Since (8.4.9) is a pointwise inequality, an analogous inequality can be shown to hold on a neigh- 



borhood of any point p belonging to any sufficiently smooth Lorentzian manifold (simply use a local coordinate 
system in which g^ u \ p = diag(-l, 1, 1, 1), define to be the g-dual of , and define X\ ocal as in (8.4.15) below). 



Thus, inequality (8.4.9) implies the hyperbolicity of the MBI equations of variation for any sufficiently smooth 
Lorentzian manifold. 



Proof. Let (E, B) and (E, B) be the electromagnetic decompositions of T and T described in Section 
simple calculations imply that the following identities hold in the inertial coordinate system: 



6.8 



Then 



.4.10) 
.4.11) 
.4.12) 
.4.13) 



T K 'ZT — T? T? a 
■Tq J-Qk, - & a h, , 



^Ok 



-B n E a , 



T kX T kX = 2B a B a - 2E a E a , 
*T kX J- k \ - 2B a E a + 2E a B a . 
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We now introduce the reciprocal Maxwell-Born-Infeld metric (b 1 ) At!/ , which has the following components 
relative to an arbitrary coordinate system: 



(8.4.14) (b- l Y u = f (g- l Y v - (1 + ^(i)^) -1 ^^. 

In a future article, we will discuss the significance of (b" 1 )^ in detail. For now, we simply make the following two 
claims. 

(1) If {g-^Ux < 0, then (b^U < 0. 

(2) If the covector £ M belongs to the characteristic subset 20 of T*M, then the vector d = f (b~ 1 ) flK £ :K belongs 
to the characteristic subset of T p M. 

The claim (1), which is easy to check in the inertial coordinate system, shows that any covector £ that is timelike 
relative to (g 1 ) pl/ is also timelike relative to (6 1 )^^. This is related to the fact that the energy-momentum tensor 



satisfies the dominant energy condition. As is illustrated below in (8.4.15) the claim (2) allows us to construct a 



suitable (for deducing local existence) multiplier vectorfield from a g-timelike covector £. For, as is discussed 
in |Chr00| . in order for / R3 Q^^X" d 3 x to be bounded from below by a multiple of J R3 {Fl 2 d?x, it is sufficient^ 
to choose £ to be a covector lying strictly inside of the characteristic subset of T*M, and X to be a vector lying 



choose £ to be any g-timelike covector, and define X^ d = f (b~ 1 ) flK £ :K . 

Rather than relying on the above abstract framework, we will instead directly show the positive definiteness of 
the quadratic form (in J 7 ) Q il v Z, il X v for a particular choice of £ and X; this pointwise positivity is stronger than 
the integrated positivity that follows from the general framework. We choose £ to be equal to the covector 
that is g-dual to the timelike Killing field T/q\. Furthermore, motivated by the above discussion, we define the 
multiplier vectorfield Xi oca i by 



strictly inside of the characteristic subset of T p M satisfying 22 £(X) d = f £, K X K < 0. By claims (1) and (2), we can 



3.4.15) 



X 11 

local 



= 2ij 



where the last equality is valid in the inertial coordinate system. We remark that the 2^ 2 MBJ ^(1 + ^(i) [.?""]) factor 
on the left-hand side of (8.4.15) is a normalization factor that was chosen out of computational convenience. Using 



the identities (8.4.10) - (8.4.13), it can be checked that 



3.4.16) 



Q^Kcai = 2(1 + \B\ 2 )l\ MBI) Q m 



9P 2 

(MB!) 



Op- 



whenever T e R. Using (8.3.4) and (8.4.10) - (8.4.13), and defining (for notational convenience) 



Our goal is to show that the right-hand side of (8.4.16) is a uniformly positive definite quadratic form in J- 



(8.4.17) 

for any pair of vectors U, V 6 St, we compute that 



(U,V) d ^U a V a 



3.4.18) 



2(: 



, mbi)Qoo - ^mbi){\E\ 2 + \B\ 2 } + {(E,E) 2 + (1 + \B\ 2 - \E\ 2 )(B, E) 2 + 2(E, B)(E, E)(B, E)} 
- {(B,B) 2 + (1 + |i?| 2 - |£| 2 )(£,i?) 2 - 2(E,B)(E,B)(B,B)}, 



20 The characteristic subset of TpM is defined in Section [l] 

21 As discussed in [ChrOO] , in general, we must choose £ to be an element of the inner core of the characteristic subset of T*M, and X to 
be an of the inner core of the characteristic subset of T P M. However, in the case of the MBI system, the inner core of the characteristic 
subset of TpM coincides with the interior of the characteristic subset of TpM, and similarly for the inner core of the characteristic 
subset of T P M. This is because the characteristic subset of T*M in the case of the MBI system is particularly simple: it is the cone 

22 Many of our definitions differ from those in |Chr00j by a minus sign. 



and 
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(8.4.19) f {MBI) QoaF ap Fo P = -(B, B) 2 {\E\ 2 + (E, B) 2 } + (B, B)(B, E)(l + \B\ 2 )(E, B) 

+ 2(B,B)(E,B)(1 + \B\ 2 )(E,B) + (B,B)(E,E)(1 + \B\ 2 ) 

- (E,B) 2 {(1 + \B\ 2 - \E\ 2 )\B\ 2 + (E,B) 2 } - (E,B)(B,E)(l + \B\ 2 - \E\ 2 )(l + \B\ 2 ) 
-(E,B)(E,E)(E,B)(l + \B\ 2 ). 

To simplify the subsequent analysis, we choose a frame of g-orthonormal vectors {e^,e L ,e x } c "E t such that 
E € spanjeii}, and such that B e span{e|| , e L }. We can therefore express (in a slight abuse of notation) the electric 
field as Ee\\ and the magnetic induction as -B||e|| +B 1 e 1 . Furthermore, we can decompose B and E relative to this 
frame as 



3.4.20) 
3.4.21) 



B = B\\e\\ + B ± e ± + B x e x , 
E = Eueu + E ± e L + E x e x , 



where 



(8.4.22) 
(8.4.23) 
(8.4.24) 
(8.4.25) 
and 



<£,e||>= f % 

(E,e [{ )^E h 
(E,e L )^E x , 



.4.26) 
.4.27) 



\B\ 2 = B 2 + Bl + Bl, 



\E\ 



Eu+Ef + E 2 . 



Let us first address the positivity of the quadratic form Q^vi^ X^ ocal in (B X ,E X ), since these components are 
very simple to analyze, as they are completely decoupled from the remaining components. The only term that 



contributes to these components is the ^ MBI ^\\ 
to the trivial inequality 



E\ 2 + \B\ 



! )j term on the right-hand side of (8.4.18), which leads 



3.4.28) 



C)V t(o) x u 

V jA/t ^lc 



^s^t ^local - ( '{MBI) 



{l + \B\ 2 ){E 2 x+ B 2 x ). 



We note that the uniform positivity of the right-hand side of (8.4.28) for T e .ft is manifest, since £ 2 mbi)[^~] ^ s 



uniformly positive whenever T e We therefore conclude from (8.4.28) that 



3.4.29) 



Q^X^yC^El + Bl). 



We now address the components (Bn , B ± , Eu , En ). After many tedious but simple calculations, and ignoring the 



components (B X ,E X ) analyzed above, it follows that the right-hand side of (8.4.16) can be viewed as a quadratic 
form in the components (Bu,B 1 ,Eu,Eu), whose corresponding symmetric matrix {^4r/}i<jj<4 has the following 
entries: 
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(8.4.30) 


A n = 


(t 2 (MBi) + E 2 Bl)(l + Bl- 


(8.4.31) 


Al2 = 


-^\mbi) + e2b I) b \\ b l, 


(8.4.32) 


A13 = 


(l + lSl 2 )^!,^, 


(8.4.33) 


A u = 


(1 + |£| 2 )(1 + B\- E 2 )EB 


(8.4.34) 


A 2 2 = 


2 ^)(l + 5jf), 


(8.4.35) 


A 23 = 


-(l + |fl| 2 )(l + Sf)^B ± , 


(8.4.36) 


A 24 = 


-(l + |i?| 2 )^l|S 2 , 


(8.4.37) 


A 33 = 


(1 + |B| 2 ) 2 (1 + B 2 ), 


(8.4.38) 


^34 = 


(1 + |5| 2 ) 2 S|]B ±J 


(8.4.39) 


^44 = 


(l + |5| 2 ) 2 (l + 5 2 -£ 2 ). 



By Sylvester's criterion, the positive definiteness of {^j}i<ij<4 is equivalent to the positivity of the following 
four quantities, which are the determinants of an increasing sequence of sub-blocks along the diagonal: 



3.4.40) 
3.4.41) 

3.4.42) 

3.4.43) 



An = (£( MBI) + E'Bf)(l + Bt-E'), 
det({Ay}i< ii;; <2j = {t\ MBI) + E2B l) 2 t 2 (MBi)> 
det({Ai Ws) = (1 + Bf)(l + |£?| 2 ) 2 (^ 2 MB/) + E 2 B 2 L )£f MBI) , 
det({A ii } 1 < i)i < 4 ) = (1 + |S| 2 ) 4 ^ MB7) . 



We remark that we calculated the right-hand sides of (8.4.40) and (8.4.41) by hand, while we used version 11 of 
Maple to compute the right-hand sides of (8.4.42) and (8.4.43). 



As before, the uniform positivity of the right-hand sides of (8.4.41 ) - (8.4.43 ) for T € .ft is manifest. The uniform 
positivity of the right-hand side of (8.4.40) follows from the fact that 



3.4.44) 



p2 {MBi)i^ = S 



E 2 = l + 



m 

1 + Bl l + B 



2 ' 



We have thus shown that Q^ffi S X\ ocal > C&{B 2 + B 2 + E 2 + E 2 ). Combining this bound with Q8.4.29| ), and using 



the relation (|9.0.9a|), we have that Q^f ] X^ ocal > Ca(B 2 + B 2 + B 2 + E 2 + Ef + El) 



This completes our proof of (8.4.9) 



C^\E\ 2 + \B\ 2 )>C^\ 2 . 

□ 



9. Norms, Seminorms, Energies, and Comparison Lemmas 

In this section, we introduce a collection of norms and seminorms that will be used in the remaining sections 
to estimate solutions to the MBI system. We also introduce the energy £/v = ^N[^(t)], a related positive integral 
quantity that is constructed via the canonical stress and the "multiplier" vectorfield K. In Section 



12 



we will study 

the time derivative of £n, and in order to close the estimates, we need to prove inequalities that relate the norms 
to the energy; we provide these inequalities in Proposition 



9.0.3 



We also introduce another norm 



I iff 



on the 



electromagnetic initial data (D,B), and for small-data, we prove that this norm is equivalent to £tv[.F(0)]. This 
allows us to express the global existence smallness condition of Theorem [l] in terms of (D,B), which are inherent 
to the Cauchy hypersurface Eo- In particular, the norm on (D,B) does not involve normal (i.e. time) derivatives. 
Finally, in (9.2.2), we provide a preliminary expression that will be needed in our estimates of J-(^[. F(t)]). This 
computation motivates Section 10 in which we provide algebraic estimates of the terms appearing in ( |9.2.2 ). 
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Definition 9.0.1. Let N > be an integer. If A is one of the sets O, or Z defined in (5.0.8a) - (5.0.8c), and 
U is any tensorfield, then we define the following pointwise norms of U ■ 



(9-O.la) \U\l A . N d = ef £ \£aU\ 2 , 

\I\<N 

(9.0.1b) P$ Ailf & E |V^£/| 2 . 

|/|<JV 

Furthermore, if U is tangent to the spheres S r ,u and 'f denotes the Levi-Civita connection corresponding to ^, 
then we define 



(9.0.2) 



\u\; 



;iV 



= f E I fo u\ 2 . 



\I\<N 



In the above formulas, the pointwise norm | • | is defined in (3.1.7), while the iterated derivatives V^, etc., are 
defined in Definition 15.0.31 

Definition 9.0.2. Let Q( Maxwel1 ) b e the Maxwellian canonical stress corresponding to the two-form T defined in 



.4.3), and let a,a,p,& be the null components of T. Let A be one of the sets ^,0, or Z defined in (5.0.8a) 



5.0.8c), and let K be the conformal Killing field defined in (5.0.6a). Then for each integer N > 0, we define the 



following pointwise norms of J- ■ 



(9.0.3) 
(9.0.4) 



+ T \ 2 d =° f q(^(^)^ = (i + q f\gf + (i + s 2)| d |2 + (2 + q 2 + s 2 )(p 2 + a 2 ), 



\I\<N 



where is the </-dual of the time translation vectorfield 7( ). 



Notice that the different components of T in (9.0.3) carry different weights. Additionally, we have that 
{l + \q\) 2 \T\ 2 <\T\ 2 . 



Definition 9.0.3. Let Q be the MBI canonical stress (8.3.4) corresponding to the "background" T and the 
variation and let K be the conformal Killing field defined in (5.0.6a). We define the energy current Jp[J~] 
corresponding to T to be the following vectorfield: 

(9.0.5) d =-Q^K u . 

We note that J^lT] depends on the background F through Q^ u , and that = Q(^°\K). 

Definition 9.0.4. Let N > be an integer, and let T and T be a pair of two-forms. We define the weighted 
integral norm \ T{t) %c z ;N of T as follows: 



(9.0.6) | T{t) % Cz ;N = f [f v \ Ht,*L) \l z] N d'x) 

Furthermore, we define the energy £jv[-^(*)] of T as follows: 



1/2 



(9.0.7) 



\\I\<N JR 



1/2 



where the component Jjr[T{t,x)^ is defined in (9.0.5) 
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Remark 9.0.4. The q, s weights under the integral in definition (9.0.6) are exactly the ones needed in the global 



Sobolev inequality (Proposition 11.0.9). 



In the next proposition, we collect together a large number of comparison estimates that will be used throughout 
the remainder of the article. 

Proposition 9.0.3. Let N > be an integer, and let U be any tensorfield. Then the following pointwise estimates 
hold: 



(9.0.8a) 
(9.0.8b) 

(9.0.8c) 
(9.0.8d) 
(9.0.8e) 



\U\cz;N ~ \U\y z; N, 

\V (N) U\ * £ WjfUl, 
|/|=JV 



N 



£(i + M) n |v (n) [/|<|c/| v , ; 



A"- 



n=0 



N 



\U\ Vz , N \<E( 1 + S ) n \^(n)U\, 
n=0 

JV 

l^lv,;ivbo-E(l + ^) n |V (n) [/| 2 | So . 



n=0 



Although |s denotes restriction to So in the above formulas, we emphasize that V( n ) denotes the full spacetime 
covariant derivative operator of order n. 

Let T be an arbitrary two-form, let a, a, p, a be its null components as defined in Section\6~3[ and let E, B, (£, 23 



be its electromagnetic decompositions as defined in Section 6.8 Let A be one of the three sets T,0,Z defined in 

def 



(5.0.8a) - (5.0.8c). Let q = r — t, s"= J r + t denote the null coordinates. Then the following pointwise estimates hold: 



9.0.9a) 
9.0.9b) 
9.0.9c) 

9.0.9d) 

9.0.9e) 

9.0.9f) 
9.0.9g) 

9.0.9h) 
9.0.9i) 

9.0.9j) 

9.0.9k) 

Although 
covariant 



\T\ 2 = 2(\E\ 2 + \B\ 2 ), 



Q 



( 2 + |a| 2 + 2(p 2 + < 7 2 ), 



\ T \ 2 = \E\ 2 + \B\ 2 + \£\ 2 + \%\ 2 

= (1 + q 2 )\a\ 2 + (1 + s 2 )\a\ 2 + (2 + q 2 + s 2 )(p 2 + a 2 ), 



i 2 



|/|<JV 

H2 



I \'2> 



|/|<iV 

^(n)^\ 2 ~ |V(„)-E| 2 + |V( n )-B 



C A *T\ 2 + \i s C A T\ 2 + \i s C A *T\ 2 ), 
C A *T\ 2 + lisC^J 7 ] 2 + lisC^J 7 ^ y 



C Z ;N 



\E 



C Z ;N 



\E 



N 

\lz:N + \B%z:N + K Z ;N + \K Z ;N * EC 1 + s) 2 ^ +1) (\V {n) E\ 2 + |V (n )i?| 2 ), 

n=0 



E {{l + \q\f\\V 1 ?T\}<\T\ Cl 
\I\<N 

N 

^{(l + M)' 1+1 |V(n)^l}<+^+^;7V, 
n=0 

N 

+ T \lz;N Ieo « E (1 + r 2 )" +1 |V (n) £| 2 | So + |V (ri)J B| 2 | So . 

n=0 

denotes restriction to £q m i/ie above formulas, we emphasize that V( n ) denotes the full spacetime 



s denotes restriction to iu$ in trie above formulas, we emphasize mat V < n \ denotes trie j 
derivative operator of order n. Furthermore, in (9.0.9i), I is a translational mutli-index. 
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Let J- be any two-form. Then for any r > 0, r = \x\, we have that 



N 



(9.0.10) 



£ (1 + M)V + VlV L T \ Cot N-k-l<\ T \ Cz ;N 



Let J- be another arbitrary two-form, and let «/^-[^ r ] be the energy current vectorfield defined in (9.0.5). There 
exists an e > such that if \ T \< e, then the following estimates hold: 



(9.0.11a) 
(9.0.11b) 



Remark 9.0.5. We remark that the estimates (9.0.8a) - (9.0.10) either were proven directly by Christodoulou 



and Klainerman in [C K90j . or follow easily from the estimates of jCK90j : we nevertheless provide proofs here for 
convenience. 



Proof of Proposition 9.0.3 



Proof of (9.0.8a): The estimate (9.0.8a) can be proven inductively using (3.2.10) and Lemma 7.0.2 



Proof of (9.0.8b) - (9.0.8e): 



The estimate (9.0.8b) is easily deduced in the inertial coordinate system {a; /i } /i= o,i,2,3- 

To prove (9.0.8c), we first note that in the inertial coordinate system, the following identity holds for \i = 0, 1, 2, 3 : 



(9.0.12) 

Therefore, we have that 



(m) 



qs 



(9-0.13) \Vt m U\ = \(qs)- 1 -\\arVn M U + x l S7 g U\ < |gr|£/| v *;i- 

Since we also have the trivial identity |Vt (m) ^1 ^ |^1v^;i) ^ follows that 



(9.0.14) 



iv T( ^i<(i + Mrviv; 



Similarly, we use Lemma 7.0.1, Corollary 7.0.10, and (9.0.12) to inductively derive the following inequality: 



(9.0.15) 



|VT (Ml) -VT (Mn) c/|<(i + |ginc/|v^. 



Combining fl9 0.15[) a nd ( |9.0.8b[ ), we deduce ( I9.0.8c[ ). 



Inequality (9.0.8d) follows trivially from (7.0.2) and the Leibniz rule. Inequality (9.0.8e) then follows from 



(9.0.8d) and the fact that s - r along Sq. 



Proof of (9.0.11a) and (9.0.11b): Since j° = Q(^°\K), (|9.0.11aD is a simple consequence of Lemma 



Inequality (9.0.11b) then follows from integrating inequality (9.0.11a) over X^. 



8.4.1 



Proof of (9.0.9a) - (9.0.9f): 



If X is any vectorfield in the plane spanned by L and L, then it can be checked that 



(9.0.16) 
(9.0.17) 



\i x T\ 2 = -(XL) 2 \a\ 2 + ~(X L ) 2 \a\ 2 + \x L X k $ (a, a) + \{{X L ) 2 + {X k ) 2 }p 2 



T| 2 = \(XrJ 2 \a\ 2 + \{X L ) 2 \a\ 2 - \x L X k J (a, a) + \{{X L ) 2 + (%) 2 }a 2 , 



where ^ (a, a) =<j) K,X 9L K a \ - (<? 1 ) KX QL K a \- Therefore, 
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(9.0.18) 



+ \ixTf - l{(X L f\a\ 2 + (X L f\a\ 2 + [(X L ) 2 + (X^) 2 ](p 2 + a 2 )}. 



Now taking first X = S, and then X = T(o)) using 5x = q, Sl - -s, T(o)L = ^(o)L = _ 1> an d recalling that 
l£ = 23 = -E = it^F, B = -iT {Q) *F, we find that 



(9.0.19) 
(9.0.20) 



£| 2 + |23| 2 = l -{q 2 \a\ 2 + s 2 \a\ 2 + (q 2 + s 2 )(p 2 + a 2 )}, 
= -{|a| 2 + |a| 2 + 2(p 2 + <r 2 )}. 



|£| 2 + \B\ 2 



Adding (9.0.19) and (9.0.20), and comparing with (9.0.3), we conclude the following: 

|2 



(9.0.21) 



\T + 2 = \Ef + \Bf + |£| 2 + 1951 



We have thus proved (9.0.9a) - (9.0.9d). With the help of Corollary 5.0.3, (9.0.9e) follows similarly. 



Inequality (9.0.9f) follows from (9.0.8b), (9.0.9a), and the fact that the electric field and magnetic induction 
one-forms associated to are respectively V^E and V^B. 



Proof of (9.0.9g): 



Using the commutation identities (7.0.12a), (7.0.13a), (7.0.14a), and (7.0.14b), it follows that for any Z-multi- 
index /, we have 



(9.0.22) 


Kt (0) £^-£^(t (0) )^I< 


E E 

»=a\J\<\i\- 


\ir M c J z n 

-1 


(9.0.23) 




3 

E E 

M=0|J|<|/|- 


-i 


(9.0.24) 


\i s L z T - L z i s T\< 


3 

E E 

M=0|J|<|/|- 


\iT^x j z n 

-i 


(9.0.25) 


\isC z *T-C z i s ^\< 


3 

E E 

H=Q\J\<\I\- 


Vt m c j z t\ 

-1 



Furthermore, it follows from (9.0.9b) and (9.0.18) that for any two- form J- and each translational Killing field 
-^(aO' (A 1 = 0' 1' 2, 3), we have that 



(9.0.26) 



(9.0.22 


) - ( 


9.0.25 


) with ( 


9.0.26 



(9.0.27) 
(9.0.28) 



\i s C z T - C z i s T\ <\T\ Cz - m -i ■ 



We will now prove (9.0.9g) by induction. The base case was established in (9.0.9c). We thus inductively assume 



that (9.0.9g) holds in the case N- 1. Using (9.0.27) - (9.0.28), it follows that 
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(9.0.29) + T \\ Z , N =\ T \\ z . N _ x + Y, (Kt (0) ^| 2 + \iT (Q) C^T\ 2 + \i s L z T\ 2 + \i S C z *^\ 2 ) 

\I\=N 

^f\I z ;n-i+ E (i4iT (0) ^i 2 + i4^ (0 /^i 2 + i4^l 2 + l4^l 2 ) 

\I\=N 

=\f\1 Z ;n-i+ E (\c z e\ 2 + \c z b\ 2 ^c 1 ^ 2 + \c z *\ 2 ) 

I T\- AT 



\I\=N 

< £ (\C Z E\ 2 + |£is| 2 + \£ z <£\ 2 + 1£^| 2 ) 

|/j<iV 

= l^lL;iV + l s lL;iV + l«lL ; iV+l®lL ; iV> 



where the induction hypothesis was used to deduce the next to last line. 

For the opposite inequality, we again use (9.0.27) - (9.0.28) and the induction hypothesis to conclude that 

(9.0.30) 



\E\l z;N + \B\l z . N + \e\l z . N + \<B\l 2 



;JV 



= E (l4^ (0) ^| 2 + \C z i T{0 *T\ 2 + \C z i s T\ 2 + \£ z i S *F\ 2 ) + \E\l z . N _ x + iBll^ + l^.^ + |<B|^ 

\I\=N 

* E (h 0) 4^l 2 + KT (0) 4^l 2 + Ks4^l 2 H^4^l 2 )+ + ^+L;7v-i 

l/l=JV 



■N-l 



Proof of (9.0.9i): 



Using the facts that Vj<. .S = and Vt^^o) = f° r M = 0)1)2,3, together with the commutation identity 
(17.0. 12bh, it follows that 



(9.0.31) 
(9.0.32) 

(9.0.33) 
(9.0.34) 



IvWco)-^- *r (0) V^| = 0, 
IvWco)*-^- *t (0) V^*J c "| = 0, 

3 _ 

|V^ s ^-i s V^|< E E l»r w V^|, 

M=0|J|<]7]-1 

|vW^-^v^|<E E Kt (m) v^t|. 

M=0 1 J|<|i"|-1 



Furthermore, from definition (9.0.3), it follows that any two-form T satisfies the following inequality: 

Hv?f + \iT (0) *f\ 2 <(i + Mr 2 + t \ 2 . 



(9.0.35) 



It thus follows from (9.0.31) - (9.0.34) and (9.0.35) that 



(9.0.36) 
(9.0.37) 



|V W 7 - isV^| 2 < (1 + \q\y 2 \ T + 2 ^.| /hl , 
WW? ~ isVW S (1 + Mr 2 + ^ + 2 v ^|/|-i • 



We now prove inequality (9.0.9i) by induction, the base case being the trivial inequality \ T \ <\ T \ . We 



assume that the inequality holds in the case N-l, and we let I be any multi-index I with |/| = N. Using (9.0.31 ) 



(9.0.32 


) and ( 


9.0.36 


)- ( 


9.0.37 
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(9.0.38) 



\ V^T\ 2 = \i T(o y?F\ 2 + Kt (0) V^T| 2 + |i<?V^.F| 2 + lisV^ 2 



'(0) 



< |v^ T(0) -F| 2 + |v^t (0) ^| 2 + |v^ s ^| 2 + IvW-^l 2 + ( 1 + k\y 2 + -F+v^jv-i 

< \V^E\ 2 + \V^B\ 2 + |V^£| 2 + |V^<B| 2 + (1 + \q\y 2N \ T ■ 



We remark that we have used the induction hypothesis to arrive at the last inequality. 

To es ti mate t he \V ^E\ 2 + | V J y -B| 2 + |V^<£| 2 + |V^*B| 2 term on the right-hand side of ( |9.0.38[ ), we use ( |9.0.8a[ ), 
(9.0.8b|), (|9.0.8cD, and fl9.0.9g|) to deduce that 



(9.0.39) 



|v^| 2 + |v^£| 2 + |v^e| 2 + |v^*B| 2 < (i + Mr 2iY (|£|L ;JV + \b\I z]N + |C|^ iiV + M 2 Cz . N ) 



{i + \ q \y 2N \t\1 Z ; N . 



Combining (9.0.38) and (9.0.39), we deduce the induction step. This completes the proof of (9.0.9i). 



Proof of fl9.0.9h[ ) and fl9lT9j| ) - ( |9.0.9k[ ): 

In the inertial coordinate system ^=o,i,2,3) inequality (9.0.9h) follows from the definition (6.8.4b) of (£ and 



53, (9.0.8d), (9.0.9f), the Leibniz rule, and the fact that the coordinate functions x^ satisfy \x^\ < s, |Va^| < 1, and 



.V( 2 )^| = 

Inequality (9.0.9j ) follows from (9.0.8b), (9.0.9i), and the fact that (1 + |<7|)|.F| <\ J- \ holds for any two-form T. 
Inequality j9.0.9k[ ) then follows from fl9.0.8a| ), p.0.9f[ ), ( |9.0.9gD , ( |9.0.9hD , fl9.0.9jP , and that fact that q = s = r 
along Sq. 



Proof of (9.0.10) 



Let k > be an integer. Then since L = T( ) - 0J a T^, d = f a^'/r, we use (9.0.9i) and the fact that T^uji = 



a; a T( a )W J = 0, (j = 1,2,3), to conclude that the following inequality holds for r > : 



(9.0.40) \V k L T\<(l + \q\)- k \F\ Z ; k . 

Next, we observe the following identity, which holds for any two- form T ■ 
(9.0.41) sVl-^ ' = 2C s F + qVLT - \T . 



Now using the commutation properties of Lemma 7.0.5, we may iterate (9.0.41) to conclude that there exist 
constants C a u such that 



(9.0.42) 



(sV L ) l F = £ £ C abl {qV k ) a C b s a F. 

b=0 a=0 



Using the facts that -V lQ - Vls = 2, it follows from (9.0.42) that there exist constants C a u such that 



(9.0.43) 



6=0 a=0 



Furthermore, if I is any rotational multi-index, then using Lemma 7.0.5, the facts that V l s - Vl9 - = - 0] 
that V L q = -2, and (19.0.401), it follows that 
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(9.0.44) 



6=0 a=0 

^EEM a (i+Mr (a+fc) + ^+^ ;fc+6+ i/i 

6=0(1=0 

<(i + M)- fe + ^+£ Zife+H /|. 



Multiplying each side of (9.0.44) by + and summing over all rotational multi-indices |/| < N-k-l, we arrive 
at 09.0.101). □ 



9.1. Norms for D and B. In this section, we introduce a weighted Sobolev norm that will be used to express our 
global existence smallness condition in terms of the restriction of (D,B) to So, which we denote by (D,B). This 
norm is computed using only quantities that are inherent to So- However, during the course of our global existence 
argument, we analyze the norm -jf T(t) %c z \N, and in particular that we need the smallness of \ ^"(0) \c z \N to 
close our estimates. Since this latter quantity involves derivatives that are normal to So, we need to use the MBI 
equations to relate the size of -jf .F(O) \c z \N to the size of (D,B) as measured by the So-inherent norm. This is 
accomplished in Lemma 9.1.3| 

We now define the aforementioned weighted Sobolev norm. 

Definition 9.1.1. Let U be a tensorfield tangent to the Cauchy hypersurface So- Let d(x) denote the Riemannian 
distance from the origin in So to the point x e So; i.e., in a Euclidean coordinate system {x J }j=i,2,3 on S, 
d 2 (x) = \x\ 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . Then for any integer N > 0, and any real number 5, we define the norm of 
Uby 



(9.1.1) 



\TTV 



def 



£ L (l + d\x)f +n % n) U(x)\ 2 d 3 x. 



As discussed in detail in [CBC81], the following norm, which is used in our proof of Lemma 9.1.2, can be 



bounded by a suitable choice of one of the above weighted Sobolev norms; see Lemma A-l 



Definition 9.1.2. Let U be a tensorfield tangent to the Cauchy hypersurface So- Let d(x) denote the Riemannian 
distance from the origin in Sq to the point x e So; i.e., in a Euclidean coordinate system on So, d 2 (x) = \x\ 2 = 



(x 1 ) 2 + {x z ) z + {x 6 ) 1 . Then for any integer N > 0, and any real number 5, we define the C£ norm of U by 



„2\2 



,.3\2 



jV 



\U\\ 



2 def 

ci 



N 



£sup(l + d 2 (x))^|V (n) t/(x)| 



(9.1.2) 

The next lemma is used in the proof of the subsequent lemma, which is the main result of this section. 

Lemma 9.1.1. Let N > be an integer. Let J~(t,x) be a two-form, and let (E(t,x),B(t,x)^ be its electromagnet 
decomposition as defined in Section 6.8 Then 



ic 



(9.1.3) 



N r 

*^(0)lL;iV- E L(l + M 2 r +1 (|V(n)-5(0,y)| 2 + |V(n)5(0,y)| 2 )dV 

n=0 JR 



Remark 9.1.1. Note that on the right-hand side of (9.1.3),V( n ) denotes the full spacetime covariant derivative 
operator of order n. In particular, S/r n \E and V^B involve both tangential and normal derivatives of E and B. 



Proof. Inequality (9.1.3) follows directly from integrating (9.0.9k) over Sq. 



□ 



We now state and prove the main result of this section. 
Lemma 9.1.2. Let N > 2 be an integer. Let T be a solution to the MBI system, and let (D,B) be its electro- 



magnetic decomposition along the Cauchy hypersurface Sq as defined in Section 



6.8 



Let -jj- • \z z \N and || • ||^at be 
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the weighted integral norms defined in (9.0.6) and (9.1.1) respectively. There exists a constant e > such that if 
\\(D, B)\\ h n < e, then 



(9.1.4) 



%H0) #C Z ;N" \\(D,B)\\ 

H 



N . 
1 



Proof. We first recall that by (6.8.8a), we have that 



(9.1.5) 



E 



D + Bx(BxD) 



(1 + \B\ 2 + \D\ 2 + \D x B\ 2 ) l l 2 



Using Corollary A-4, it follows from (9.1.5) that if e is sufficiently small, then 



(9.1.6) 



:d,b) 



Iff* 



:e,b) 



We now introduce the abbreviation V d =* (B,D) and re-write equations (6.8.11a) and (6.8.11b) as 



(9.1.7) V T(0) V = £(VV) + 3(V,VV), 

where £(■) is a constant linear map (i.e., a constant-coefficient matrix), and $(V, VV), which consists of "cubic- 



order" error terms, is a smooth function V, VV. We will make use of the following consequence of Lemma A-l 



(9.1.8) 



IUII (jN-2 < ||U||j^JV. 



We now repeatedly differentiate (9.1.7) with Vj^ and V, using the resulting equations to replace Vt (0) deriva- 



tives with So -tangential derivatives, and using (9.1.8) to inductively conclude that if e is sufficiently small and 
< n < N, then 



(9.1.9) 



V(n)^ll^o<|V (n) U| + EniV (m /l 



8=1 



where the sum is taken over all non-negative integers mi,m2,---,m n such that Y]l=\ m i = n - We remark that the 
implicit constant in (9.1.9) depends on 5 and its first n derivatives with respect to V, VV, and that we have used 



( |9.1.8[ ) to bound non-differentiated factors of V in L°° from above by Ce. Multiplying each side of (9.1.9) by 



A-3 



(1 + |x| 2 )( n+1 )/ 2 = (1 + |x| 2 )( 1+ ^ti m ^/ 2 , squaring, integrating, and using Corollary 

r N r 

(9-1.10) JJi + \yJ 2 ) n+1 \V(n)V(o,y)\ 2 d 3 y_< £ JJi + \y\ 2 V +1 \v (a) v(y)\ 2 d 



it follows that 



a=o ■ 
< \\(D,B) 



We have thus shown that for sufficiently small e, we have 



(9.1.11) 



N r 

E L(l + H 2 ) n+1 (|V(„)^(0,y)| 2 + \V (n) B(0,y)\ 2 ) d 3 y M \\(D,B)f HN 
n=0 JK 1 



the > direction inequality being trivial. We remark that the left-hand side of (9.1.11) involves normal derivatives, 
while the right-hand side involves only So-tangential derivatives. 



Applying similar reasoning to equation (6.8.8a), we derive the following inequality, valid for all sufficiently small 



e : 



(9.1.12) 



f R3 (l + M 2 r +1 (|V( n) £(0,j/)| 2 + \V (n) B(0,y)\ 2 ) d 3 y < ^(1 + M 2 )" +1 (|V (n) D(0,y)| 2 + \V {n) B(0,y)\ 2 ) d 3 y. 
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(9-1.13) £ £ (l + |y|Y + H|V (n) £(0,^ 

n=0 JR 1 

The desired inequalities ( |9.1,4[ ) now follow from (9.1.3) and (9.1.13). 



□ 



9.2. An expression for V^( Jiir[J-]). In Section 12, we will bound (from above) the time derivative of the energy 



£n[F~\ defined in (9.0.7). By the divergence theorem, the analysis amounts to estimating the L 1 norms of the 
quantities V^( jjr[J~]), where j- A = CJ Z T. In this section, we take a preliminary step in this direction by providing 



an expression for V At ( Jjr[^]), where J- is a solution to the equations of variation (8.2.1a) - (8.2.1b). 
To begin, we compute that for any vectorfield X, 

(9.2.1) = (v„qw + \cr {X) ^ + h^cr - Q u nv»x v , 

where ^'itpv is defined in (3.2.12). Now according to (5.0.7c), we have that = 4tg^ u . Therefore, by (8.3.5), 

the Q* 1 ^iiv term on the right-hand side of (9.2.1) vanishes. If we also make use of equations (8.3.6) and (8.3.7), 
then we easily arrive at the following lemma. 

two-form T be a solution to the equations of variation (8.2.1a) - (8.2.1b) associated to the 



Lemma 9.2.1. Let the 

background two 



Let the two-form T be a solution to the equations of variation (8.2.1a) - (8.2.1b) associa~> 
-form T. Let = -Q^K be the energy current vectorfield defined in (|9.0.5[). Then 



(9.2.2) 



V M (j^[^]) = \H^ x F Cn Z VKX K v - K V T V7] V - {y li H^ x )T KX P u< K v + -^k" V v H^ x )f Cn P KX 



Remark 9.2.1. Because of (8.2.2a), the first term on the right-hand side of (9.2.2) is for all of the variations 
of interest in this article. 

10. The Null Condition and Geometric/ Algebraic Estimates of the Nonlinearities 

In this section, we provide a partiaf^l null decomposition of the terms appearing on the right-hand side of equa- 
tion (9.2.2 ), where T is equal to one of the iterated Lie derivatives C^J- of a solution T to the MBI system. These 
geometric/algebraic estimates form the backbone of the proof of Proposition 12.0.1, which is our main energy 
estimate. It is in this section that the null condition is revealed; as we will see, and as is described at the end of 



Section 1.3, the worst possible combinations of terms are absent from the right-hand side of (9.2.2) 



We begin with the following simple lemma, which show that covariant and Lie derivatives of null form expressions 
can also be expressed in terms of null forms. 



Lemma 10.0.2. Let m > be an integer, and let Qu\(-,-), % = 1,2, denote the null forms defined in (6.4.1a) 
(6.4.1b). Then for any vectorfield X and any pair of two-forms Q, and i = 1,2, we have that 



23 



We can prove our desired estimates without the use of a fully detailed null decomposition. 
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(10.0.1) 



a+6=m V a 1 



Furthermore, under the convention of Remark 3.2.1, for any Z -multi-index 1 ', there exist constants Ci lt j 2 such 
that 



(10.0.2) 



h+h<I 



Proof. To prove (10.0.2) in the case i = 1, we use the Leibniz rule and the fact that S7x9 = 0- m the case i = 2, we 



use similar reasoning, plus Lemma 5.0.1 



To prove (10.0.2) in the case i = 1, we use the Leibniz rule and the fact that (5.0.9b) holds for any Z e iJ. In 
ie case i = 2, we use similar reasoning, plus Corollary 5.0.3 

□ 

We now state a lemma concerning the null structure of some of the factors appearing in the terms in braces on 



the right-hand side of (9.2.2). This lemma is in the spirit of Lemma 6.4.1 



Lemma 10.0.3. If T and Q are two-forms, then 

(10.0.3) l^g^R"] < s(\t\cu\G\ + \HG\cu + \f\tt\Q\tt)- 

Proof. Lemma 



10.0.3 



follows from (6.2.4) together with the null decomposition of VK given in (6.6.3a) - (6.6.3i) 



□ 



The next lemma is a technical precursor to the subsequent one. 



Lemma 10.0.4. Let N > be an integer. Let J- be a two-form, and let H^ kX be the (J- -dependent) tensorfield 



defined in (4.3.1b). Suppose that J, J' are Z- multi-indices, and that \ J\ < N. There exists an e > such that if 
\J 7 \c z ;\N/2\ ^ e > then the following pointwise estimate holds for £ = 0, 1, 2, 3 : 

(10.0.4) 

{ci(H^ x ))v,C J zF KX < E \C J iF\{\C J inv{C J zF)\cu + |4 2 ^|^|V(4>)| + |£f ^|rr|V(£i>)|rr}- 

\Jx\+\H<\A 

Remark 10.0.2. In the above inequality, [iV/2j denotes the largest integer less than or equal to N/2, and the 



seminorms | • \cu, |V ■ \cu, etc. are defined in Definition 6.1.2 

Proof. We begin by recalling that 

(10.0.5) 



H 



1 



(mb;)^(^a + \2)t(MBI) 



Note that to avoid the possible confusion described in Remark 3.2.1 we have lowered all of the indices on T in 
preparation for Lie differentiation. We now claim that 



(10.0.6) 



\{C j z {H^ X ))v»C J zT kX <£ E \C J iT\\Q {i) {C J i^,V{C J ^))\. 

**\Ji\+\M<\J\ 



Inequality (10.0.4) then follows from (10.0.6) and Lemma 6.4.1 



To obtain (|10.0.6[), we first split H^ KX into the following four pieces: 
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(10.0.7) 




(10.0.8) 


(U)^ kX 


(10.0.9) 


(iii)^ KX 


(10.0.10) 


(iv)^ KX 



-1 \ KK ( 



-1\AA /?-2 



(^ 1 )^(5- 1 ) CC (5- i ) KK (^ i ) A ^(2)^B/)-^c^v 



-1 \KK 



1 \ AA/ 
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Since the analysis is the roughly the same for each piece, we will focus only on the term (iv). 

Differentiating the term (iv) with C z , contracting with V ^C J z T K x, and using ( 5.0.9b| ) plus Corollary 
see that a typical term that arises after expanding via the Leibniz rule is of the form 



5.0.3 



wc 



(lo.o.ii) GfWif^V^Gf 1 )^^^ 



J2 T~ \ r-7 r J ■ 

Z - 

where \J[\ + \J 2 \ = |J|. Note that the factor (g- l )^(g- 1 ) xx \C J z 2 T~^)V ^C j z T kX is equal to Q (2) (£| 3 ^, V^C J Z F). 
Now among the remaining factors / d = f (g 1 ) fM ^(g (l + ^2)^(MBl))(^(MBI)*^pZ) }' a fter fully expanding 

via the Leibniz an d chain rules, there is at most one factor of T with more than N/2 derivatives falling on it. 
Thus, by Corollary 5.0.3 and the assumptions, we have that |/| < Eijxkij'i \^z^\- Thus, in total, we have shown 
that such a typical term (10.0.11) is bounded in magnitude from above by the right-hand side of (10.0.6). This 



completes the proof for the term (iv). Terms (i) - (Hi) can be handled similarly. 



□ 



The next lemma (more precisely, its corollary) will be used to control the K j- U n^ v term on the right-hand side 



^- dcf 



of (9.2.2), where T = CzJ~, and 3 U is the inhomogeneous term from Proposition 



8.2.1 



Lemma 10.0.5. Let N > be an integer, let T be a two-form, and let H^ kX be the (J- - dependent) tensorfield 



defined in (4.3.1b). There exists an e > such that if \ J T \c z -.[N/2\ ^ e i an d I is an V Z- multi-index satisfying \I\ < N, 
then the following pointwise estimate holds: 



(10.0.12) 



< CKA v /t (4^ A ) - C z {h£ kX V ,T kX )}(C z T uC )K v 
<{(1 + S ) 2 |4^ + (1 + M) 2 |4^|} 



£ \C J z ^\{\£%nv(£ J z 3 F)\cu + l4 2 -F|a/|V(4 3 ^)| + \C% F\ TT \v(C J iF)\ TT \. 

l-WI-i 

\Jl\ + \J2\+\JzW\ 



Proof. Using (10.0.5), the definition (5.0.10) of L z , the null decomposition (6.6.2a) - (6.6.2c) for K, and Lemma 
7.0.3 it follows that 



(10.0.13) 



< (1 + sf\L l z T\cu £ £ \{{£z(h£ kX ))v»C z 't k> \ 
|J'|<]/]-i C=o" 
\J\+\J'W\ 

+ (i + E £ \{{ci(H^ x ))v,c J z F KX } 

\J'W\-i C=o 
|J|+|J'|<|/| 



Inequality (10.0.12) now follows from applying Lemma 10.0.4 to (10.0.13). 



□ 



Jared Speck 

60 

Corollary 10.0.6. Under the assumptions of Lemma 10. 0.5, we have that 



(10.0.14) 



\I\<N 

*{ E (l + S ) 2 l^lL + (l + «) 2 l4^r+l^| 2 l 

l[J|<LJV/2j J 

4 E ( 1 + *) 2 l4^lL + (1 + s) 2 \£zF\tt + (i + M) 2 I^I 2 }- 



Proof. Inequality (10.0.14) follows from (10.0.12), from the facts that 



(10.0.15) \vT\cu + |V^|rr S E l^lrw + |£z-F|rr 

ZeZ 

holds for any two- for and from simple algebraic estimates of the form ab < a 2 + b 2 . 



□ 



The following lemma wi ll be u sed to control the terms (V »H^ kX )F kX F vC R V -\(K V V„H^ kX )F^F kX appearing 
on the right-hand side of (9.2.2). 



Lemma 10.0.7. Let T,f be a pair of two-forms, and let H^ kX be the (T— dependent) tensorfield defined in 



(4.3.1b). There exists an e > such that if \T\, \\JJ- \ < e, then the following pointwise estimates hold: 



(10.0.16) 



(10.0.17) 



\{V,H^ x )T kX T\K v \ < (1 + s) 2 E {\cWcu)?f + ^zFWcu + \£>zF\tt\F\tt\ 

[J[<1 

+ (1 + M) 2 E|£^| 2 |^| 2 , 

I/I<1 

(k v v u h^ x )^ kX <(i + s ) 2 E {l^^^ + ^n^ + l^^l^r}. 

[JI<1 



Proo/. Consider the decomposition V M #£ = V m (0 KkA h-V^m)**** +V m (?^) mCkA +V^(w) KkA implied by fll0.0.7[) 



- (10.0.10). We will focus only on the case of term (z); terms (ii) - (iv) can be handled similarly. We now further 



decompose V At (i) /i ^ KA into three pieces, which up to constant factors can be expressed as follows: 



(10.0.18) 
(10.0.19) 
(10.0.20) 



V,(i"T CsA "(vAn)^ 



Using (6.2.4) and the null decomposition (|6.6.2a) - (|6.6.2c) for K, it follows that if e is sufficiently small, then 
(10.0.21) 



[V M (*0^] * + s ) 2 \F\cu\VF\ + (1 + |g|) 2 |^||V^|}|Q (1) (^,^)|, 



(10.0.22) 



(10.0.23) 



[V^ii'Y^]T KX T\K v 
[V^i!"Y^ X ]T KX T\K v 



< {(1 + s) 2 \P\ cu \F\ + (1 + M)W^|} E |Q(i)(V^,^)|, 

< {(1 + s) 2 \T\cu\T\ + (1 + \q\) 2 \T\\T\)\Q {l) {T,T% 



10.0.2 



to 



where the Q(i)(v) terms arise from the k, A indices. Also applying the null decomposition of Lemma 
the Q(i)(v) terms, the fact that (V^lvw ~ E|/|=i l^-^lvW) an d using simple algebraic estimates of the form 



ab\ < a 2 + b 2 , we conclude that each of the right-hand sides of (110.0.211) - (110.0.231) are < the right-hand side of 



Q10.0.16D . Consequently, the same is true of [V ^{iY^^T ' K \F\K 



. Note in particular that our estimates for the 
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terms in (10.0.21) - (10.0.23) involving the (1 + \q\) 2 factor are not optimal; we have simply bounded them by 
(1 + |g|) 2 Em<i \C l z T\ 2 \T\ 2 . The cases (ii) - (iv) can be handled similarly; this completes our proof of (10.0.16). 



To prove fll0.0.17D, we first consider the decomposition (iT W v H^ ikX )T^J 7 k x = [K" V v (i)& KX ]F Cri r, 



kA 



+ [K v XJ v (ii)^ KX ]T Cri T KX +[K v \J v {iii)^ KX ]T Cv T KX +[K v XJ u (iv)^ KX ]T Cri T KX implied by p0!0J| | - dl0.0.10[ ). We 
will focus only on the case of term (i); terms (ii) - (iv) can be handled similarly. We now further decompose 
[K U \7 u (i)^ KX ] = (I')^ kX +(j")Ci7«A +(jW)f»J«A ) w here 

(l')^ KXd = £f MBI) (K V Vv? Cr >)F K \ 



(MBI) 



(10.0.24) 
(10.0.25) 

(10.0.26) (/'")C^A dcf ^K v V v r^ BI) )T^T KX . 

Using the null decomposition K v Vv = |{(1 + s 2 )Vl + (1 + <? 2 )Vl}, it follows that if e is sufficiently small, then 



(10.0.27) 1(1')^^^! S {(1 + s) 2 \Q (1) (\7 L F,F)\ + (1 + \q\) 2 \Q(i)(VLT,T)\}\Q (1) (T,T)\, 

(10.0.28) i(n c ^ A ^.^ 

(10.0.29) | (/ '")C^a^- c ^ a | < ^| (1 + s) 2| Q( . )(v ^^)| + (1 + | g |)2| Q( . )(v ^ ^iji^^^ ^ } | 

i=l 



where the Q(i)(Vz,.F, F), Q(i)(V lJ 7 , F) terms in (10.0.27) arise from the n,( indices, the Q^(!F,!F) terms in 
(10.0.27) arise from the n,\ indices, the estimate (10.0.28) follows from (10.0.27) by interchanging the roles 
of and ft, A, the |Q(i)(^ r ,^")[ term in (10.0.29) arises from both the (,7] and the k,X indices, and the 
Q(i)(y , , Q(i)(V lJ" F) terms in (10.0.29) arise from the fact th at I(mb i) can De expressed as a function of 
null forms; see (4.2.2). Also applying the null decomposition of Lemma 10.0.2 to the Q(j)(v) terms, using the fact 
that (V^lvw ^ £]J|=i l^z-^lvWi and using simple algebraic estimates of the form |a£>| < a 2 + b 2 , we conclude that 
each of the right-hand sides of (10.0.27) - (10.0.29) are < the right-hand side of (10.0.17). Therefore, the same 
is true of \K V v (i)'' r,KX \. We remark that we do not need the full structure of the right-hand sides of (10.0.27) - 
(10.0.29) to conclude the desired estimates; rather, we only need the fact that the right-hand sides of (10.0.27) - 
(10.0.29) are at least quadratic in the Q(j)(-,-)- The cases (ii) - (iv) can be handled similarly. 

□ 

Finally, the last lemma in this section will be used to control the terms inside braces on the right-hand side of 



(9.2.2). 



Lemma 10.0.8. There exists an e > that if \J-\, \VJ-\ < e, then the following estimates hold: 



(10.0.30) 



(MBI) -OeA^ -r ^ J- J- + *(2) l (MBI)) r S K \{ S / £ J- J-f-j 



+ WiMBI)*^^^ - + *< 



□ 



S s{m 2 cu\F? + I^I'I^IL + I^IttI^Itt)- 
Proof. Inequality ( |10.0.30| ) follows from (6.3.5a) - ( |6.3.5d ), Lemma 6.4.1 and Lemma 10.0.3 

11. Global Sobolev Inequality 

In this section, we recall a version of the global Sobolev inequality that was proved in [CK90]. This fundamental 
inequality allows us to deduce weighted L°° bounds for a two-form T from weighted L 2 bounds of the quantities 
C^T. It provides the mechanism for deducing the factor in our estimate (12.0.1). Although many of the 
estimates in this section were proved in |CK9 0l . we reprove some of them for convenience. However, in order 
to derive the improved decay estimates ( 11.0.3a[ ) - (11.0.3b) for the a component of a solution T to the MBI 
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system, we will have to make use of the null decomposition equation (6.7.7a). The fact that inequalities (11.0.3a) 
- (11.0.3b) hold is another manifestation that the nonlinearities in the MBI system satisfy a version of the null 
condition. 



Proposition 11.0.9. Let T be a two-form, and let a, a, p, a be its null decomposition as defined in (6.3.1a) 



(6. 3. Id). Let M > 2 be an integer. Then in the interior region {(t,x) \ \x\ < 1 + t/2}, we have that 



(11.0.1) 



S/ (m) Ht,x)\<(l + S )- 5 ' 2 - m %T% C2 



;M, 



m = 0,-,M-2. 



Ln the exterior region {(t,x) [ \x\ > 1 + t/2}, we have the following estimates: 

(11.0.2a) |VM t(m) a(t,x)\ < (1 + S )- W " m (l + M)" 3 / 2 "* f J= %c z;M , < k + I + m < M - 2, 

(11.0.2b) |v|vi ? (m) (a(t,x),p(t,x),a(t,x))| < (1 + S )- 2 ^ m (l + |^| > 1/2 ^ f ^ ftc z ;M, < fc + I + m < M - 2. 
Furthermore, if T is a solution to the MBL system ( |4.3.2a ) - (4.3.2b), then we have the following improved 



estimates for a ■ 

(11.0.3a) \V l L f(m) a(t,x)\ < (1 + S )- 5/2+m | f | £z;M , 0<l + m<M-2, 

(11.0.3b) |v! + V^ (m) a(t,x)|<(l + S r 3 ^ m (l + l9ir 1/2 ^l^l^;M, 0<fc + / + m<M-3 (i/M>3). 
The following corollary follows easily by using C Z T in place of T in Proposition 



role in our proof of Proposition 12.0.1 



11.0.9 



It plays a fundamental 



Corollary 11.0.10. Let T be any two-form, and let L be any Z-multi-index such that \I\ < M — 2. Let the 

pointwise seminorms \ ■ |yyy be as in Definition 



6.1.2 



Then with r = q = r -t, s = r + 1, we have that 



(11.0.4a) 
(11.0.4b) 
(11.0.4c) 



\£ Z F\<(1 + s)-\l+\q\y 3 / 2 % T U Z ;M, 
\£ z F\cu < (1 + sY\\ + kl)- 1 / 2 | T \ Cz , M , 
\L z ?\tt £ (1 + sy 2 (l + \ q \)-^ 2 | T \ Cz , M . 



□ 



The proof of Proposition 11.0.9 is heavily based on the next lemma, which was proved in [CK90J. 



Lemma 11.0.11. [CK90, Lemmas 2.2 and 2.3] Let U(x) be a tensorfield defined on Euclidean space M 3 . Then ft 
any real number t > 1, we have that 



or 



(11.0.5) 



sup \U(x)\<(l + t)-H Y,t 2m [ 

x\<l+t/2 \m=0 J \y 



W(m)U(y)\ 2 d 3 y 



\x\<l+t/2 
def. 



<l+3t/4 



1/2 



For all iel with \x\ = \r\ > 1, we have that 



(H-0.6) \U(x)\ < r-^H I , \U(y)\ 2 + \y\ 2 \V j,U(y)\ 2 d 3 y 



\y\>r 



1/2 



where N = d r is the radial derivative. 



For all real numbers t>0 and all x e M 3 with \x\ d = r > 1 + 1/2, we have that 

(n.o.7) |tf(^|s(i + *)-Hi + M)-W f , l^(y)lL;2 + (i + ||y|-f)|v^(y)|| o;1 d 3 y) ' , 



1/2 
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def 



where N - d r is the radial derivative. 
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□ 



Before proving Proposition 11.0.9, we prove two additional technical lemmas. 



Lemma 11.0.12. [CK90, Corollary of Lemma 3.3] Let T he a two-form, and let a, a, p, a be its null decomposition 



as defined in (|6.3.1al) - (|6.3.1d|) . Then with r d = \x\, s "= r + t, q"= J r -t, the following pointwise inequality holds: 



def 



def 



(11.0.8) 

M 

£ C 1 + (1 + \q\)\V k LV l L^\Co;M-k-l + s\vlV l L a\ Co -M-k-i + s\v1V L p\c ;M-k-l + s|v|vH,Co;M-jfe-; \ <\ T \c z ;M 

k+l=0 



Proof. Inequality (11.0.8) follows from Corollary 7.0.7, (9.0.9c), and (9.0.10). 



□ 



Lemma 11.0.13. Let T be a solution to the MBI system, and let a, a, p, a be its null decomposition. Then there 
exists an e > such that if \J~\c z -\mI2\ - e i then the following pointwise inequality holds: 



M 



Af-1 



fll.O.c 



1=0 



Proof. To deduce that the first sum on the left-hand side of (11.0.9) is <\ T \c z -Mi we simply apply Lemma 



11.0.12 To estimate the second sum, we begin by recalling the equation (6.7.7b) satisfied by the null components 



of T . We write the equation relative to an arbitrary coordinate system, instead of a null frame: 

(11.0.10) 
Vlc^ - r" 1 ^- W p - t£ f u a 

- 7^MB/){(V^(i) " 2^(2)Vl^( 2 )K " 2(^" % - 2^ (2) ^ f u * (2) )cj + (Vz^(i) - 2* (2) V L ^( 2 ))a M } 

" ^(mb^(2){(VlV) - 2^(2)V^ (2 ))^a, - 2{tf % - 2* (2) ^ f u <, {2) )p - (V^ (1 ) - 2^(2)V^(2))fe} 

+ \{&L\2))tZ<*» - 2(^" % \2))P~ (V^(2))^«,} 
= 0. 

We then differentiate each side of (11.0.10) with V^,V^, and jCq, where / is a rotational multi-index satisfying 



\I\ < M - k - I - 1, multiply by r 2 r l (l + \q\) k , and use Lemma 7.0.5 , L emma 7.0.6, Corollary 7.0.18 (to exchange the 
y derivatives for r _1 - weighted rotational Lie derivatives), Lemma 10.0.2 the fact that V l t - -Vz,r - lj the fact 
that (1 + \q\) < r in the exterior region, and the assumption that \ T \c z ;[ M l 2 \ * s sufficiently small to derive the 
following inequality: 



(11.0.11) rV{l + \q\) k \vfW L a\ Co - M - k -i-i<r £ /(l + \q\f\V% via\ Co , M -k'-V-l 

0<k'<k, 
0<l'<l 

+ rr\l + \q\) k \v]y l L p\co;M-k-l +rr\l + \q\) k \V l L a\ Co -M-k-l 
+ r 2 r l (l + \q\) k ^nonlinear terms. 

After fully expanding via the Leibniz rule and using the smallness assumption on |^ r |/; z ;| m/2Jj it follows that up 
to order 1 factors, each nonlinear term on the right-hand side of (11.0.11) is of one of the following three types: 
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(11.0.12) 



(0 



Q(0(4 1 vf +1 V^^gvJ , vS^)|}|4-vi'vg(a ) (7 J p)| n I4 C V^V^|, 

4<a,6,c, 
0<fc a +ii,+|/ c | 



(11.0.13) («') = {|S(0(4 1 V^V I L 1 ? ^£§vjv^)|}|£§vjv^(«,a,p)| [I I4 c vtv^|, 

4<a,6,c, 
0<fc a +Z b +|7 c | 

(11.0.14) (m) = {l^oCrgv^V^^gvJv^lJl^gvJvJal fl 14^^, 

4<a,6,c, 
0<fc a +Z b +|J o | 

where the fc a are non-negative integers such that k\ + ■■■ + kk = k, the lb are non-negative integers such that 
l\ + ■■■ + li = I, the I c are rotational multi-indices such that |ii| H — + \lM-k-i\ < M — k — I, and the Q(j)(-, •) are null 
forms arising from the We remark that the type (i) terms arise from e.g. the derivatives of ^/^(Vl^i))^ 



on the right-hand side of (11.0.10), the type (ii) terms arise from e.g. the derivatives of f-n^Bi)^^ ^ v ^(l))* 7 ' an( ^ 
the type (Hi) terms arise from e.g. the derivatives of ^fgj)(VL ^(i))g ^- 

Each linear (in J-) term on the right-hand side of (11.0.11) is manifestly bounded by the left-hand side of 
( |11.0.8 ). Therefore, in order to prove (11.0.9), what remains to be shown is that the nonlinear terms of type 
(i) - (in) are each < r~ 2 r~ l (l + \q\ )~ k \ T \c z \M ■ 

For the type (i) terms, we use (11.0.8) and the smallness assumption on \^F\c z -[m/2\ t° deduce that 



< r -2 r -(Ji+-+J«)( 1 + | ? |)-(fei+-+fc t ) 



+ T \C Z ;M= T-\-\\ + \q\y k \ T \ Cz ;M, 



(11.0.15) 

where the r~ 2 arises from the fact that at least 2 of the factors on the right-hand side of (11.0.12 ) involve derivatives 



of the more rapidly decaying terms a,p,a; one of the factors is explicitly written, while the second arises from the 
fact that each Qu\ is a null form. The type (ii) terms can be handled similarly to the type (i) terms. In fact, we 
note that they have even better decay since they have an angular derivative *f T in place of one the Vl; we do 
not make use of this fact. 



To bound the type (in) terms, we first note that only one factor on the right-hand side of (11.0.13) involves 
the fast-decaying terms a,p, a; it arises from the Qu)- However, there is an additional power of Vl available to 
compensate. Therefore, we again use ( 11.0.8| ) to deduce that 



■ M =r- 2 r- l (l + \q\)- k \F\ CzM ■ 



|(m)| < r - 1 r-< ,1+ ~ +, ' +1 >(l + | g |)-(*i + "-**> \ T \ Ci 
We now set k = and combine the estimates for the linear terms and the type (i)-(iii) nonlinear terms, arriving 



at the estimate (11.0.9). We remark that we will use the expressions for the terms (i) - (Hi) in the case k > 1 later 



in the article, during our proof of (11.0.3b). 



□ 



Armed with the previous estimates, we are now ready for the proof of the proposition. 



11.1. Proof of Proposition 11.0.9 (global Sobolev inequalities). Most of these estimates were proved as 



Theorems 3.1 and 3.2 of [CK90J. In particular, we do not repeat the proof of (11.0.1). However, our proofs of 



(11.0.3a) - (11.0.3b) involve modifications of the arguments that take into account the special nonlinear structure 



of the MBI system. Therefore, we prove (11.0.2a) - (11.0.3b) in complete detail, and supply some additional details 
not contained in [CK90j. 



The arguments we give concern the exterior region {(t,x) \ \x\ > l + t/2}. To begin, we square inequality (11.0.8) 
and integrate over the exterior region, thereby obtaining the following inequality: 



r M or 9 

(ii-i-i) /, V (l + ||y| - tr^^ICi + - *| a )|vivL«(t, y)!^.^-*-! + ly^lviv^aCt, v)!^^-^ 



+ \y\Viy l LP(t,y)\ 2 c -M-k-i + \v\ 2 \^]yi^y)\lo-,M-k-i}d z y <! Fit) \\ z -m 
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Now for any k + l + m = 0,-,M-2, we define U(t,x) = r m+l (^l + q 2 ) m+1 V k L V l L f (m) a 
or U(t,x) d = f r l+m+1 (\fl + g 2 ) fc V^V^ f( m ) (a,^/ 3 )- If N denotes the outward normal to the sphere S rit , then the 

Corollary 7.0.18, (111. 1. lb, 



7.0.5 



fact that V^r - d r implies that V ^(^(l + q 2 ) 1 ^ 2 ^ < I. Using this estimate, Lemma 
and the fact that V ^ - ^(V l - Vl)> we arrive at the following inequality: 

(11.1-2) L , {\U(t, y)\l ;2 + + f)|v^f/(t,y)li o;1 } d 3 y <| ^(t) f ^ 



M • 



From (11.0.7) and ( |11.1.2 ), we conclude that in the exterior region, we have that 



(11.1.3) 

which implies that 



\U(t,x)\ < (1 + S)-\1 + M)- 1 / 2 I Ht) %C Z ;M, 



(11.1.4) 
(11.1.5) 



|V|V^ y (w) a(t, x)| < (1 + S )- W " m (l + M)- 3 / 2 - fc | ^(t) |^ ;M , 
|v£vi t (m) {a(t,x),p(t,x),a(t,x))\ < (1 + S )- 2 ^ m (l + \q\y l l 2 - k f ^(t) \ Lx . M . 



This proves Qll.0.2a[ ) and ( |11.0.2b[ ). 



To prove (11.0.3b), we apply the estimates ( 11.0.2a) - (11.0.2b) to the right-hand side of (11.0.11), using ( 11.0.12 ) 



(11.0.14) to estimate the nonlinear terms, and Corollary 7.0.18 to translate Lie derivative estimates into intrinsic 



(to the S r> t) covariant derivative estimates. 



To prove dll.0.3a| ), for each l + m = 0,-,M-2, we define U(t,x) d = r l+m+1 \7 l L f (m) a. Using Corollary 



7.0.18 



and 



arguing as above, we integrate the inequality (11.0.9) over the exterior region to obtain the following inequality 



r M M-1 

(n.i.6) / l|y| 2 Elyl 2i lvWi,y)li ;M-^|y| 4 £ \v\ 2l \VLvia(t,ml -, M -i-i}d 3 y<^m \\ z , 

J\y\>\+tft «■ , =0 - , =0 > 



\M ■ 



Using Corollary 7.0.18 and (11.1.6), it then follows that 



(11-1-7) f m>i+i/2 W, y)\l ;2 + |y| 2 | V„U(t, y)\lo;i) ^ ^ ^(0 

Therefore, using ( 11.0. 6[ ), we conclude that in the exterior region, we have that 



(11. li 



\U(t,x)\<r-V 2 ttHt) Wc z ;M- 



From (11.1.8), we conclude that the following inequality holds in the exterior region: 



(11.1.9) 



|Vl t(m) a(t,x)\ < 



-5/2-i-m 



This proves inequality (11.0.3a). 



12. Energy Estimates for the MBI System 

The goal of this section is to prove the most important estimate in the article: an integral inequality for the 
norm \ J~(t) \c z \N ■ The inequality, which is the conclusion of the next proposition, is the crux of the proof of 
Theorem [TJ 

Proposition 12.0.1. Let N > 3 be an integer. Assume that T is a classical solution to the MBI system 
(4.3.2a) - (4.3.2b) existing on the slab [0,T] x M 3 . Then there exist constants e' > and C > such that if 



su Pte[o,T] W \c z \N- e 'i then the following inequality holds for t e [0,T] : 
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(12.0.1) 



t m fc z ;N< C{ | HO) #L;iV + f r l ^2 + ^t) f Cz . N dr}. 



The proof of Proposition 12.0.1 will follow easily from the following lemma and its corollary. 



Lemma 12.0.2. Let N > 3 be an integer. Assume that T is a classical solution to the MBI system (4.3.2a) 



(4.3.2b) existing on the slab [0,T] x R 3 . For each \I\ < N, let jj\C l z T\ = -Q^ U K U be the energy current current 



(9.0.5) constructed out of the variation T = C^J- and the background T . Then there exists a constant e > such 
that if jf- T \ c z ;[N/2\+2^ e ) then the following pointwise estimate holds on [0,T] x R 3 : 



(12.0.2) 



|v M (j£[4^(*^)])| 



f 2 

°[A r /2j+2 



[Ht)] 



1 + s 2 



Corollary 12.0.3. Assume the hypotheses of Lemma 12. 0.2\ Then there exists a constant e > such that if 
-ff T ^c z ;[N/2\+2^ e > then the following estimate holds for t e [0,T] : 



(12.0.3) 



^ 3 V M (4[£^(t,x)]) 



d s x £ W2 l+ 2 



IHt)] 



1 + t 2 



Proof. Corollary 12.0.3 follows from integrating inequality (12.0.2) over S 



□ 



We will now prove the proposition; we will subsequently provide a proof of the lemma. 



Proof of Proposition 12.0.1 



Using the definition of £]^[J-(t)], the fact that [iV/2j + 2 < N, the divergence theorem, (9.0.11b), Corollary 
and the smallness assumption on \ T \ c z ;[N 721+2) we have that 



12.0.3 



(12.0.4) 



d 
dt 



{£ 2 N [Ht)])= E [d t (j^T(t,x)])d 3 x = £ £ 3 vJj^[4^(i,x)])d 3 x 



\I\<N- 

f 2 

c [A r /2j+2 



iTt 2 — * m w^-ir^ 



sum] 



Inequality (12.0.1) now follows from integrating inequality (12.0.4) from to t and using (9.0.11b) again. 



□ 



We now return to the proof of Lemma 12.0.2 
Proof of Lemma I12.0.2t 



First, we note that by the smallness assumption jf- T %c z ;[N/2\+2^ e > together with Corollary 
that 



11.0.10 



we have 



(12-0.5) + T \C Z ;[NI2\ <C\T % Cz ;[N/2\ + 2< Ce(l + ^(l + | | > 1 7 2 - 

The above inequality is more than sufficient to guarantee that if e is sufficiently small, then the hypotheses of 



Proposition 9.0.3 and of all of the lemmas and corollaries of Section 10 are satisfied; we will make use of these 



results in our argument below. 



By (8.2.2a) - (8.2.2b) and (9.2.2), we have that 
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(12.0.6) 



'(v^ Ka )4^A4^| + |^(v^ kA )4^4^a 

MBly -TnXy -r ^ -T J- ^ ) 



(1 + ^(^^^^^(T^ - ^*j* c ) 



Inequality (12.0.2) now follows from (12.0.6), Corollary 10.0.6 Lemma 10.0.7 Lemma 10.0.8, and Corollary 11.0.10 



As an example, we describe the estimate of the first term on the right-hand side of (12.0.6) in more detail. To 



estimate this term, we first note that by Corollary 10.0.6, we have that 



(12.0.7) 



|/|<iv 



<\ £ (l + s ) 2 \£ z ^ + (l + s ) 2 \£ z F\ 2 TT + \C J z P 2 

.\J\<[N/2\ 



x i £ {l + s f\C^\l u + {l + s) 2 \C I z ^ T + {l + \ q \f\C I z ^ 

\I\<N 



By Corollary 



11.0.10 



while by definition, 



we have that 



£|j|<Uv/ 2 j(l + s) 2 I4^IL + (1 + s) 2 \£ J zF\tt + \Lz? 



2\ r J ti2 



.J -ri2 



<(l + S )- 2 ^ /2j+2 [^], 



v|/|<JV 



(l + s ) 2 l4^l£ W + (i + 'YWzfVrr + (1 + M) a l^l 



21/./ 



\2|/./ 



the left-hand side of (12.0.7) is bounded by the right-hand side of (12.0.2). All other terms can be estimated in a 
similar fashion, and we omit the details. □ 



<\ T \\ . N . It thus follows that 



13. Local Existence and the Continuation Principle 



In this section, we briefly discuss the initial value problem for the MBI system. With the exception of the 
availability of Proposition 8.4.2, the material presented here is very standard. For the purposes of our global 



existence theorem, which is proved in Section 14 the most important fact presented is the continuation principle: 
it shows that a-priori control over the norm -jf F{t) %£ zi 3 is sufficient to deduce global existence. 



Proposition 13.0.4. Let N > 3 be an integer, and let the pair of one-forms (B, D) be initial data that are tangent 
to the Cauchy hypersurface £rj> that satisfy the constraints (6.8.9a) - (6.8.9b), and that satisfy \\(B,D)\\ h n < oo. 
Here, is the weighted Sobolev norm defined in (9.1.1). Then these data launch a unique classical solution T 
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to the MBI system existing on non-trivial a maximal spacetime slab of the form [0,T max ) x R 3 . 
the following regularity properties: 



The solution has 



(13.0.1a) 
(13.0.1b) 



^eC iV - 2 ([0,r wa:r )xR 3 ), 



k=N-2 



(B,D)tC N - 2 ([0,T max )xR 3 )n f| C k ([0,T max ),H^ k ). 



k=0 



Furthermore, either T max = oo, or one of the following two breakdown scenarios must occur: 

(1) lim 4tTmax | T{t) §c z ;N= °° 

(2) There exists a sequence (t n ,x n ) with t n < T max such that Hindoo ii M s/)[^ r (tn, % n )] = 0, 

is the norm defined in ( 9.0.6[ ). 
Remark 13.0.1. The classification of the two breakdown scenarios is known as a "continuation principle." 



where l^MBT) = ^(mbi)[^] - (l + ^(i)[^ r ] ~ ^(2)[^ r ]) ^ function of T defined in (4.2.2), and % F(t) %c z ;3 



Since Proposition 13.0.4 is rather standard, we don't provide a full proof, but instead refer to the reader to e.g. 
|H5r971 Ch. VI], |Maj84], [SS98] . |Sog95| , |Spe09b| , and |Tay97j Ch. 16] for the missing details concerning local 
existence, a nd e.g . |Spe09a| for the ideas behind the continuation principle. The crucial point is the availability of 
Proposition 8.4.2, which can be used to derive estimates for solutions to the linearized MBI system, that is, 



the equations of variation. More specifically, in constructing the local solution of Proposition 13.0.4, one typically 
uses an iteration argument or a contraction mapping argument. Both methods involve an analysis of solutions to 
the equations of variatiorj^J and in particular, they require uniform estimates of their weighted Sobolev norms; 
these uniform estimates can be derived using energy currents and the ideas contained in the proof of Proposition 



12.0.1, In particular, suitable energy estimates for solutions to the equations of variation (8.2.1a) - (8.2.1b) can 



be derived by using energy currents Ji oca ijrY^ defined by 



(13.0.2) 



v locals 



8.4.2 



Now by 

Proposition |8.4.2| if T e .ft, where .ft is a compact subset of the domain ffl of state space where the MBI theory is 



where Q u is the canonical stress from (8.3.3), and XY , is the vectorfield defined in Proposition 



defined, then we have that 



(13.0.3) 



if+ 2 <(i + /)|f| 2 <c,jL#]. 



By Remark 6.8.2 in terms of the state-space variables (B,D), this domain comprises the set finite values of (B, D). 
On the other hand, using the simple inequalities (1 + q 2 )\J- 2 \ <\ J- \ 2 and 1 + s 2 < (1 + i 2 )(l + g 2 ), together with 
the fact that \Q\ < |^"| 2 , we deduce that 

(13.0.4) j&cois^ < C-\l + S 2 )|.F| 2 < C~\l + t 2 )(l + q 2 )\T 2 \ < (1 + t 2 ) \ T \ 2 . 

Consequently, if we define the energy £i oca i-N [-^X*)] by 

(13-0.5) £local;N[Ht)] d = E £ 3 i^L^] < 



d 3 x, 



then (13.0.3) and (13.0.4) imply that 



(13.0.6) 



I Ht) #C Z ;N< Slocal AHt)] $ (1 + *) % H^) \ C z ; 



N 



We remark that the implicit constants in (13.0.6) depend on .ft. 



We now illustrate the fundamental energy estimate that can be used to deduce the desired local existence result. 



We set N = for simplicity, and consider a solution T to the MBI equations of variation (8.2.1a 



2.1b) with 



: In the equations of variation, one can think of the background T as the "previous" iterate, and T as the "next" one. 
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), (8.3.7), (9.2.1 


), (13.0.3 


), ( 


13.0.6 
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Cauchy- Sch war z inequality for integrals, it follows (as in our proof of Proposition 12.0.1) that 



(13.0.7) 

j^LaAHt)}) * m mm*--, fji + t 2 + n 2 )^,*)! 2 + a + * 2 + n 2 )i^ki^,*)i + w,x)\) d\ 

< /(£; \\F(t)\\ L ~; \\VF(t)U~){£? ocal [F(t)] +E loC rt ;N [F(t)]\\#(t,x)\ + \3(t,x)\\\ H o}, 



where is the weighted Sobolev norm defined in (9.1.1), Zx^v^ 1 ' are the inhomogeneous terms on the right- 
hand sides of (8.2.1a) - (8.2.1b), and /(.£;■) can be chosen to be a positive, increasing, continuous function of its 
arguments. We remark similar inequalities can be deduced for N > 1, and that the inhomogeneous terms would 
be measured using the norm. 



The availability of inequality (13.0.7) for solutions to the equations of variation is the fundamental reason that 
Proposition 13.0.4 holds. From (13.0.7), Gronwall's inequality, and appropriate weighted Sobolev estimates for the 
inhomogeneous terms, it can be shown that £ 2 oca i[3~{t)] can be uniformly bounded in terms of \\(B,D) 



sufficiently small. Similar inequalities can be deduced for the higher order energies £ 2 oca i-N[^(t)]- As mentioned 
above, this is the main step in deducing local existence for the nonlinear equations; the remaining details can be 
found in the aforementioned references. 

There is one additional step in the proof of local existence that we will comment on, namely the issue of showing 
that £z 2 oca ;.Ar[-^"(0)] is uniformly bounded by \\(B, D)\\ h n whenever N > 3 and T is a solution to the equations 
of variation. To accomplish this rather tedious step, one can first express the equations of variation in terms of 
(E,B) and (E,B), and inhomogeneous terms, where (E,B) and (E,B) are the electromagnetic decompositions 
of T and J- respectively. One would then use weighted Sobolev multiplication estimates, as in our proof of Lemma 



9.1.2 to deduce that 



(13.0.? 



:b,d) 



JV < oo 



1^(0) %L z ;N<f(KB,D)\\ H »), 



where / can be chosen to be a positive, increasing, continuous function of its argument. By (13.0.6), the desired 



uniform bound for Sf ocal . N [^"(0)] then follows from ( |13.0.8 ). To deduce ( 13.0.8[ ), we have assumed that both T and 



j- have the same initial data, and that £^o(^ r2 ) n+1 (l^(n)-^| 2 |£o + l^(n)-£>| 2 |so)' the relevant Sobolev-Moser norm 

for T during a proof of (13.0.8), can be bounded by a positive, increasing, continuous function of \\(B,D)\\ h n. 
In practice, during an iteration scheme, this argument would need to be slightly modified; for technical reasons, 
typical iteration schemes involve a slightly different smoothingp'|of the initial data at each stage, so that the initial 
data change slightly from iterate to iterate. 

□ 



14. Global Existence for the MBI System 
In this section, we provide a proof of our main theorem. The global existence aspect of our result will be an easy 



consequence of the energy inequality of Proposition 12.0.1 and the continuation principle of Proposition 13.0.4 
while the decay aspect will follow directly from Proposition 11.0.9, which is the global Sobolev inequality. 



Theorem 1. Let N > 3 be an integer, and let the pair of one-forms (B,D) be initial data that are tangent to the 
Cauchy hypersurface So, and that satisfy the constraints (6.8.9a) - (6.8.9b). There exists an eo > such that if 
\\(B,D)\\ h n < eo, then these data launch a unique classical solution T to the Maxwell-Born-lnfeld system (4.3.2a) 



The data are smoothed because for several reasons, one reason being that during the iteration process, it is convenient to work with 
classically differentiable functions, rather than distributions. 
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- (4.3.2b) existing on the spacetime slab [0,oo) x IR 3 . Furthermore, there exists a constant C* > such that 

%Ht) U z ;N<C4(B,D)\\ H x 



14.0.1) 



holds for all t > 0. Here, || • \\ h n is the weighted Sobolev space defined in (9.1.1), while jf- • \c z \N is the weighted 



integral norm defined in (9.0.6). In addition, the null components a,a,p,a of T , which are defined in Section 6.3 



decay according to the rates given by Proposition 11.0.9 



Proof. We will show that if \\(B, D)\\ h n is sufficiently small, then neither of the the two breakdown scenarios from 
Proposition 13.0.4 occur. To this end, let e' > be the small constant from the conclusion of Proposition 12.0.l| 
Choose a positive constant e" such that < e" < e' and such that % F(t) \c z \N^ e" = > hd\re]R3 ^a/b/)^*)?)] ^ 



1/2, where £(mbi) is defined in (4.2.2); this is possible by Sobolev embedding. Define 



(14.0.2) 



sup{i > | The solution exists on [0,t] and \ F(t) \c z \N^ e }■ 



By the local existence aspect of Proposition 



13.0.4 



we have that T max > if ||(-E>, D)\\ h n is sufficiently small. 



Applying Proposition 12.0.1, we conclude that the following inequality holds on [0,T max ) : 



(14.0.3) 



t Ht) e Cz ;N< C{ | -F(0) e Cz ;N + f\ ~7~ 2 * ^ dr\. 

^ J T=U 1 + T 3 



Applying GronwalFs inequality to (14.0.3), and using Lemma 9.1.2, we conclude that the following inequality holds 
on [0,T max ) : 



(14.0.4) 



n 7-^ )<C 2 4(B,D)\\* 

T=[J 1 + T ' 1 



Now if C*\\(B, D)||^iv < e", then the continuation principle of Proposition 



13.0.4 



imply that T max = oo. Furthermore, inequality (14.0.1) is a direct consequence of (14.0.4). 



and inequality (14.0.4) together 

□ 
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A. Appendix 



In this Appendix, we state the lemmas and corollaries that are used in Section 9.1 to relate the smallness 
condition on the inherent data (B,D) to a smallness condition on -jf ^"(0) \c z \N ■ The lemmas were essentially 
proved as Lemmas 2.4 and 2.5 of |CBC81j . while the corollaries are easy (and non-optimal) consequences of the 

lemmas; we leave their proofs as exercises for the reader. Throughout the appendix, we abbreviate d = C^(M. 3 ), 

Hf d = f (R s ), etc. Furthermore, V denotes the Levi-Civita connection corresponding to the standard Euclidean 
metric u on R 3 , and we equip M 3 with standard rectangular coordinates x. 



Lemma A-l. [CBC81, Lemma 2.4] Let N,N' > be integers, and let 5,5' be real numbers subject to the constraints 
N' < N - 3/2 and 5' < 5 + 3/2. Assume that v e Hf . Then v e Cf,' , and 

(A.l) \\v\\nN> < C\\v\\ 



C 6> 



Iff*- 



□ 
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Lemma A-2. [CBC81, slight extension of Lemma 2.5] Let Ni,—,N p > be integers, and let 5\,---,5 p be real 
numbers. Suppose that Vi e H^ 1 for i = l,---,p. Assume that the integer N satisfies < N < minjiVi, N p } and 
N < T.ti Ni-(p- 1)3/2, and that 5 < £f =1 Si + (p- 1)3/2. Then 

(A.2) ft^eiff, 



and the multiplication map 



(A.3) flf x ... x iff* -> iff, («!,-, «„) - f[vi 

is continuous. 

□ 

Corollary A-3. Let N > 2 be an integer, and let 5 > 0. Assume that v% e iff /or i = 1, and i/iat m; > are 
integers satisfying JTtj < iV. T/ien 

(A.4) (l + tf) ( ^ m ' )/2 nv (m ,^L 2 



and 



(a.5) i(i+i£i 2 ) (5+E - mi)/2 nv (mj) ^iL 2 <nii^-ii^- 



□ 



Corollary A-4. Let N > 2 be an integer, and let 5 > 0. Assume that R is a compact set, that 5 e C N (&) is a 
function, and that v\ is a function satisfying v\{E?) c .ft. Assume further that V\, V2 £ iff • fTien (i 7 o v\)v2 £ iff , 
and 

iV 

(A.6) KyowO^l^^ccJV)!!^!^ £|S 0) MM;l4- 

In the above inequality, denotes the array of all j th order partial derivatives of$ with respect to its arguments, 
and \$^\si = sup v€A \^ j \v)\. 



□ 
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